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ELECTROSTATIC POTENTIAL GRADIENTS IN A PENNING DISCHARGE 


F. Salz, R. G. Meyerand, Jr., E. C. Lary, and A. P. Walch 
Research Laboratories, United Aircraft Corporation, East Hartford, Connecticut 
(Received April 10, 1961) 


It is well known that electrostatic voltage gradi- 
ents over dimensions appreciably greater than a 
Debye length do not exist in a plasma in thermal 
equilibrium.’»? However, it has been shown that 
self-consistent potential gradients may be calcu- 
lated for appropriate nonthermal velocity distri- 
butions of ions and electrons.* The purpose of 
this Letter is to show that the Penning or oscil- 
lating-electron discharge produces such distri- 
butions together with potential gradients in the 
discharge. The mechanism by which the gradients 
are generated and maintained is the net volume 
production of ion-electron pairs. 

The existence of a potential gradient can be 
justified by the following physical argument. The 
density of particles of a given type is 


n={%*-o/a, (1) 
v 


where @ is the scalar particle flux and d the mag- 
nitude of the average velocity. If a potential gra- 
dient is assumed to exist initially (such as the 
vacuum electric field), then ions will be accel- 
erated in the direction of the electric field and 
electrons will be decelerated. Consequently, if 
the flux remains constant for both, the density 

of ions will decrease and the density of electrons 
will increase. However, if a fraction of the 
electron flux is electrostatically reflected by 

the potential gradient, then the flux of electrons 
can decrease sufficiently to maintain the charge 
density equality between ion and electrons,‘ and 
the appropriate potential gradient may exist over 
arbitrary length scales. (A mathematical treat- 





ment of this situation in one dimension is given 
in reference 3.) 

This argument may be generalized to show that 
a potential gradient not only can exist but may be 
self-generated and maintained in the presence of 
ion-electron pair production. To show that po- 
tential gradients can exist, the following physical 
argument may be made. Starting with an assumed 
potential gradient (or the vacuum electric field), 
the volume production of charged particles results 
in an increase of the total ion flux in the direction 
of the field, and a decrease in the electron flux. 
Charge density equality* requires that 


o,/8,-26,/ (2) 


so that the average velocity of the ions has to 
increase and the average velocity of the electrons 
decrease, consistent with the potential distribu- 
tion assumed. The potential distribution required 
to maintain equilibrium by satisfying the density 
equality may then be computed using Eq. (2) and 
the appropriate boundary conditions. This method 
of solution will be illustrated for a specific case. 
The electrode configuration of the Penning dis- 
charge is such that the vacuum electric field 
traps electrons and ejects ions which are produced 
in the discharge (see Fig. 1). For simplicity it 
may be assumed that: (1) ions and electrons are 
essentially constrained to motion in one dimen- 
sion (due in this example to the magnetic field); 
(2) the potential V(x) is symmetric; (3) the elec- 
trons produced in the electrostatic well are 
trapped longitudinally and lost only by lateral 
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diffusion across the magnetic field; (4) ions and electrons have mean free paths that are large com- 
pared to the length scale of the discharge; (5) electrons born at x make R(x) electrostatic reflections 
before being lost by diffusion. The equality of electron and ion density, based on Eq. (1), takes the 





form 


dp x’) 


m\N? e7L M. ee 
e = ia 
n_ +34 J [va)- ven" “laze Jy Ve)-Ve)” 


where n,,(x) is the density of electrons emitted 
from the cathode and the ion flux produced in dx 
is 


do (x) = [¢, +n, ()0, (X)IN(x)o, (x)dx. (4) 


N is the neutral gas density, o; the ionization 
cross section appropriate to the local velocity 
distribution, v,, the average speed of the cathode 
electrons, and ¢, the trapped-electron flux, given 
by 


L 
$,(*) -2f R(x')[o, &") +n, @ 0, &)] 


x N(x')o.(e"\dx’. (5) 


The factor 2 takes into account the electrons 
produced between -x and -L. A self-consistent 
potential distribution V(x) is in principle de- 
termined by the above equations but the solution 
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FIG. 1. Penning source electrode configuration and 
potential distribution. 
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do x’) 
(3) 








is difficult to obtain analytically. Recourse is, 
therefore, made to numerical solutions obtained 
by an iterative procedure using finite differences. 

The actual procedure is as follows: The region 
of interest is divided into a number of cells. 
Since the mean free path of particles is large 
compared to the dimensions of the plasma, the 
flux of trapped electrons characteristic of a given 
cell (i.e., those electrons that electrostatically 
reflect in that cell) is essentially constant through- 
out the discharge between the points of reflection. 
Moreover, the local potential determines the 
speed of these electrons and hence their density 
may be deduced. The production, the number of 
reflections, and the density of ions and electrons 
may then be computed for each cell. The produc- 
tion in the center cell is computed first. The po- 
tential required to balance the ion and electron 
charge density at the cell interface is determined. 
On the basis of the computed interface potential, 
the electrostatic attenuation of the electron flux 
is computed. This procedure is repeated for. 
each subsequent cell. The trapped-electron flux 
may then be redetermined from the computed 
production rates in each cell and the number of 
reflections each group of electrons makes by the 
evaluation of the appropriate loss rates. This 
provides the input for the next iteration which 
can then be carried out in the same manner. A 
stable, self-consistent solution is obtained when 
the results of two succeeding iterations agree 
within established tolerances. 

This method has been used to calculate the 
potential gradient in a Penning discharge. The 
resulting voltage distribution for a particular 
case (assuming only a given value of cathode 
current J, anode voltage V,, and neutral gas 
density N,, together with known collision and 
ionization cross sections) is shown in Fig. 2. 
The potential distribution measured experimen- 
tally with both emitting and conventional Lang- 
muir probes (and verified by a momentum anal- 
ysis of the ejected ions) is also plotted in Fig. 2. 








961 


ed 
ces. 
tion 


iven 


ough - 
ion. 


ty 
of 


ons 
luc- 


he 


on 





VoLUME 6, NUMBER 10 


PHYSICAL REVIEW LETTERS 


May 15, 1961 








290 





280 F INITIAL CONDITIONS 











” 

a ana 

5 @ Vo = 300 VOLTS 

~~ 969 No = 2.25% 10" MOL/CC 

2 Ie = 0.125 AMPS 

5 250 F 

~ x EXPERIMENTAL 
40 — 

& 2 © THEORETICAL RESULTS 

4 ook (CELL SIZE = 1CM) 

< 

2 220F 7 
20 F 7 
200 F 7 
e ee a 

-- oOo ! 2 3 4 5 6 


AXIAL POSITION, CM 


FIG. 2. Comparison of experimental and theoretical 
plasma potentials in the Penning discharge. 


As can be seen, there is satisfactory agreement 
between theory and experiment. 
In conclusion, it may be said that stable elec- 


trostatic potential gradients may exist in a 
plasma over a length scale, determined by the 
production mechanism, which for most cases of 
interest is large compared to the Debye length. 
These gradients are established by the nonthermal 
velocity distributions generated by the volume 
production of ion-electron pairs. The potential 
gradient may be determined, in general, bya 
simultaneous solution of Poisson’s equation and 
the appropriate equations of motion. However, 
in most cases of interest, the assumption of 
equal ion and electron densities leads to a con- 
siderable simplification and allows the determi- 
nation of the potential gradient from the flux- 
velocity quotient. 





11. Spitzer, The Physics of Fully Ionized Gases 
(Interscience Publishers, Inc., New York, 1956), p. 16. 
*L. Tonks and Irving Langmuir, Phys. Rev. 34, 876 

(1929). 

3p. J. Rose, Massachusetts Institute of Technology 
Quarterly Progress Report, No. 53, April, 1959 
(unpublished) . 

‘In general Poisson’s equation must be satisfied (i.e., 
the difference in the densities of ions and electrons is 
the source of the electric field), but for most plasmas 
of laboratory interest the approximation of equal ion and 
electron densities is sufficiently accurate and leads to 
a considerable simplification. 








EXPERIMENTS ON THE ENERGY BALANCE AND CONFINEMENT OF A MAGNETIZED PLASMA 
J. Bergstrom, S. Holmberg, and B. Lehnert 


Royal Institute of Technology, Stockholm, Sweden 
(Received April 20, 1961) 


The possibility to confine a plasma in the mag- 
netic field of a current loop has been discussed 
in a number of earlier reports.'~* It has been 
suggested’ to supply energy and to heat the plas- 
ma in crossed electric and magnetic fields by a 
method used earlier in rotating plasma devices.*»* 
In a mirror machine with a rotating plasma the 
confinement is improved by the centrifugal 
force.*;® The large radial extensions of the cur- 
rent loop configuration are of special advantage 
when a strong centrifugal confinement has to be 
realized. This is indicated both by a theory on 
particle motion® and by a macroscopic approach 
including the energy balance.‘ 

The purpose of the present experiments is to 


study how the energy balance and the confine- 
ment of a rotating plasma are influenced by par- 
ticle losses along the magnetic field lines to a 
nonconducting wall. Figure 1 demonstrates the 
experimental arrangement which consists of a 
cylindrical vacuum vessel with an interior, ring- 
shaped magnetic coil (main coil) and two external 
auxiliary coils. The lines of the magnetic field 
B are sketched in the figure. A discharge is 
produced by a transverse electric field E aris- 
ing from the voltage applied between the shield 
surrounding the main coil and the walls of the 
vacuum chamber. The energy is supplied from 

a condenser bank by means of a timed ignitron 
switch. A second ignitron can be used to study 
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back- pulses by applying a short circuit (“crow- 
bar”) across the device. 

As shown in the figure, a wall consisting of 
a movable, ring-shaped plate (a) can be placed 
across the field of the magnetic bottle at dif- 
ferent distances from the equatorial plane. The 
mean field strength in this plane is about B=0.44 
weber/m’. When the plate is in its upper position 
(¢ =165 mm) the mean values of the radial ratio 
r9/Yw and the mirror ratio B,,/By become about 
3.4 and 7.6, respectively. Consequently, by a 
variation of the plate position conditions can be 
obtained which closely correspond to the homo- 
polar ,® Ixion,® and current loop configurations 
within the same experimental device. 

The experimental results can be summarized 
as follows: 

1. Voltage and current characteristics. When 
the plate is in its upper position [Fig. 2(a)], the 
discharge starts with a large current pulse 
during which ionization and acceleration mainly 
take place and the voltage reaches a maximum 
about 0.2 millisecond from the start. Then, a 
period follows during which the voltage is de- 
caying slowly and the current remains at a low 
level. A short circuit of the device at 0.5 mil- 
lisecond from the start shows a large back- 
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pulse, indicating that a considerable amount of 
energy is still present in the plasma at this time. 
The back-pulses increase rapidly with the 
strength of the applied magnetic field. The ki- 
netic energy of the rotating plasma was found 

to decay from about 135 joules at 0.1 milli- 
second from the start to about 15 joules at 0.7 
millisecond later. The equivalent capacity of 
the rotating plasma seems to be nearly constant 
at the same time. 

With the plate in lower position [Fig. 2(b)] the 
acceleration period requires a somewhat larger 
current which slowly returns to zero at the same 
time as the energy of the plasma is dissipated. 
At 0.5 millisecond from the start, Fig. 2(b) 
shows no trace of any remaining energy in the 
form of a back-pulse; only a discharge of low 
ionization degree is still burning at this time. 

The voltage level of the decay period is much 
lower in Fig. 2(b) than in Fig. 2(a) and reaches 
a constant value after 0.2 millisecond which 
agrees well with the values observed earlier 
by Fahleson.” The same levels are also ob- 
served in a low-current discharge with the plate 
in its upper position. 

If the series resistor Rs is reduced to very 
low values, the discharge breaks down across 
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I (a) Plate in upper position. 
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(b) Plate in lower position. 


FIG. 2. Discharge current/ and voltage V as func- 
tions of time ¢ for hydrogen at an initial density cor- 
responding to 45 X10~° mm Hg. The condenser bank, 
which has a capacity of 330 microfarads and an initial 
voltage of 4.5 kilovolts, supplies the energy to the dis- 
charge through a series resistor, R, =0.4 ohm. 

(a) Plate in upper position. Back-pulse of current 
is obtained at 0.5 millisecond from start. (b) Plate 
in lower position. No back-pulse is observed. 


the plate when it is in the lower position and 
the current starts to oscillate violently. How- 
7 ever, with the plate in the upper position the 
voltage still remains at a high level and the cur- 
rent drops to a low value after the acceleration 
period has been finished. 

When the plate was placed in an intermediate 
position, electrical characteristics were ob- 





tained about halfway between those correspond- 
ing to the upper and lower positions. 

2. Light intensity. Photometric measure- 
ments on a hydrogen discharge indicate a strong 
pulse of light during the early stages of the dis- 
charge when the plate is in its lower position. 
When the plate is in the upper position, the 
light intensity decreases considerably and its 
duration is extended at the same time. 

3. Gas production and wall constitution. For 
a nitrogen discharge with a density correspond- 
ing to 9x10-* mm Hg at the start, the relative 
gas production was found to decrease gradually 
from 190% to about 3% when the plate was 
moved from the lower to the upper position. 

After having been exposed to more than 1000 
discharges the insulator surfaces (b) and (c) of 
Fig. 1 showed very small signs of contamination 
but plate (a) in the same figure was heavily dam- 
aged by the discharge. 

The present results indicate that the inter- 
action between a rotating plasma and the sur- 
rounding walls of its container can be reduced 
considerably by using a magnetic field with 
large radial ratios and mirror ratios. Further, 
there are apparently no severe losses of energy 
due to instabilities since the plasma contains a 
considerable amount of energy for more than 
0.7 millisecond. 

A detailed report is under preparation. 





'B. Lehnert, Nature 181, 331 (1958). 

2B. Lehnert, J. Nuclear Energy, Sec. Cl, 40 (1959). 

*B. Bonnevier and B. Lehnert, Arkiv Fysik 16, 
231 (1960). 

‘B. Lehnert, Arkiv Fysik 18, 251 (1960); Revs. 
Modern Phys. 32, 1012 (1960). 

50. A. Anderson, W. R. Baker, A. Bratenahl, 
H. P. Furth, H. P. Ise, W. B. Kunkel, andJ. M. 
Stone, Proceedings of the Second International Con- 
ference on the Peaceful Uses of Atomic Energy, Geneva, 
1958 (United Nations, Geneva, 1958), Vol. 32, p. 155. 

°k. Boyer, J. E. Hammel, C, L. Longmire, D. 
Nagle, F. L. Ribe, and W. B. Riesenfeld, Proceedings 
of the Second International Conference on the Peaceful 
Uses of Atomic Energy, Geneva, 1958 (United Nations, 
Geneva, 1958), Vol. 31, p. 319. 

"U, V. Fahleson, Phys. Fluids 4, 123 (1961); see 
also H. Alfvén, Revs. Modern Phys. 32, 710 (1960). 
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MICROWAVE MODULATION OF THE ELECTRO-OPTIC EFFECT IN KH,PO, 





I. P. Kaminow 
Bell Telephone Laboratories, Holmdel, New Jersey 
(Received April 17, 1961) 


KH,PO, (KDP) is a colorless transparent crys- 
tal that belongs to the tetragonal group 42m (V4), 
which lacks a center of inversion, and exhibits 
a linear electro-optic effect. It undergoes a 
ferroelectric phase transition at 120°K.' Accord- 
ing to a model of KDP and isomorphous com- 
pounds, the dielectric polarization depends upon 
the motion of the protons in a double minimum 
potential.',? These protonic motions can alter 
the symmetry and magnitude of the electronic 
polarizability and, hence, the optical properties 
of the crystal. 

In the paraelectric phase, KDP is optically uni- 
axial with the optic axis along the tetragonal Z 
axis. Light propagating along Z travels with the 
same velocity irrespective of the plane of polar- 
ization. The application of an electric field E 
along Z reduces the crystal symmetry to ortho- 
rhombic and the optical symmetry to biaxial.® 
Hence, the index of refraction for light prop- 
agating along Z is different for light polarized 
in the two principal planes. After passing through 
a length L, the phase retardation [ between the 
principal waves is 


[=7EL/V,, (1) 


where V, is the (zero-strain) value of EL for 
half-wave retardation. 

Measurements of the dielectric constant at a 
frequency of 25 kMc/sec (at 25°C) agree with the 
low-frequency dielectric constant of a clamped 
crystal.'»* Since the electro-optic effect is pro- 
portional to the polarization,’ one would not ex- 
pect any appreciable reduction in V, below 25 
kMc/sec.® Previous measurements have shown 
V, to be frequency independent to 500 Mc/sec.® 
We have observed the electro-optic effect in 
KDP at room temperature and 9.25 kMc/sec 
and find V, to be of the same order as the low- 
frequency value. 

The apparatus | Fig. 1(a)] consists of a 25-watt 
zirconium concentrated arc lamp with a Corning 
7-69 infrared filter peaked at 8000 A, a Glan- 
Thompson polarizer and crossed analyzer, an 
RCA 7102 (S-1) photomultiplier, and a pulsed 
X-band magnetron which feeds a cylindrical 
cavity containing the KDP rod. The cavity is 
filled with polystyrene and the dimensions are 
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FIG, 1. 
apparatus. 


(a) Transmission apparatus. (b) Reflection 


adjusted so that the microwave phase velocity 
approximates the light velocity when the cavity 
is excited in the “TM,,,-like” mode. If the cavity 
standing wave is decomposed into a forward and 
a backward travelling wave, it can be shown that, 
to good approximation, the backward wave pro- 
duces no net retardation when the velocities are 
matched, while the forward wave provides a 
constant E at a given point in the light wave 
train as it passes through the crystal.’ The 
light velocity in the crystal is c/1.47; the micro- 
wave velocity at 9.25 kMc/sec is c/1.38. 

With the polarizer and analyzer crossed, pulses 
from the photomultiplier are observed coincident 
with the microwave pulses (3 usec at 60 pulses/sec). 
In order to obtain direct evidence of microwave 
modulation, a movable mirror and half mirror 
are placed as shown in Fig. 1(b). When the trans- 
it time for light leaving the end of the crystal, 
reflecting from the movable mirror, and re-enter- 
ing the crystal is equal to an odd multiple of the 
microwave half-period, light pulses are observed 
in the photomultiplier. When the mirror is moved 
a distance D equal to one-quarter the free-space 
microwave wavelength, no pulses are observed 
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FIG. 2. Photomultiplier response with mirror set 
for maximum signal (top). With mirror displaced A/4, 
the response is nearly coincident with that obtained by 
covering the source (bottom). (3 wsec/div. Circuit 
time constant ~ 4 psec.) 


(Fig. 2). In the first case, the retardation is 
doubled, while in the latter the net retardation 
is zero. . 

The time-average relative intensity for mod- 
ulated light with crossed polarizers, p,, and 
with parallel polarizers, py, is® 


p,=all-Jo(t)], py=all+do(r,)), 2) 


where J, is the zero-order Bessel function and 
ly is the peak retardation. With a peak power 
P of 760 watts absorbed in the cavity, p, =p i 
(Fig. 3) and Py = 2.40. 

The peak field intensity in the forward wave 
E may be estimated by assuming a uniform 
longitudinal E across the rod with sinusoidal 
variation along the rod and assuming all the 
power P to be dissipated in the rod. Then 


E ~(QP/wee,v) ~70 volts/mm, (3) 


in which € and @ are the dielectric constant and 
dielectric Q for the rod at angular frequency w 





FIG. 3. Photomultiplier response with parallel 
polarizers (upper) and with crossed polarizers (lower) 
for p= py 


and v is the rod volume (L =35.5 mm, diameter 
4 mm). In the absence of 9-kMc/sec data, we 
have used the 25-kMc/sec values* « ~20, Q ~30. 
For ry = 2.40, 

Vo" sEaL/r,~7 kv. (4) 
This estimate may be compared with the low- 
frequency clamped value of 11 kv obtained by 
extrapolating measurements of W. L. Bond to 
8000 A. 

KDP is transparent between 4000 A and 
13000 A. Since V, decreases with decreasing 
wavelength,® the effect may be enhanced by op- 
erating at shorter wavelengths. Further en- 
hancement, with increased loss, may be ob- 
tained by increasing L and by operating nearer 
to the Curie temperature, either by reducing 
the temperature, or by choosing an isomorphous 
compound with a higher Curie temperature.’ 

The tunnelling frequency’ for protons in the 
double minimum potential occurs at ~200 cm™ 
(6000 kMc/sec). Therefore, it may be possible 
to extend the modulating frequency to much 
higher frequencies than the 9.25 kMc/sec used 
here. 

It may also be possible to construct a broad- 
band low-power microwave light modulator, 
using more sophisticated structures or other 
materials, for use in a practical optical com- 
munication system.® The bandwidth for the 
present structure is limited by the cavity Q: 
Af~f/150~60 Mc/sec. In a pure travelling wave 
structure, a limitation would be the band over 
which the microwave and light velocities can be 
matched. 

A microwave light modulator might also be 
used as an instrument for studying short relax- 
ation times and other effects in optically pumped 
systems or for refined measurements of the 
velocity of light. 

It is a pleasure to acknowledge instructive dis- 
cussions with R. Kompfner, A. G. Fox, E. H. 
Turner, and R. F. Trambarulo, and the assist- 
ance of F. A. Dunn. 





'W. Kinzig, in Solid State Physics, edited by F. Seitz 
and D. Turnbull (Academic Press, Inc., New York, 
1957), Vol. 4. 

*R. Blinc, J. Phys. Chem. Solids 13, 204 (1960). 

3B. H. Billings, J. Opt. Soc. Am. 39, 797 (1949). 

‘w. A. Yager, quoted in Piezoelectric Crystals and 
Their Application to Ultrasonics, edited by W. P. 
Mason (D. Van Nostrand and Company, New York, 
1950), p. 258. 
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5N. Bloembergen, P. S. Pershan, and L. R. Wil- 
cox, Phys. Rev. 120, 2014 (1960). 

®R. O’B. Carpenter, J. Opt. Soc. Am. 40, 225 
(1950). 

"The use of a travelling wave structure has been 
suggested by N. Bloembergen (talk, NEREM, Boston, 
No. 1960] and the possibility of using the standing wave 
structure by R. Kompfner. 

8clark, Holshouser, and von Foerster, Technical 


Note 1-2, Engineering Experimental Station, Univer- 
sity of Illinois, 1957 (unpublished). 

°E. H. Turner attempted to observe microwave light 
modulation (1952) using NH,H,PO, (ADP) in a rectangu- 
lar TE, cavity but was unsuccessful because of the 
short length of ADP permitted by the structure. Bloem- 
bergen, Pershan, and Wilcox® have suggested, inde- 
pendently, the use of dihydrogen phosphates for micro- 
wave modulation of light. 





EVIDENCE FOR FOCUSING COLLISIONS IN IRRADIATED PLATINUM 


E. Ruedl, P. Delavignette, and S. Amelinckx 
Solid State Physics Department, Centre D’ Etude de L’ Energie Nucléaire, Mol, Belgium 
(Received April 21, 1961) 


Focusing collisions along close-packed di- 
rections in face-centered cubic metals were 
first predicted by Silsbee. A more elaborate 
theory has recently been developed by Leibfried 
and co-workers,’ for the case of (110) focusing. 

Gibson, Goland, Milgram, and Vineyard® found 
evidence for this process by means of machine 
calculations on a dynamical model. They 
showed, moreover, that focusing also occurs 
along (100) and (111) directions; this was at- 
tributed to the influence of the neighboring rows. 

Direct experimental evidence is practically 
limited to the sputtering experiments of Thomp- 
son.* He found that preferential ejection of at- 
oms occurs along (110), (100), and (111) direc- 
tions. The focusing effect along the non-close- 
packed rows was explained using an optical an- 
alogy.°® 

The rapid decrease in internal friction on ir- 
radiation has been attributed to the preferential 
formation of defects at dislocations as a conse- 
quence of focusing collisions.® 

In this note we wish to present some rather 
direct evidence of a completely different type, 
which was obtained during the electron optical 
examination of irradiated platinum foils. 

The specimens consisted of about 0.1-micron 
thick beaten foil, annealed for one hour at 800°C 
in order to eliminate part of the dislocations and 
to obtain a reasonable grain size. The anneal 
also resulted in the formation of annealing twins, 
generally limited by coherent twin boundaries. 

Such specimens were mounted in contact with 
uranium foils and irradiated at reactor temp- 
erature (= 80°C) in the BR-1 reactor over peri- 
ods of $, 2, and 5 hours. The number of fission 
fragments that penetrated into the platinum foil 
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under these circumstances were, respectively, 
2.5x10°, 1x10*°, and 2.5x10'° cm~?. After the 
irradiation the specimens were examined with 
the Philips electron microscope operated at 
100 kv. Figure 1 shows the increasing concen- 
tration of black dots with increasing dose. Such 
“measles” have been observed previously in 
other irradiated materials.’~*® 

By pulse annealing it was verified that the 
“measles” disappear at about 400-500°C in ac- 
cord with the annealing peak found by Piercy’® 
and attributed by him to the motion of vacancies. 
From this evidence and from the approximate 
proportionality of their concentration with dose, 
it is reasonable to accept that the black dots cor- 
respond to defect clusters. The most significant 
feature of this investigation is, however, that 
measles appear preferentially along coherent 
twin boundaries as shown in Fig. 2; and much 








less pronounced along ordinary boundaries. 

This fact can, of course, not be attributed to 
preferential impinging of fission fragments along 
the coherent twin interface. 

It is therefore difficult to avoid the conclusion 
that momentum transfer has taken place result- 
ing in the preferential creation of defects along 
the twin boundary. Since coherent twin bound- 
aries do not act as sources or sinks for point 
defects, as shown by Barnes,’ the defects would 
not easily disappear there as would be the case 
if they were formed at ordinary boundaries. 

The preferential formation of defects at coher- 
ent twin boundaries can easily be understood on 
the basis of the theory of focusing collisions. 
Suppose that momentum transfer occurs along 
a (110) direction, as would be the case in the 
Silsbee type of focusing. The (110) direction 
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FIG. 1. Increasing concentration of defects in plati- 
num irradiated in contact with uranium foil. The doses 
were, respectively, (a) 2.5 <10" n/em?, (b) 1x10" 
n/em*, and (c) 2.5 <10"* n/cm?. 


in one of the crystals of the twin becomes a non- 
focusing (114) direction in the other. Coherent 
twin boundaries are in fact effective barriers 
for focused collision chains, without being sinks 
for the defects. 

Primary collisions occurring in a strip of 
material on both sides of the twin boundary and 
causing focusing along (110) rows, intersecting 
the twin interface, will produce defects on pas- 
sing the latter. The width of this strip is deter- 
mined, either by a purely geometrical factor, 
namely, the possibility of reaching the twin in- 
terface along (110) directions, or by the range 
of the focusing chains. This width is clearly a 
function of the foil thickness, at least for thin 





FIG. 2. Preferential formation of defects along 
coherent twin boundaries (Tw) in irradiated platinum. 


foils. In very thin foils one would not expect to 
find a detectable effect since the strip becomes 
too narrow. This may explain why, e.g., Gren- 
all® did not find any such effect in very thin gold 
particles, although they also contained twins and 
stacking faults. 

It is clear that similar considerations apply to 
other focusing directions. 

A more detailed account concerning the elec- 
tron microscopic observation of dislocations 
and radiation effects in platinum is being pub- 
lished elsewhere. 

This work is part of a research program un- 
dertaken at Centre d’Etude de l’Energie Nu- 
cléaire, Mol, Belgium, on behalf of the joint 
research and development program of the 
agreement for cooperation between Euratom and 
the United States. 

We would like to thank Mr. J. Finet for skill- 
ful preparation of specimens and Mr. H. Beyens 
for careful photographic work. 
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SUPERCONDUCTIVITY OF Nb,Sn IN A PULSED MAGNETIC FIELD* 


J. O. Betterton, Jr., R. W. Boom, G. D. Kneip, and R. E. Worsham 
Oak Ridge National Laboratory,f Oak Ridge, Tennessee 


and 


C. E. Roos 


Vanderbilt University, Nashville, Tennessee 
(Received April 24, 1961). 


Kunzler, Buehler, Hsu, and Wernick’ have 
shown that quenching of superconductivity ina 
niobium-clad wire with a Nb,Sn core occurs at 
magnetic field strengths greater than 88 kilogauss 
and at core current densities in excess of 10° 
amp/cm?. Samples of this wire were prepared by 
methods suggested in reference 1. Critical cur- 
rent as a function of field was measured by a 
pulsed-field, pulsed-current method at 4.2°K in 
the Vanderbilt University 200-kilogauss magnet. 
The transition time from superconducting to 
normal state and the resistance were also de- 
termined. 

A longitudinal specimen (0.38-mm diam and 
2.5 cm long) and a constantan standard resistor 
were mounted at the end of two pieces of 3-mm 
diam, thin-wall Monel tubing as shown in Fig. 1. 
The Monel tubes served as the outer shields of 
coaxial lines which carried the sample voltage 
and current signals. A pickup coil to measure 
the field strength was placed adjacent to the 
test samples. The coils of the Vanderbilt mag- 
net are arranged in a Helmholtz configuration to 
provide a field uniform to 2% over an experi- 
mental volume of 30 cm*. Current pulses of 30 
to 150 microseconds were used at the peak of 
the magnetic field (half-period 1.2 milliseconds). 
The magnetic field was constant to within 2% 
during the longest current pulses. By mounting 
the sample as noninductively as practicable, a 
resistive sensitivity of 10 yohm was achieved, 
as checked with copper samples at a series of 
temperatures. A transverse specimen was 
mounted in a similar way. 

Typical oscilloscope traces of the measured 
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current and voltage are shown in Fig. 2. In the 
upper set, the voltage across the sample illus- 
trates a purely inductive response for a nonre- 
sistive superconductor. In the lower set, the 
voltage across the sample increases abruptly 
when the sample becomes resistive. The critical 
current is determined by measuring the current 

at the point where a resistive component is first 
detected. After the inductive voltage is subtracted, 
the normal-state resistance is approximately 
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FIG, 2. Upper set: superconducting sample, Ldi/dt 
response. Lower set: superconducting-to-normal 
transition, L(di/dt)+iR response. Top traces: sample 
current, 89.5 amp/cm, positive direction downward. 
Bottom traces: sample voltage, 100 mv/cm, positive 
direction upward. 


one milliohm. Both the transition time and the 
normal resistance depend on the previous history 
of the sample. Resistances as small as 10 pohm 
and transition times as short as 2 usec were ob- 
served. These effects are being investigated and 
will be reported in a later communication. 

The critical current was plotted as a function of 
field for two samples; see Fig. 3. Sample No. 1 
was measured both parallel and transverse to the 
field, while sample No. 2 was oriented parallel 
to the field only. The curves show that all sam- 
ples carry currents of 210 to 260 amp in a zero 
external field. Sample No. 1 carries 100 amp 
at 100 kilogauss with the field parallel to the wire 
axis. A much smaller critical current occurs 
with the field oriented in a transverse direction. 
Both the current-carrying capacity and the max- 
imum field at which current is carried are lower 
in sample No. 2. The difference between speci- 
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FIG. 3. Critical current vs applied magnetic field 
for Nb-Nb;Sn wires at 4.2K. (Wire diameter, 0.038 
cm; core diameter, 0.015 cm.) 


mens, shown in Fig. 3, is not entirely unexpected; 
the properties of hard superconductors are known 
to vary with slight differences of strain and de- 
fect structure. 

Domain structure would appear to explain the 
differences obtained for parallel and transverse 
orientations. Only a small volume of the Nb,Sn 
is superconducting at higher fields according to 
the magnetization curve of Bozorth, Williams, 
and Davis.” Studies of the intermediate state of 
various superconductors by Schawlow,® DeSorbo,* 
and Haenssler and Rinderer® suggest that the 
internal structure tends to become arranged in 
domains with interfaces parallel to the field and 
normal to the current. A domain structure of 
this type in Nb-Nb,Sn could account for the much 
smaller critical currents for transverse orienta- 
tion. This can explain the differences found be- 
tween the transverse and parallel orientations 
in soft superconductors. In a hard superconduc- 
tor such as Nb,Sn, however, the superconducting 
regions should be determined primarily by the 
macroscopic physical structure of the specimen. 
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Since this structure is oriented along the axis of 
the wire, it is somewhat surprising that such 
large differences between the parallel and trans- 
verse orientations were found. 

The critical current at zero external field was 
also measured with direct currents in the sam- 
ples. The direct-current values were 1/3 lower 
than those obtained with pulsed techniques. This 
is attributed to heating at the current joints or 
in the cladding which surrounds the Nb,Sn core. 
This heating is substantially reduced with pulse 
techniques. 





*Work supported in part by the National Science 


Foundation and T.E.A.S.P. 

tOperated for the U. S. Atomic Energy Commission 
by the Union Carbide Corporation. 
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FERMI SURFACE TOPOLOGY OF EVEN-VALENT METALS 
FROM THEIR MAGNETORESISTANCE ANISOTROPY 


E. Fawcett 
Physics Department, Royal Radar Establishment, Malvern, England 
(Received March 27, 1961) 


The anisotropy of the magnetoresistance of Au 
and Cu has recently been successfully interpreted 
by employing the concept of open orbits on a mul- 
tiply-connected Fermi surface in a periodically 
extended zone scheme.’ Open orbits occur for a 
limited range of field directions, and the trans- 
verse magnetoresistance Pyq in a sufficiently 
high field H with the current making an angle a 
with the direction of the open orbits is then of 
the form,’ 


‘on —— 
Phas A +BH? cos’a. (1) 


A plot of the resistance in a fixed high field whose 


direction varies is characterized by sharp maxima 


for field directions permitting some open orbits 
(with a #27/2), which rise above the low resistance 
obtained for most field directions for which only 
closed orbits occur and the resistance saturates. 
This behavior is in marked contrast with the 

anisotropy of the magnetoresistance observed for 
Mg,° Zn,* Cd,* Pb,® and Sn° (the transition metals 
are not considered here and insufficient data are 


available for other even-valent metals). For these 


metals the resistance varies quadratically with 
field and is uniformly high for most field direc- 
tions, falling to sharp minima for certain field 
directions where it varies less rapidly. 

The proposed explanation of this behavior in the 
even-valent metals is as follows. Consider the 
range of field directions for which closed cyclo- 
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tron orbits occur over the whole Fermi surface, 
and fields so large that for all such orbits w 7 
>» 1, w, being the cyclotron frequency and 7 the 
relaxation time averaged over the orbit. If the 
component o,, of the conductivity tensor with 
field H along the z axis is expanded in ascending 
powers of 1/H, the leading term is of the form,’ 


_ e@2 
Cy HRD J S(0,)ab, (2) 


Here the area S;(p,) enclosed by the intersections 
with the sheet i of the Fermi surface by the plane 
b, =constant will, in general, be the sum of pos- 
itive terms for closed orbits containing occupied 
electron states of lower energy within and imme- 
diately adjacent to each orbit, and of negative 
terms for orbits containing hole states. For the 
case where all sheets of the Fermi surface are 
simply connected, the range of integration need 
be only over a single Brillouin zone (with suita- 
ble remapping), and Eq. (2) results in the well- 
known expression,?~* 


e 
i" y n, “nN (3) 
Here N e is the number of atoms per unit volume, 
and mg and np, are the numbers of electron and 
hole states per atom, which are well defined 
since each sheet consists of closed surfaces con- 
taining either electrons or holes. 
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But Eq. (3) also applies to the case where one 
or more sheets of the Fermi surface are multiply 
connected, provided no open orbits occur in the 
range of field directions considered and w,7>1 
for all closed orbits. This may be seen by ex- 
tending the range of integration to sufficient peri- 
odically extended zones to include the most ex- 
tended orbits of any practical significance. The 
value of the integral for each sheet in Eq. (2) is 
still independent of the field direction: Its sign 
determines the electron or hole character of the 
sheet, and its numerical value is accordingly the 
electron volume ngNgh* /2 or the hole volume 
nj,N gh’ /2 enclosed by the sheet. This is true 
despite the fact that on the same multiply-con- 
nected sheet both electron and hole orbits may 
occur, since in the case of an electron sheet 
every hole orbit will be surrounded by an elec- 
tron orbit and vice versa for a hole sheet.’ This 
provides a useful criterion for determining the 
electron or hole character of a multiply-connec- 
ted sheet. 

Moreover, inspection of the proposed Fermi 
surface topologies for the metals considered here 
according to this criterion shows that mp and np 
are equal (and it seems a reasonable hypothesis 
to suppose this to be true for any even-valent 


metal). It follows that the leading term in oy, 
vanishes for all field directions for which only 
closed orbits with w,7 >1 occur, and the com- 


ponents p,, and py, of the resistivity tensor then 
rise quadratically with the field.® 

The transverse magnetoresistance in a fixed 
high field will therefore be uniformly high ex- 
cept for minima in field directions where there 
are open orbits making an angle a ~7/2 with the 
current direction, or where there are closed 
orbits for which w,7<1. In the first case the 
field dependence of the resistance is given by 
Eq. (1), saturating at the highest fields when 
a=7/2, while in the second case the resistance 
will saturate at intermediate fields but again rise 
quadratically at fields sufficiently high to make 
wT > 1 for all orbits. 

The multiply-connected hole surface for Mg is 
shown in Fig. 1; another hole surface is split off 
by the spin-orbit coupling,® but this need not be 
drawn since it is simply connected, as are the 
electron surfaces. Measurements of the de Haas- 
van Alphen effect in Mg and Zn are consistent 
with this picture of the Fermi surface, while sim- 
ilar measurements for Cd suggest that the arms 
of the “ring” parallel to the basal plane are 
pinched off.'° 











FIG. 1. The multiply-connected Fermi surface con- 
taining holes in the first Brillouin zone for Mg® (and 
probably Zn’), drawn in a periodically extended zone 
scheme. aa—open orbit along the hexagonal axis; 
bb, cc—closed orbits around the “ring”; dd—open orbit 
parallel to the basal plane. 


The connectivities of these hole surfaces are 
consistent with most features of the observed 
anisotropy and field dependence of the transverse 
magnetoresistance. For example, the experi- 
mental results illustrated in Fig. 2, with the 
current J near the basal plane and the field H 
making an angle @ with the basal plane, show a 
sharp minimum in the resistance at @ =0° in both 
Zn and Cd, which is due to open orbits along the 
hexagonal axis (for example orbit aa in Fig. 1). 

A minimum is observed in Zn also at 6 =90°, 
which is attributed to the lower value of w, for 
orbits around the circumferénce of the ring (or- 
bits bb and cc in Fig. 1) than for orbits embrac- 
ing its arms; the absence of this minimum con- 
firms that the ring is pinched off in Cd. Measure- 
ments in this orientation have not been published 
for Mg, but the minimum observed when H and J 
are coplanar with the hexagonal axis and J makes 
an angle of 25° with the basal plane (and in sever- 
al other specimens of unspecified orientations*) 
may be due to open orbits parallel to the basal 
plane (orbit dd in Fig. 1). For Mg and Cd (and 
probably also for Zn) with J along the hexagonal 
axis and H in the basal plane, the resistance is 
roughly isotropic and varies quadratically with H 
due to open orbits along the hexagonal axis.® 


535 





VotuME 6, NuMBER 10 


PHYSICAL REVIEW LETTERS 


May 15, 1961 






































FIG. 2. Polar diagram of the transverse magneto- 
resistance, Apy/po, in Zn and Cd at 4.2°K. a is the 
angle between the current J and the hexagonal axis, 
and ry the resistance ratio, P73/p,4_ >, in zero field 
(after Justi et al.*). 


A-Cd2, H=12.9 kgauss, 
B-7Zn10, H=12.9 kgauss, 
C-Cd47, H=14.8 kgauss, 


r= 800. 
vr =2100. 
r=2000. 


a = 87°, 
a =88°, 
a =84°, 


A multiply-connected electron surface in the 
third Brillouin zone for Pb, consisting in a peri- 
odically extended zone scheme of thick cylinders 
along the edges of the zone intersecting at its 
corners, was proposed by Gold." With J along 
both [111] and [001] the transverse magnetore- 
sistance of Pb is uniformly high and varies quad- 
ratically with the field, except for sharp minima 
where it saturates when H is along [110]. It is 
difficult to explain this result employing open 
orbits since, for any direction of the orbit, the 
angle a in Eq. (1) must differ from 1/2 for one 
or the other of these current directions. But with 
H along [110] a large closed orbit exists around 
the inside of the multiply-connected surface in a 
single zone, of which there is some evidence from 
magnetoacoustic measurements.’* One might ex- 
pect a large cyclotron mass and hence small we 
for this orbit, and if w,7<1 for the fields and 
specimens used, one has a ready explanation of 
the observed anisotropy of the magnetoresistance. 
But this interpretation is necessarily tentative 
since, although the de Haas—van Alphen measure- 
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ments indicate that the hole surface in the second 
zone is closed,” the behavior of the magneto- 
acoustic attenuation in high fields suggests that 
this surface also is multiply connected.” In any 
case open orbits on the surface in the third zone 
probably occur for some field directions. 

The Fermi surface of Sn is so complex, con- 
sisting of several sheets one or more of which 
may be multiply connected, as to have evaded 
attempts at its determination using several dif- 
ferent experimental methods.'* The anisotropy 
of the magnetoresistance suggests that open or- 
bits occur [for example, the resistance at the 
minimum for one field direction changes accord- 
ing to Eq. (1) as the current direction is varied 
in a plane normal to the field"*], and the predom- 
inantly quadratic field dependence of the resist- 
ance is consistent with the hypothesis that n, =n; 
for all even-valent metals. 

The author has benefited from stimulating dis- 
cussions with his colleagues at the Royal Radar 
Establishment, particularly Dr. J. M. Lock and 
Dr. M. D. Sturge; he is also indebted to Dr. D. H. 
Parkinson for his encouragement ana to D. W. 
Charles for drawing Fig. 1. 
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"This statement is valid even for those field direc- 
tions in a plane of high symmetry for which open 
orbits of type A ® occur, if such an orbit is regarded 
as the limiting form of a highly extended orbit. It is 
not true when the field is in the isolated direction of 
high symmetry which lies at the center of a solid 
angle of field directions for which open orbits of type 
B*® occur (for example, the [100], [110], and [111] 
axes in Au and Cu’). In fact the integral in Eq. (2) 
cannot be defined where orbits of type B occur, and 
one may envisage the transformation of an electron 
sheet into a hole sheet as the volume occupied by 
electrons is expanded to fill more than half the zone 
by supposing the sheet to pass through intermediate 
stages for which orbits of type B occur for all field 
directions. 
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RESOLVED ISOTROPIC HYPERFINE STRUCTURE OF THE ELECTRON 
PARAMAGNETIC RESONANCE ABSORPTION OF F CENTERS IN NaHt 


W. T. Doyle and W. L. Williams 





Physics Department, Dartmouth College, Hanover, New Hampshire 


The wide, inhomogeneously broadened electron 
paramagnetic resonance lines of F centers in 
alkali halides are due to the hyperfine interactions 
between the center electron and the surrounding 
nuclei.’»? These lines are commonly structure- 
less and the hyperfine interactions must be stud- 
ied using the electron nuclear double resonance 
(ENDOR) method.* When the first shell contact 
interaction predominates, however, and the 
nuclear spin of the alkali ion is 3/2, one expects 
a nineteen-line spectrum with relative intensities 
of the components of 1, 6, 21, 56, 120, 216, 336, 
456, 546, and 580, corresponding to total shell 
nuclear magnetic quantum numbers of +9,---, 0, 
respectively.’ In a few crystals, viz., LiF and 
NaF,*»° and RbCl,° a partially resolved spectrum 
has been observed. The structure of the LiF 

line, originally reported to be isotropic,*’® is 

now known to possess many more than the pre- 
dicted nineteen lines’ and to be anisotropic.® 
ENDOR studies® have shown that the first shell 
interaction does not predominate in this case, 

but subsequent calculations’® have established 

that the resolved spectrum is consistent with 

the ENDOR results and with the de Boer model. 

It seems likely that in RbCl, too, the spectrum 

is due to more than the first shell nuclei alone, 


(Received April 21, 1961) 





since the spectrum could not be reconciled with 
the known isotopic abundances of the alkali ions.® 
In NaF, on the other hand, we believe that the 
resolved spectrum is due to the first shell inter- 
action as originally reported.*»> We have found 
that the resolution of the NaF spectrum is unim- 
paired when the sample is in the form of a powder, 
whereas in LiF the resolution is washed out when 
the sample is powdered. NaF is thus the only 
alkali halide in which the hyperfine interaction 
with the first shell nuclei has been seen directly. 
However, even in NaF not all of the nineteen 
lines have been observed, and the signals ob- 
served were neither well enough resolved nor 
strong enough to permit a direct comparison 
with the predicted relative strengths of the nine- 
teen lines. Since the favorable circumstances 
in NaF may be traced in part to the relatively 
large hyperfine interaction constant of sodium, 
and in part to the low atomic number of fluorine, 
we were led to examine the spectrum of the very 
similar crystal NaH, in which it was to be ex- 
pected that the situation should be even more 
favorable for the appearance of a resolved first 
shell structure. 

It is known that the hydrides may be colored 
by radiation and that in LiH the principal product 
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is colloidal lithium when the irradiation is car- 
ried out at room temperature.’ In contrast with 
this, x-ray induced F centers in sodium hydride 
are stable at room temperature, although they 

are thermally bleached within minutes at temper- 
atures around 100°C. The rate of production of 

F centers, as well as the optical and thermal 
bleaching properties, monitored using the electron 
spin resonance signal, were all found to parallel 


closely the corresponding properties in the halides. 


In view of the difficulty of producing single crys- 
tals of NaH, the observations reported here were 
made using commercial sodium hydride powder. 
This has the advantage of tending to smear out 
the small anisotropic contribution. This material 
usually contains an appreciable amount of excess 
sodium, so the samples used were additively 
colored at the start. All of the samples exhibited 
a strong, narrow electron paramagnetic resonance 
(EPR) line which was identified by its shape, g 
value, and temperature independence, as conduc- 
tion electron spin resonance in sodium metal. 

The narrow resonance was not changed by irradi- 
ation and, since it could not be saturated with the 
microwave power available, it was a simple mat- 
ter to subtract it from the easily saturated F- 
center resonance developed by the irradiation. 

In Fig. 1 the F-center resonance alone is shown. 
This curve was obtained by taking the difference 
between a high- and low-power spectrum. The 
shape of the spectrum agrees in detail with the 
predictions of the de Boer model. There are 19 
well-resolved lines. The spacing is determined 
by the first shell hyperfine interaction constant 
A(1). All of the relative intensities agree with 
the theoretical distribution to within the experi- 
mental error of a few percent, as judged by com- 
parison with computed curves of sums of Gaus- 
sian derivatives. It appears from such curves 
that the ratio of peak separation A to the full 
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FIG. 1. Derivative of EPR absorption curve for 
F centers in NaH, showing a well-resolved 19-line 
spectrum of de Boer shape. 
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Table I. Experimental results. The g value has 
been corrected to second order in A/H, using the 
Breit-Rabi equation. The uncertainty quoted for g 
is the rms deviation from the mean of 15 determina- 
tions. The hyperfine interaction energy, A, and the 
modulus of the wave function at the nucleus, ||, for 
the first and second shells were obtained from the 
peak separation and width, respectively. 





|¥(1)|?=4.3 x10" em 
|¥(2)|?=0. 26 x10 em-? 


g&= 1.9979 +0. 0004 
A(1) = 26.5 +0.2 gauss 
A(2) =6.1 +0.3 gauss 





width, 20, between points of maximum and mini- 
mum slope (rms width for a Gaussian shape) is 
about 1.25+ 0.1. Since the ENDOR studies on 
LiF have shown that only the first two shells 
interact strongly with the F-center electron, it 

is probably safe to assume that the width of the 
individual peaks is due primarily to the twelve 
second shell nuclei. When, as in NaH, these have 
spin 1/2, the relationship of the hyperfine con- 
stant to the linewidth is A =o/V3.5 The observed 
values of A as well as the modulus of the F-center 
wave function at the nuclei, |~|?=3/JA/8au, for 
the first and second shell are given in Table I. 
Because they are based upon peak widths rather 
than resolved spacings, the second shell values 
are somewhat less reliable. The g value tabulated 
was corrected to second order in A/H using the 
Breit-Rabi equation. For this purpose the first 
shell is treated as a single entity with shell spin 
I=9.° As expected, the g shift is negative and 
agrees in order of magnitude with the shifts found 
in the halides. 

It appears that sodium hydride is a particularly 
favorable crystal for viewing the hyperfine inter- 
action of the F-center electron with the first shell 
nuclei. The validity of the de Boer model is im- 
mediately apparent in the resolved spectrum. 
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MAGNETIC COOLING WITH PARAMAGNETIC METALS* 
R. D. Parks and W. A. Little 


Department of Physics, Stanford University, Stanford, California 


In this Letter we report on an investigation of 
the magnetic properties at low temperatures of 
a series of thorium-erbium alloys. These alloys 
are unique in that even for high spin densities 
magnetic ordering occurs only at temperatures 
well below 1°K. For this reason it is possible to 
use such alloys as the cooling agents in adiabatic 
demagnetization experiments. 

Conventional paramagnetic salts used for this 
purpose such as potassium chrome alum and fer- 
ric ammonium alum, being dielectric crystals, 
have a thermal conductivity due to the lattice 
which varies as T°. Consequently, at temper- 
atures below 0.1°K this becomes exceedingly 
small, resulting in very low rates of heat trans- 
fer between the paramagnetic salt and the sample 
being cooled. In a metal, however, the thermal 
conductivity which is due to the conduction elec- 
trons varies linearly with T. For temperatures 
below 0.1°K this greatly exceeds that of a dielec- 
tric crystal. In addition a metal is more durable 
than conventional cooling salts which decompose 
on exposure to air or heat. For these reasons 
we have looked for a paramagnetic alloy which 
could be used instead of a paramagnetic salt for 
magnetic cooling. This search led us to the 
study of thorium-—rare earth alloys. As reported 
earlier,’ we found some magnetic cooling below 
1°K with thorium-dysprosium alloys. We have 
now investigated the thorium-erbium system and 
have found excellent cooling properties to tem- 
peratures below 0.1°K. 

The alloys were prepared by melting thorium 
and erbium metals together under argon in an 
arc furnace. This method is described in more 
detail elsewhere.’ A series of Th-Er samples 
were adiabatically demagnetized from initial 
temperatures of 0.73°K and varying magnetic 
fields up to 8400 gauss. The cooling curves 
appear in Fig. 1. When demagnetized from 
fields of about 8000 gauss, all the samples cooled 
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FIG. 1. Magnetic cooling curves of different Th-Er 
alloys. T,is the final temperature reached by the 
samples when demagnetized from an initial temperature 
of 0.73°K and from a magnetic field of H. A residual 
field of 500 gauss (Hy) was left on the sample to sup- 
press either superconductivity or antiferromagnetic 
ordering. The complete removal of H resulted in a 
small increase in temperature which is characteristic 
of the onset of either phenomenon. 


to temperatures below 0.2°K. The 2.3% Er sam- 
ple reached a final temperature of 0.07°K. From 
these curves one should expect that for larger 
values of H/T; (T; = initial temperature) the 
samples should reach even lower final temper- 
atures. In the most dilute sample (0.82% Er) 
the electronic entropy is an appreciable fraction 
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of the spin entropy removed and therefore limits 
T,. In the most concentrated samples 7, is lim- 
ifed by the larger internal interactions which 
lower the entropy. The cooling curves suggest 
that Th-Er alloys with concentrations of Er in 
the range 1-4 atomic percent would be most use- 
ful for cooling below 0.1°K. 

We believe that Th-Er should exhibit better 
magnetic cooling properties than other thorium— 
rare earth alloys. The requirements of a metal 
for magnetic cooling are that it have a large total 
spin entropy per unit volume, N,,,k 1n(2J +1) 
(where N,,, is the number of magnetic ions per 
unit volume), and that it have a low magnetic 
transition temperature, T. (the temperature at 
which magnetic ordering occurs). In the rare 
earth metals the indirect exchange interaction 
between the 4f electrons and the conduction elec- 
trons determines T,.” Liu* and de Gennes* have 
shown that T, should vary as N»»,(J-S)?/J(J+1). 
Of the 14 rare earth elements, these require- 
ments are best satisfied by cerium, ytterbium, 
and erbium. However, Ce5 and Yb® exhibit anom- 
alous behavior which indicates that they would 
not be suitable. 

We have made specific heat measurements on 
the Th-Er alloys which give further evidence that 
the indirect exchange interaction is small. Fig- 
ure 2 shows the specific heat of Th-19% Er in 
the temperature range 0.3-4.2°K. From this 
curve we find that 80% of the total spin entropy 
remains below 0.4°K.”? This is quite remarkable 
since in pure erbium magnetic ordering occurs 
at 84°K.® 

To determine the usefulness of the Th-Er al- 
loys for magnetic cooling, a comparison can be 
made between Th-3 % Er and ferric ammonium 
alum with respect to cooling power, which is 
proportional to the total spin entropy, and ther- 
mal conductivity. Th-3% Er has a total spin 
entropy per unit volume of 33 x10~° joule/cm* 
deg,® compared to that of 53 x 1075 joule/cm' deg 
for ferric ammonium alum. From measure- 
ments of the residual resistivity the value of the 
thermal conductivity, k, due to the electrons, 
was calculated. This gives a value of k =2.8 
x10-* T joule/cm sec for the Th-3% Er alloy. 

At 0.05°K this is about 100 times larger than 
the thermal conductivity reported for compressed 
powdered ferric ammonium alum.’° 

In a cyclic magnetic refrigerator, described 
elsewhere,"* we have used a Th-3.5% Er alloy to 
reach a final reservoir temperature of 0.2°K 
with a starting temperature of 0.73°K. This is as 
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FIG. 2. Specific heat of a thorium-19% erbium alloy. 


low a final temperature as has been obtained by 
using a salt." The refrigerator is now being 
modified for He* pumping which will give a start- 
ing temperature of 0.35°K. By using a Th-Er al- 
loy instead of a salt in this new refrigerator, we 
expect to reach a reservoir temperature below 
0.1°K and obtain high heat extraction rates at 
this temperature. 
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ANOMALOUS LATTICE SPECIFIC HEAT OF SUPERCONDUCTORS™ 


Richard A. Ferrell 
University of Maryland, College Park, Maryland 
(Received April 25, 1961) 


Bryant and Keesom’ have discovered that the 
lattice specific heat of indium is significantly 
less in the superconducting than in the normal 
state. A similar effect has been detected in 
niobium by Boorse, Hirshfeld, and Leupold. 
Daunt and Olsen® have proposed to account for 
the difference in terms of a temperature depend- 


ence of the zero-point energy in the superconduct- 


ing state. But according to the second law of 
thermodynamics, a physical system can only 
take on heat by means of an increase in its en- 
tropy, or disorder. If the high-frequency lattice 
oscillators remain in their ground states, their 
quantum numbers do not change and they cannot 
directly contribute any disorder to the system. 
Their energy is indeed temperature dependent 
because of their interaction with the low-lying 
electron excitations,* but this interaction energy 
is customarily included as part of the electron 
excitation energy and must not be counted twice. 
Thus the zero-point energy does not constitute a 
new and independent source of heat, but instead 
simply contributes to the “renormalization” of 
the single-electron excitation energies. In any 
case, it cannot produce a specific heat anomaly 
at very low temperatures (T<T7,) because of 
the “freezing-out” of the electrons in the super- 
conducting state.® 

The purpose of the present note is to emphasize 
that the specific heat anomaly is clear and un- 
mistakable evidence for an anomalous dispersion 
in the phonon spectrum of the superconducting 
State. Although a more complete report is in 
preparation, we shall also indicate here the a 
priori basis of such a dispersion, and show that 
it gives a good quantitative account of the data. 





In addition we propose several straightforward 
experimental tests of the theory. 

Consider a constant shift upwards of Aw, in 
the range of the angular frequencies involved. 
The specific heat of any given lattice mode is 
kf(iw/2kT), where k and 27h are Boltzmann’s 
and Planck’s constants, T is the temperature, 

w is the angular frequency of the oscillator, and 


F(x) = &/sinhx)?. (1) 


The fractional decrease in the total lattice spe- 
cific heat is 


ac /C, = T. /T, (2) 


where 
T. = 0.556hAw/k, (3) 


and it is assumed that T>T,. The numerical 
factor is one half of the ratio of the sum of the 
inverse cubes of the integers to the sum of the 
inverse fourth powers. To first order in the 
anomaly, the lattice specific heat in the super- 
conducting state is 


C. =aT* - BT, (4) 


where the coefficient of the T* anomalous term 
is B=aT,. 

Figure 1 shows the data of Bryant and Keesom!’ 
divided by T* and replotted vs T. The slope of 
the lines is not adjusted but is instead fixed by 
the value computed by Chandrasekhar and Rayne® 
from their measurements of the elastic constants 
of indium. The dashed line passing through the 
origin represents the Debye specific heat in the 
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FIG. 1. Specific heat divided by the square of tem- 
perature vs absolute temperature, for the supercon- 
ducting state of indium. (The nuclear quadrupole 
specific heat C, is subtracted from the measured total.) 
The solid line is shifted downwards by the dispersion 
resulting from the self-energy process of Fig. 2. 


absence of the dispersion. It is seen that the 
displaced line, which contains only one adjustable 
parameter, its displacement, gives a fit at the 
very low temperatures which could hardly be ex- 
celled by any temperature dependence other than 
that of Eq. (4). The vertical and horizontal inter- 
cepts are B=0.09 millijoule mole deg~* and 
T, =0.065°K, which inserted into Eq. (3) give 
Aw = 1.610" sec™, or a frequency shift of 
2.5x10° cycles per sec. At the higher tempera- 
tures (0.55°K< 7<0.85°K) the difference between 
the heavy line and the experimental points is 
fitted well by the exponential ae~2Tc/T, where 
T, is the transition temperature of 3.37°K and 
the coefficients are a=7.3 and b=1.43. These 
values can be compared with the predictions’ 
of the BCS theory of a=8.5 and b=1.44. The 
good agreement with the BCS exponent is con- 
sistent with the similar agreement with BCS 
theory in the temperature dependence of ultra- 
sonic attenuation found by Morse and Bohm® for 
indium. Turning to the niobium anomaly, we 
find that the data are not fitted quite as well by 
our T* dependence. If we assume, however, that 
the corrections referred to in reference 2 will 
lower the lowest temperature points the most, 
the fit becomes satisfactory. 

The microscopic basis for the shift in the fre- 
quency of the phonons in the superconducting 
state is the well-known result of Bardeen and 
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FIG. 2. Phonon self-energy Feynman diagram. The 
absorption of the phonon produces an electron-hole 
pair which can then re-create the phonon, resulting 
in a frequency shift. The same electron excitations 
are involved in both the real and imaginary parts of 
the frequecny shift. Thus the dispersion in the super- 
conducting state is related to the ultrasonic attenuation 
in the normal state [Eq. (8) of the text]. 


Pines® and others that the frequency of a sound 
wave is very much dependent on the response of 
the conduction electrons. This is illustrated by 
the phonon self-energy Feynman diagram in Fig. 
2. The phonon (represented by a wavy line) con- 
tinually causes electrons to be excited out of 
their normal state. The electrons then act back 
on the phonon, changing its energy. (Time is 
imagined to flow from right to left.) The disper- 
sion of the sound waves depends upon the real 
part of the matrix element corresponding to Fig. 
2. We want to calculate the change in the real 
part when the metal passes from the normal to 
the superconducting state. It is useful, however, 
first to consider the imaginary part, which rep- 
resents the damping calculated by Kittel,’® and 
which can be written as w/47Q, where w is the 
angular frequency of the sound wave, and Q is 

a dimensionless number measuring the phonon 
mean free path in units of the wavelength. Kittel’s 
result is 


Q=hpc /(m*cV, (5) 


where 27h is Planck’s constant, p is the mass 
density of the metal, c, is the velocity of sound, 
m* is the electron effective mass, and C is the 
usual electron-phonon coupling constant." It is 
clear that there is not only a contribution to the 
imaginary part of the phonon frequency, but that 
a real contribution must also result, from virtual 
excitation of the same type of electron states as 
contribute to the damping. 

The general features of Kittel’s formula have 
been verified by Morse and Bohm,® including the 
frequency independence of @. They have also 
shown that at temperatures very much below the 
transition temperature the Kittel damping disap- 
pears completely. This result follows readily 
from an energy gap picture of superconductivity: 
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If the phonon energy is less than the gap, there 

is no continuum of excited states into which a 
transition can take place. But this same de- 
ficiency of states in the gap will also lead to a 
loss in the virtual excitations, and thereby pro- 
duce a change in the real part of the phonon fre- 
quency. Introducing #@,, as the effective energy 
gap within which the normal density of electron 
excitations per unit excitation energy «, p(e), is 
replaced by zero in the superconducting state, 

we can write the difference in the complex phonon 
self-energy between the normal and superconduct- 
ing state as the following integral: 


Pe 2% 
haw is)-2 | 


H,, is the matrix element for the production of 

an electron-hole pair by the absorption of a sound 
wave of wave number g and is independent of ec, 
while 6 is an infinitesimally small positive quan- 
tity. Since p(e) vanishes for « =0, we can use the 
Taylor series expansion p(e)=«p’, where p’ is 
the derivative. [We now restrict our attention to 
values of q sufficiently large that the maximum 
value of « for which p(e)+0 is much larger than 
hG,; i.e., we consider wavelengths considerably 
smaller than the coherence length.] Both |H q y 
and p’ can now be taken outside the integral sign, 
and elementary integration yields 


hog lH, |* plede 


Gia -te-10r 2 ©) 





Aw + (iw/4 7Q) = Ho |?p"(2 @ + inw). (7) 


We verify that insertion of the usual expressions 
for H gq and p’ reduces the imaginary part of this 
equation to Eq. (5). The real increase in the 
sound-wave frequency in the superconducting 
state can now conveniently be compared with the 
imaginary part, giving us our basic equation, 


Aw = @ /2n°Q). (8) 


This is a type of dispersion relation, relating the 
dispersive behavior of the metal to its absorptive 
properties. 

It is now desirable to see if Eq. (8) yields a 
reasonable value for @. For our present pur- 
poses, we can approximate the effective gap by 
the actual energy gap, h@ =hw o* 4kT,. Sum 
rule arguments indicate that this is probably an 
overestimate, but not by more than 50%.” With 
this sort of accuracy understood, we can write 


Q=0.117 /T., (9) 


from which we obtain for indium, Q~=5.8. Equa- 
tion (5) enables us to go further, and deduce the 
mean electron-phonon coupling strength. From 
the Debye temperature® of 111.3°K we infer c, 

= 0.83 x 10° cm/sec, while the normal state 
specific heat of indium’ gives m* equal to 1.21 
times the free electron mass. We find C=7.9 ev, 
which happens to be just equal to the Sommer- 
feld-theory Fermi energy. This result is not 
unreasonable for a metal with a complicated 
energy band structure sensitive to lattice 
changes, and indicates that a given percentage 
elastic strain must change the individual single- 
electron energies in about the same proportion. 
Since all superconductors do not necessarily 
have such large deformation potentials, it is 
clear that the specific heat anomaly cannot al- 
ways be expected to appear. Furthermore, the 
lowest shear mode, with which the anomaly is 
primarily associated, has a longitudinal admix- 
ture for nonsymmetry directions. The resulting 
longitudinal deformation potential and electric 
field could in some cases interfere destructively 
with the shear coupling constant for the low- 
energy electron excitations. Therefore, it is 
clearly not feasible to predict a priori which 
superconductors will exhibit the anomaly. Never- 
theless, it is possible to measure the quantity 

Q directly from ultrasonic absorption (for the 
lowest velocity mode averaged over propagation 
directions) and thereby employ Eq. (9) to pre- 
dict the value of T;. This would yield a quanti- 
tative prediction of the heat anomaly in advance 
of actual thermal measurements. Conversely, 

it would be highly desirable to have a check on 
the present considerations by means of an actual 
ultrasonic measurement of Q for indium (and, of 
course, also niobium for which T, ~0.2°K and 
Q~=5). Such a measurement has already been 
made on the (001) shear wave in tin,’* yielding 
Q~=60. Assuming that this is representative of 
the other propagation directions, we can predict 
that the specific heat anomaly is an order of 
magnitude smaller in tin than in indium. A fur- 
ther indication of an especially strong coupling 
of the shear waves to the electrons in indium is 
the strong temperature dependence found by 
Simmons and Slichter™ in the nuclear quadrupole 
resonant frequency. 

An alternative ultrasonic experimental check 
on the present considerations would be a direct 
measurement of the dispersion in the supercon- 
ducting state, rather than of the attenuation in 
the normal state. Wavelengths \ smaller than the 
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coherence length ~, are required, because of the 
general principle that the bulk elastic constants 
must give a correct description of sound waves 
for \> t,. (The very small normal-supercon- 
ducting difference in the elastic constants results 
in turn from the fact that the transition energy 
depends parametrically on the lattice strains 
and forms only an extremely small fraction of 
the total ground-state energy.) Since the co- 
herence length can be estimated for indium at 
é, = 6300 A, the minimum frequency for detect- 
ing appreciable dispersion is about 1.3 x 10° 
cycles per second. As a consequence of this 
minimum frequency, the lattice specific heat 
must return to its undisplaced Debye value in 
the limit T~0. That is, if it were possible to 
extend the measurements to such low tempera- 
tures, the points in Fig. 1 would leave the solid 
line and rejoin the dashed line. 

It should also be noted that the constant fre- 
quency shift follows from Eq. (6) only for w«<@ 
At higher frequencies a more complicated dis- 
persion results, leading for w > vg to a very 
small shift [reduced by the order of (Gg/w)*}. 
Consequently, inelastic neutron scattering or the 
Mossbauer effect, for example, does not offer a 
feasible alternative method of detecting the dis- 
persion. 

An additional experimental test of Eq. (8) is 
suggested by Pippard’s*® study of the effect of 
the electron mean free path / on the ultrasonic 
attenuation. He finds that the Kittel damping is 
reduced by a factor proportional to //, for l< ), 
in agreement with similar studies by Kittel,’° 
Mason,’* and Morse.*” This suggests that if a 
metal showing the heat anomaly were heavily 
enough doped with an alloying impurity to reduce 
significantly the coherence length (as would be 
indicated by a decrease in the superconducting 
skin depth), then the anomaly would disappear. 





*Research supported in part by the Air Force Office 
of Scientific Research and by the Office of Naval Re- 
search. 
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ANOMALOUS CAPACITANCE OF THIN DIELECTRIC STRUCTURES 
C. A. Mead 


California Institute of Technology, Pasadena, California 
(Received March 9, 1961; revised manuscript received April 10, 1961) 


The classical analysis of a plane parallel capac- 
itor gives a capacitance which becomes infinite 
as the spacing between the metal electrodes ap- 
proaches zero. The fact that real capacitors 
show a marked deviation from this behavior has 
not generally been recognized.'»? This communi- 
cation describes experiments in which the devia- 
tions have been measured with some care. 

In order to determine the dependence of capaci- 
tance upon thickness it is necessary to have a 
reliable measure of the relative thickness of the 
device and also a method for producing very thin, 
uniform dielectric layers. Normal methods for 
measuring thickness are quite inconvenient for 
the thin layers used here. Consequently the 
voltage required to obtain a given tunneling cur- 
rent density*»* was used as a measure of the 
thickness. This procedure is considered to be 
a good approximation if the difference in metal- 
insulator work functions of the two metals under 
the conditions used is small, since the tunneling 
current density is only a function of the electric 
field at the interface, and not of the thickness of 
the dielectric. Investigation of the slight asym- 
metry of the volt-ampere characteristic of the 
structures used indicated that differences in the 
metal-insulator work functions were indeed 
small. 

The present experiments were all carried out 
on a tantalum sheet which had been recrystalized 


“a 


0.15- 


FIG. 1. Inverse capacitance 
as a function of tunnel voltage 
for Ta-Ta,O,-Au structures. 


in a dry argon atmosphere at 2700°C for 1/2 hour. 
The sheet was then composed of many crystals 
approximately 1 mm across, and was subsequently 
electropolished to a mirror finish. Tantalum 
oxide was formed on the surface anodically in 
dilute ammonium citrate solution. The samples 
were then placed in a vacuum system where they 
were baked at 300°C at residual pressure of 
3x10-° mm Hg while counter electrodes of gold 
were evaporated through a mask. The resulting 
structure was a large number of gold squares 
approximately 5 mils (0.005 in.) across on the 
tantalum oxide surface. Measurements of ca- 
pacitance and tunnel voltage were made by stand- 
ard techniques on a large number of squares on 
each sample. Only squares lying wholly on one 
crystal were used. The tunnel current density 
used was approximately 0.05 ma per square mil. 
Previous experience had indicated that a very 
small amount of contamination of the surface re- 
sulted in large variations of capacitance between 
squares on one sample. On the samples tested it 
is believed that surface contamination has essen- 
tially been eliminated since variations were in 
many cases less than 1 % and in all cases were 
well within the limits placed upon the measure- 
ment of the variation of area between squares. 
The results of this experiment are shown in Fig. 
1. It is seen that instead of having zero intercept, 
the plot of inverse capacitance versus tunneling 
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FIG. 2. Inverse capacitance 
as a function of tunnel voltage 
for Ta-Ta,O;-Bi structures. 
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voltage intercepts at approximately 0.05 square 
mil per puf. This intercept corresponds to an 
€, layer approximately 2.75 A thick. It was felt 
that perhaps variations in the capacitance with 
applied de electric field could be observed. The 
experiment was performed but no detectable var- 
iation was observed up to electric fields where 
tunneling current made the measurements unre- 
liable. 

Theoretical work® indicates that the dielectric 
constant of an insulator should be preserved down 
to a few atomic layers. If this is true, we should 
expect the basic limitation of the capacitance to 
be electric field penetration into the metal. A 
crude classical calculation using a method simi- 
lar to that of Mott and Jones® gives a penetration 
depth of the order of one A for typical metals, 
varying inversely with the square root of the 
density of states in the metal. Although no real 
confidence can be placed in such a calculation, 
if this is indeed the mechanism of the capacitance 
limitation, the 1/C intercept should vary with the 
density of states in the metal. 

This dependence was investigated by another 
set of experiments using bismuth counter elec- 
trodes. Because of its low density of states, 
bismuth should have a larger value of penetration 
depth and hence a higher intercept on the plot of 
1/C versus tunneling voltage. That this is indeed 
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the case is shown in Fig. 2. Here the intercept 
corresponds to an €, layer approximately 5.5 A 
thick. 

At present it has not been possible to achieve 
truly symmetric structures and hence the series 
capacitance contribution for each interface is not 
known, but only the sum of those for two inter- 
faces. A symmetric structure plus a knowledge 
of the true dielectric constant within the metal 
would be necessary for a full determination of 
the parameters involved. 

The author is indebted to H. M. Simpson for 
preparing the samples, and to C. H. Wilts and 
R. W. Gould for helpful suggestions and comments. 
The work was supported in part by grants from the 
International Telephone and Telegraph Company, 
the General Electric Company, and a contract 
with the Office of Naval Research. 
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AMPLIFICATION OF 9.3-kMc/sec ULTRASONIC PULSES BY MASER ACTION IN RUBY 
E. B. Tucker 


General Electric Research Laboratory, Schenectady, New York 
(Received April 27, 1961) 


The observation of the attenuation of 9.3-kMc/ 
sec ultrasonic pulses by the spin-phonon inter- 
action! in ruby implies that amplification should 
be observable in the same system if the spin 
population is inverted. We have now observed 
the amplification of 9.3-kMc/sec pulses of ultra- 
sonic energy in passing through a ruby rod in 
which the spin system was inverted. This is the 
first observation of amplification of energy other 
than electromagnetic in such a system. 

The conditions for the observation of maximum 
gain from the inverted spin system are identical 
with the conditions for the observation of the max- 
imum attenuation except for the inversion, rather 
than the normal thermal equilibrium, of the spin 
system. The published results for attenuation in 
ruby’ show that one of the highest ultrasonic 
interactions occurs, for longitudinal waves, at 
an angle of 60° between the magnetic field and 
the ruby c axis (@ = 60°) for the 2-3 transition 
at high field (3700 gauss for 9.3 kMc/sec). For- 
tunately, the conditions for use of push-pull 
pumping” (saturation of both the 1-3 and the 2-4 
transitions) at a single frequency occur at 0 
= §5°30’. The actual orientation used for the 
experiment is a compromise between the above 
requirements. Since the ultrasonic transition 
probability is a slow function of angle near 
6=60°, it was found that the best inversion was 
obtained at an angle of 56° using pump power of 
23.3 kMc/sec. 

The experimental observations were made 
using the same methods previously described’ 
except that the ruby was contained in a K-band 
pump cavity. The crystal configuration consists 
of a 15-cm X-cut quartz crystal bonded, with 
indium, to a 1}-cm length of Linde pink ruby 
(Cr*** in ALO,). Both crystals are in the form 
of cylindrical rods 3 mm in diameter and the c 
axis of the ruby is parallel to the rod axis to 
within } degree. The experiments were carried 
out at 1.5°K. 

The traces of Fig. 1 illustrate the effects of 
the spin system on the first few echoes observed. 
The echo to the right of the center graticule has 
made a single return trip through the ruby while 
the echo to the left of center has made two passes 
through the ruby. The larger amplitudes of these 
two echoes in the top picture as compared to the 
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FIG. 1. Echo patterns observed for 9.3-kMc/sec 
ultrasonic pulse of 0.5-microsecond duration. Large 
echoes to the right and the left of the center graticule 
have passed through the ruby once and twice, respec- 
tively. Time increases from right to left. Top picture: 
magnetic field on resonance line and pump on. Center 
picture: magnetic field on resonance line, no pump. 
Bottom picture: magnetic field off resonance line. 


bottom trace is due to the amplification through 
interaction with the spin system. Similarly, the 
smaller amplitudes of the center trace as com- 
pared to the bottom are the result of attenuation 
due to transfer of energy to the spin system in 
thermal equilibrium at 1.5°K. The gain observed 
is 0.12 per centimeter of travel in the ruby.® 
This is somewhat higher than the attenuation of 
0.09 per centimeter and is an indication of the 
efficiency of inversion. 

The slower traces of Fig. 2 illustrate much 
more graphically the magnitude of the effect. 
The upper picture shows clearly a number of 
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FIG. 2. Echo patterns observed for 9.31-kMc/sec 
ultrasonic pulse of 0.5-microsecond duration. The 
change in frequency accounts for the amplitudes of the 
first few echoes being different from those of Fig. 1. 
Sweep speed of oscilloscope is about 5 microseconds 
per centimeter. A round trip through the ruby requires 
2.1 microseconds. Time increases from right to left. 
Top picture: magnetic field on resonance line and pump 
on. Bottom picture: magnetic field off resonance line. 


echoes not visible without the aid of the amplifi- 
cation. From the spacing of these echoes it can 
be surmised that they are due to pulses which 
have traveled back and forth in the ruby for as 
many as eighteen round trips; and that their only 
traversal of the quartz has been the 3 cm travel 
to get to the ruby initially and, after the reflec- 
tions in the ruby, to get back to the transducer 
end again. Echoes from the bond are almost non- 
existent after 10 microseconds. The beating due 
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to phase interference caused by the nonparallel 
ends and bond as well as crystal inhomogeneities 
and misalignment is very evident in the top trace 
of Fig. 2. 

The physical perfection required to determine 
the true propagation characteristics of the quartz- 
indium-ruby combination is lacking, and hence 
conclusion of whether or not the gain obtainable 
from the spin system is sufficient to overcome 
the system losses is impossible. It is fairly ob- 
vious what conditions must be satisfied in order 
to realize over-all gain, i.e., a phonon maser. 
Just as in the case of the optical maser‘ the gain 
from the inverted population must be sufficient 
to overcome the losses due to attenuation in the 
materials and reflections at the ends of the cavi- 
ty, in this case the optically polished ends of the 
crystals. If one assumes that an attenuation of 
0.01 per centimeter (such as is observed in 
quartz*) is attainable in a good ruby crystal, the 
gain of 0.12 per centimeter is certainly sufficient. 
For appreciably higher losses either higher spin 
populations or spin systems with larger magneto- 
elastic coupling constants, or perhaps both, will 
be required. 

This experimental demonstration of phonon 
amplification was undertaken after discussions, 
on the possibility of a phonon maser, with Dr. 

E. H. Jacobsen and Dr. N. S. Shiren. It isa 
pleasure to acknowledge their contributions. 
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'The pictures illustrate amplitudes proportional to 
voltage whereas the attenuations and gains are referred 
to power. 
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QUADRUPOLE MOMENT OF Li’ AND QUADRUPOLE COUPLING CONSTANT OF Li,t 


S. L. Kahalas and R. K. Nesbet* 
Department of Physics, Boston University, Boston, Massachusetts 
(Received April 17, 1961) 


Recently, the quadrupole coupling constant of 
Li’ in LiH has been measured to high accuracy 
by Klemperer, Gold, and Wharton.’ When com- 
bined with our calculated value of the electric 
field gradient at the Li’ nucleus, this measure- 
ment leads to a definitive value for the quadrupole 
moment of Li’, Q(Li’) = -3.56x1077* cm?. We 
have also calculated the electric field gradient in 
Li,. Using Q(Li’) calculated from LiH, we find 
the quadrupole coupling constant for the Li, mole- 
cule to be +0.049 Mc/sec compared with the value 
of +0.060 Mc/sec estimated by Logan, Coté, 
and Kusch.? This would indicate that the Logan, 
Coté, and Kusch value may be too large and that 
it would be desirable for the quadrupole coupling 
constant of Li, to be remeasured with higher ac- 
curacy in order to check the prediction. The 
remainder of this Letter presents some details 
of the molecular calculations. 

Our electronic computer calculations of the 
molecular properties of the LiH and Li, mole- 
cules have been carried out by the Roothaan* 
(Hartree-Fock, self-consistent field, molecular 
orbital, linear combination of atomic orbitals) 
procedure as modified by Nesbet.*»> Extended 
sets of basis functions of the Slater type were 
employed and configuration interaction correc- 
tions to the energy and electric field gradient 
were calculated. For LiH, a basis set of 13 
Slater-type orbitals were used, including pz 
functions on both the lithium and hydrogen cen- 
ters. For Li,, 18 Slater-type orbitals were used 
(9 on each center), including pz functions. In 
both cases the pm functions were found to give 
large contributions in configuration interaction. 
For configuration interaction calculations in LiH, 
the 27 most important configurations are included 
in the energy and 19 configurations in the electric 
field gradient. For configuration interaction cal- 
culations in Li,, the 28 most important configura- 
tions are included in the energy and 24 configura- 
tions in the electric field gradient. These calcu- 
lations were carried out at three internuclear 
distances, including the internuclear distance 
corresponding to the minimum energy. 

The accuracy to which the molecular wave 
function has been calculated may be indicated 
from the following results. For LiH, a total 


molecular energy of 99.3 % of the experimental 
value, and a dissociation energy of 89.3 % of the 
experimental value were calculated. Similar cal- 
culations in Li, resulted in 99.4% and 89.9%, 
respectively, for these two quantities. The en- 
ergy of atomic Li, which enters the dissociation 
energy calculation, was calculated to comparable 
accuracy (99.44). 

Recently Klemperer, Gold, and Wharton have 
measured the quadrupole coupling constant of Li’ 
in the LiH molecule’; the value eqQ/h =+346+ 2 
kc/sec is found. Our calculated value for the 
electric field gradient at the Li nucleus is g/2e 
= -0.02070 ay;~°, evaluated at the internuclear 
distance corresponding to the minimum of the 
energy. é is the magnitude of the electron charge 
and ay is the first Bohr radius. This leads to a 
value of -3.56x10~** cm? for Q(Li’), in agree- 
ment with predictions from theoretical models 
of the nucleus.® This value is estimated to have 
an error associated with it of less than 10%. As 
distinct from the Li, case,’ the electronic and 
nuclear contributions to the electric field gradient 
are not nearly equal, the latter being about 0.65 
of the former. 

For Li,, q/2e is calculated to be -0.0029 ay~. 
Together with Q(Li’) found above, this leads to 
a quadrupole coupling constant of +0.049 Mc/sec. 
The experimental value of Logan, Coté, and 
Kusch,’” which is +0.060 Mc/sec, would thus ap- 
pear to be too large. 

A complete account of the molecular calculations 
will be presented later. 

We wish to thank Professor W. Klemperer for 
permission to quote the LiH quadrupole coupling 
constant result prior to its publication. One of 
us (S.L.K.) gratefully acknowledges several help- 
ful discussions with Dr. R. E. Watson. We wish to 
thank Professor N. F. Ramsey for a discussion 
concerning some measurements made in his lab- 
oratory. Professor P. Kusch is to be thanked for 
his correspondence concerning the Li, quadrupole 
coupling constant. 
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chusetts Institute of Technology, and RIAS and the 
Martin Company at Baltimore, Maryland. 
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INELASTIC SCATTERING OF 18.9-Mev NUCLEONS FROM THE 9.6-Mev STATE OF C”” 


E. Bradford and B. A. Robson 
Research School of Physical Sciences, The Australian National University, Canberra, Australia 
(Received April 17, 1961) 


Recently it has been suggested’ that the 9.6-Mev 
state of C’* has spin and parity J7=3~ and not 17 
as previously supposed.” This Letter reports a 
new analysis of Peelle’s* 18.9-Mev inelastic pro- 
ton scattering data, which definitely favors 3~ 
for the 9.6-Mev state. The calculation uses a 
direct volume interaction with spin-dependent 
distorted waves.‘ 

Peelle’s® analysis of his data is based upon the 
direct surface interaction theory, which predicts 
a jx” (RR) angular distribution. Using an interac- 
tion radius R =3.3 fermis, he found a best fit for 
K=1 giving J"=0", 1°, or 2°. The same theory 
applied to inelastic scattering from the 4.4-Mev 
2* state of C’* gave very poor agreement with 
experiment. 

Levinson and Banerjee® showed that the direct- 
surface-interaction theory is probably inadequate 
for a nucleus as small as C, and Robson and Rob- 
son* found spin-orbit effects to be important for 
12-Mev nucleons inelastically scattered from the 
4.4-Mev level of C’*. Thus it was considered es- 
sential to take both these effects into account be- 
fore using Peelle’s data to determine the spin and 
parity of the 9.6-Mev state. In the following, only 
J"=1~ or 3” are considered because there exists 
a reasonable amount of evidence’ for eliminating 
other values. 

The distorting potential is taken to be the usual 
Woods -Saxon optical potential plus a Thomas- 
Fermi spin-orbit potential,® 

2 
Vy”) =-(V+iW)f(r) + V> + A a) . os. L, 


f(r) ={1 + exp[(r -R) /a]}". (1) 
The values of V, W, Vg, and Wg used are 45, 12, 
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FIG. 1. Inelastic nucleon scattering from the 9.6- 
Mev state of C'*. The points are the experimental re- 
sults of Peelle for 18.9-Mev protons. The theoretical 
curves are for 18.9-Mev neutrons. The normalization 
is arbitrary and both curves are fitted at 75°. 
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15, and -4 Mev for the incident energy of 18.9 
Mev and 45, 8, 17, and -4 Mev for the emergent 
energy of 9.3 Mev, with a=0.4 fermi and R 

= 2.75 fermis in both cases. 

The direct-interaction two-body potential is 
assumed to be zero ranged and spin independent, 
and all forms of exchange are neglected. 

Following Barker et al.,’ the C’ ground state 
is taken as 


%o= W(1s* 1p*[4, 4]000, 0), (2) 


where the numbers following the configuration are 
values of [A]TSL,J. The 9.6-Mev state is taken 
as 


Ys = W((1s* 1p7[4, 3341, 1d)003, 3) (3) 
for J"=3", and 


3 
Y= 2 ads; (4) 
i=1 


for J"=1~, where 
¥,, =V((1s* 1p"[4, 3}541, 2 s)001, 1), (4a) 
ro = ¥((1s* 1p"[4, 3}331, 1.d)001, 1), (4b) 
drs = V((1 s°[3 40, 1 p°[441]$41)001, 1). (4c) 


For the Elliott and Flowers interaction, Barker 


et al.’ found a, = 0.916, a, =-0.399, and a, = -0.028. 


For simplicity, the values used in the present 
calculation are a, = 0.869, a,=-0.494, and a, 

=0. The 1p, 1d, and 2s wave functions are taken 
for an ideal harmonic oscillator with length para- 
meter 6 =2.0 fermis, obtained by fitting the 1p 


harmonic-oscillator wave function to the 1p 
square-well solution of the Schrédinger equation 
for a radius R =3.66 fermis, which gives a rea- 
sonable fit to the 4.4-Mev data.‘ 

Figure 1 shows the result of the calculation. It 
is seen that the theoretical curve for the 3° as- 
sumption is a much better fit than for the 1° as- 
sumption (both curves are normalized at 75°). 
Greater emphasis should perhaps be placed on fit- 
ting the scattering for angles @<90°, since other 
processes like heavy-particle stripping’ may be- 
come important for @>90°. On account of the very 
lengthy calculations involved, no attempt was 
made to vary the parameters, which were se- 
lected by a consideration of the elastic scatter- 
ing and inelastic scattering from the 4.4-Mev 
state data.*»® 
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In this note we present preliminary data ob- 
tained by exposing the magnet cloud chamber 
to the beam of K,° particles from the internal 
target of the Joint Institute for Nuclear Research 
synchrophasotron. Neutral particles emitted 
from the target at an angle of 97° to the direc- 
tion of the circulating proton beam are used in 
the experiment. The cloud chamber is placed 
at a distance of 8 m from the target. 

We have recorded so far more than 500 V° 
events interpreted as decays of K,° mesons. 
Based on the selective analysis of experimental 
data, an estimation was made of the relative 
probability of a Kg3 decay and quantitative 
data on the charge ratio in three-body decays 
were obtained. Only 225 V° events are measured 
completely. There is not a single 27 decay a- 
mong them. We have found one four-prong event 
and four electron-positron pairs (with compara- 
tively large opening angles) emitted at a large 
angle to the direction of K,° flight. 

(1) The passage of the charged decay products 
through a lead plate 5.8 g/cm’ thick located in 
the operating volume of the chamber was studied 
in order to estimate the relative probability of 
K,3 decay. Out of 112 events we have found 24 
events in which a particle either stops in the 
plate or undergoes a momentum loss of more 
than 30%. Making use of the Bethe-Heitler equa- 
tion’ (improved by Eyges”) which provides the 
probability of the given energy loss by an elec- 
tron for radiation in traversing the definite lead 
thickness, we have found that the fraction of K,3 
decay with respect to all the decays with charged 
products constitutes 38%.*° Taking into consid- 
eration the corrections for K,°-particle motion 
which leads to a relative increase in the number 
of pions and » mesons traversing the plate, we 
determined the fraction of K,3° to be equal to 42 
+12%. 

If one considers that the relative probability of 
3 n° decay of K° mesons is not more than 20% for 
K,° mesons,* then the probability of K,3 decay 
for these particles turns out to be equal to the 
twofold probability of K,3 decay for positive K 
mesons’ within experimental error. Such a ratio 
between the two probabilities should occur if 
the isotopic spin selection rule® AT =} is valid. 
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DECAY PROPERTIES OF K,° MESONS* 


D. Neagu, E. O. Okonov, N. I. Petrov, A. M. Rosanova, and V. A. Rusakov 


Joint Institute of Nuclear Research, Moscow, U.S.S.R. 
(Received April 20, 1961) 

























































Therefore, the experimental data obtained point 
out, within statistical deviations, the extension 
of the above-stated selection rule to the decay 
processes involving leptons. 

(2) Quantitative data on the charge ratio (i.e., 
the ratio of the number of decay events with the 
emission of 7- mesons to the number of events 
accompanied by the emission of 7+ mesons) are 
obtained both by identifying the charged decay 
products having momentum not larger than 120 
Mev/c and by studying the passage of decay par- 
ticles through the lead plate. In measuring the 
ionization, the electrons can be reliably selected 
when the momentum of decay particles does not 
exceed 100 Mev/c. When the momentum is be- 
tween 100 and 120 Mev/c, it is only possible to 
determine whether a particle under study is a 7 
meson or not. 

The number of decays obtained which were ac- 
companied by the emission of 7- mesons and m+ 
mesons is equal to 46 and 51, respectively. They 
provide the value 0.90+ 0.18 for the charge ratio. 
Within experimental errors this ratio does not 
differ from unity and is in disagreement with the 
analogous ratio 0.42+0.27, estimated from the 
data of Bardon et al.’ 

Neither the authors of the quoted paper nor we 
have found any event of K,°-meson decay into two 
charged 7 mesons. 

Combining our data with those obtained in refer- 
ence 7, we set an upper limit of 0.3 % for the rel- 
ative probability of the decay K,°—1-+m+. Our 
results on the charge ratio and the degree of the 
2-decay forbiddenness are in agreement with 
each other and provide no indications that time- 
reversal invariance fails in K° decay. 

(3) Up to now only three V° events’~® which 
might be interpreted as the decay of K,° mesons 
according to the mode K,°~1-+7++m° have been 
recorded. We have found one four-prong event’® 
which can be considered as the decay of a long- 
lived K° meson into two charged and one neutral 
m meson. In view of the presence of the Dalitz 
pair, in our experiment the nature of the neutral 
product of the decay is determined more direct- 
ly than in the case of two-prong events; the ob- 
served four-prong decay is a direct proof for 
the existence of K,°-meson decay according to 
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Table I. Characteristics of observed electron- 
positron pairs. 











Momentum (Mev/c) Opening angle 
(+) (-) 
1 55 42 7°’— 9° 
2 10 : 43 99° 
3 111 103 18°—20° 
4 26 79 25° 








the above-mentioned mode. 

(4) Among all the recorded V° events there were 
found 4 electron-positron pairs with comparative- 
ly large opening angles; the directions of three of 
these pairs make angles greater than 40° with that 
of the K,°-meson beam. The characteristics of 
these pairs are summarized in Table I. 

There is not a single star with an electron-posi- 
tron pair from the stars found in cloud-chamber 
gas. On the other hand, the probability of large 
opening angles for electron-positron pairs aris- 
ing in the external conversion of y rays is very 
small. These are the reasons to consider the 
pairs obtained to be Dalitz pairs resulting from 
K,°-meson decay into neutral 7° mesons. There- 
fore, the fact of their observation can be treated 
as experimental evidence for the existence of 
the long-lived K°-meson decay into three 7° mes- 
ons. The observed quantitative ratio (1:4) of four- 
prong decay and single Dalitz pairs can also be 
considered as a support of this treatment. 

The authors are grateful to Professor B. M. 
Pontecorvo for his constant interest in the work 
and for the discussion of the results. They would 





also like to thank Professor V. I. Veksler and 
Professor V. P. Dzhelepov for their interest in 
the experiment and their support. 
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exactly as it was submitted to the Rochester Conference, 
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distorted. See D. Neagu et al., Proceedings of the 1960 
Annual International Conference on High-Energy Physics 
at Rochester (Interscience Publishers, Inc., New York, 
1960), p. 603. 

1H. A. Bethe and W. Heitler, Proc. Roy. Soc. 
(London) A146, 83 (1934). 

*L. Eyges, Phys. Rev. 76, 264 (1949). 

3Here the correction for 7-meson nuclear absorption 
is already taken into account. 

‘L. M. Lederman, 1958 Annual International Con- 
ference on High-Energy Physics at CERN (CERN 
Scientific Information Service, Geneva, 1958), Vol. 5, 
p. 275. 
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Stevenson, Phys. Rev. Letters 2, 361 (1959)] with the 
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ELASTIC SCATTERING OF MUONS IN NUCLEAR EMULSION* 


P. L. Connolly,t J. G. McEwen, and J. Orear 


Laboratory of Nuclear Studies, Cornell University, Ithaca, New York 
(Received April 24, 1961) 


In the past 10 years several experiments have 
been performed to measure the nuclear scatter- 
ing of cosmic-ray muons.’ Many of these exper- 
iments are claimed to show, in addition to the 
Coulomb interaction, an anomalous muon-nuclear 
interaction. These cosmic-ray experiments suf- 
fer from several difficulties: (1) They measure 
inelastic as well as elastic scatterings; (2) the 
pion background is uncertain and difficult to cor- 
rect for; (3) the incoming muon energy is not 
known; and (4) the use of thick scattering plates 
requires serious corrections for multiple scat- 
tering (also for successive single scatterings 
combined with multiple scattering). At high mo- 
mentum transfer these scattering corrections 
require knowledge of the elastic and especially 
the inelastic cross sections. Such theoretical 
knowledge is yet only approximate.” Masek, 
Heggie, Kim, and Williams® have recently per- 
formed an experiment which overcomes most of 
the above difficulties. They measured a mixture 
of elastic and inelastic scattering of 2-Bev y~ 
on carbon using a Bevatron separated muon beam; 
most of their scatterings with high momentum 
transfer should be inelastic, and their results 
are in agreement with the inelastic cross section 
estimates of Drell and Schwartz.‘ 

Our experiment, on the other hand, measures 
only elastic scatterings as do the electron scat- 
tering experiments of Hofstadter.° The corre- 
sponding calculations can be done precisely so 
that detailed information on the charge distribu- 
tion of the nucleus can be obtained from the elas- 
tic scattering measurements. A departure of the 
muon elastic scattering from these precise cal- 
culations would reveal a nonelectromagnetic 
muon-nuclear interaction, or a difference be- 
tween the mean radii of the muon and electron 
charge distributions. 

In this Letter we report the final results of 
our u* scattering experiment and preliminary 
results of our u~ scattering experiments. The 
u* stack was exposed to 3.5x10° muons in the 
43-Mev separated muon beam of the Carnegie 
Tech synchrocyclotron and the y~ stack was 
exposed to 2.510" muons in the 60-Mev sepa- 
rated muon beam of the CERN synchrocyclotron. 
In the u.* experiment the small pion background 
is eliminated by observation of the 7-: decay. 
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In the .~ experiment the ratio of muon endings 
to pion endings is greater than 1000 to 1. Nega- 
tive pions are also eliminated by observation of 
prongs at either the scattering or the ending. We 
have observed no inelastic scatterings. Most in- 
elastic scatterings in our energy region would be 
expected to excite the giant resonance and thus 
show up as a large increase in grain density at 
the scattering.® Other inelastic scatterings could 
be distinguished by proton or conversion-electron 
prongs at the scattering. 

Large effective path lengths are achieved by 
area scanning for muon endings in the region 
where muons are not supposed to end and tracing 
these endings back to a possible large-angle scat- 
tering. Since the geometry of the stack is known 
precisely, it is possible to make numerical cal- 
culations of the effective muon path scanned. 

In our + experiment we have restricted our 
data to 14- to 40-Mev muons scattering at angles 
between 80° and 180°. Our energy and angle lim- 
its correspond to momentum transfers of 80 to 
160 Mev/c. Figure 1 shows the number of events 
observed having a momentum transfer greater 
than the value g. The curve is the number of 
events predicted from the elastic scattering cal- 
culations of Rawitscher,’ who used the nuclear 
charge distributions obtained by Hofstadter.°® 
Figure 2 shows our experimental cross sections 
plotted against momentum transfer gq, together 
with Rawitscher’s calculated cross sections for 
27-Mev u* scattering elastically on Ag (Z =47) 
and Br (Z =35) averaged together. Rawitscher’s 
cross sections when plotted against g are virtu- 
ally independent of muon energy in the energy 
region studied. The statistical standard devia- 
tions are shown in Figs. 1 and 2, and it appears 
that the data are consistent with Rawitscher’s 
calculations. The predicted number of events 
depends both on our geometry calculation and 
interpolations of Rawitscher’s cross sections. 
We estimate an over-all accuracy of +10% for 
our calculation of the predicted number of events. 
This possible systematic error is not shown in 
Figs. 1 or 2. 

Early scanning on our ~ experiment has yield- 
ed 10 events so far with g>100 Mev/c and this 
also appears to be consistent with Rawitscher’s 
calculations. 
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FIG. 1. The number of events having a momentum 
transfer greater than g is plotted vs g. The curve is 
the predicted number of events based on the elastic 
scattering cross sections calculated by Rawitscher . 


We conclude that our positive muons behave as 
heavy positrons in the same momentum transfer 
region as was studied by the cosmic-ray experi- 
ments. We have no indication of an anomalous 
muon-nuclear interaction. If we make allowance 
for a possible finite size of the muon by multiply- 
ing Rawitscher’s cross sections by f*(g), where 
f(g) =(1+ 8Ry av?) is the form factor due to 
the muon charge distribution,® we find that our 
data determine an upper limit of 1.7107 cm 
for the rms muon radius ((R py av)”. 

We wish to thank Professor R. Siegel and Pro- 
fessor R. Sutton of Carnegie Tech and Dr. A. Ci- 
tron and Dr. D. Michaelis of CERN for their gener- 
ous help with the exposures. Professor G. Raw- 
itscher of Yale University has kindly performed 
and forwarded to us additional calculations more 
appropriate to our experiment than those which 
appeared in his original paper.’ 





*Research supported by a National Science Founda- 
tion research grant. 
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FIG. 2. One-half the sum of the differential cross 
sections for the elastic scattering of 27-Mev u* on Ag 
and Br plotted against the momentum transfer g. The 
experimental points are based on a total of 138 observed 
scatters in nuclear emulsion. The curves for extended 
and point charge distribution are those calculated by 
Rawitscher. 


TtNow at Brookhaven National Laboratory, Upton, New 
York. 

‘For the most recent cosmic-ray experiments see 
R. L. Sen Gupta, S. Ghosh, A. Acharya, M. M. Bis- 
was, and K. K. Roy, Nuovo cimento 19, 245 (1961). 
For a summary of the experiments up to 1958, see 
G. N. Fowler and A. W. Wolfendale, Progress in 
Elementary Particle and Cosmic-Ray Physics (North- 
Holland Publishing Company, Amsterdam, 1958), 
Vol. 4, p. 123. 

°. N. Cooper and J. Rainwater, Phys. Rev. 97, 
492 (1935). 

3G. E. Masek, L. D. Heggie, Y. B. Kim, andR. W. 





555 











VoLuME 6, NUMBER 10 


PHYSICAL REVIEW LETTERS 





May 15, 1961 











Williams, Phys. Rev. 122, 937 (1961). 

4s. D. Drell and C. L. Schwartz, Phys. Rev. 112, 
568 (1958). 

5R. Hofstadter, Ann. Rev. Nuclear Sci. 7, 231 (1957). 





SHans Bethe (private communication). 

'G. H. Rawitscher, Phys. Rev. 112, 1274 (1958); 
also Bull. Am. Phys. Soc. 6, 301 (1961). 

85. D. Drell, Ann. Phys. 4, 75 (1958). 










PROPOSED METHOD OF MEASURING THE SPIN OF THE K’ MESON 


M. Schwartz* 


Columbia University, New York, New York and Brookhaven National Laboratory, Upton, New York 
(Received April 27, 1961) 


The recent discovery’ of the K’ meson with a 
mass of about 880 Mev has ied to considerable 
speculation about its role in strange particle 
reactions. A knowledge of its spin is hence of 
obvious importance. It has already been verified 
that the K’ has isotopic spin 1/2 and decays via 
the modes K’+K+z. At the moment, one possi- 
bility for determining whether its spin is 0 or 1 
is by observing the occurrence (or nonoccurrence) 
of the decay mode K’+K+y. It is the purpose 
of this note to point out another method of meas- 
uring the spin of the K’, which, although no bet- 
ter than the above, is rather unique and perhaps 
somewhat elegant. 

The reaction we shall consider is p +p~+K°+K°’ 
or K°’+K°, where the antiproton capture takes 
place at rest. Only one assumption will be made 
here—namely that the capture proceeds from an 
S§ state. The justification for this assumption 
rests on the recent arguments of Day, Snow, and 
Sucher.’ At any rate, there is an inherent check 
of this assumption in the experiment as will be 
seen below. 

Subject to the above assumption, the initial 
state can be either the triplet or the singlet 
state. The triplet state is odd under charge 





K 9 °° om” 0 0.0 
Ki +K Kia 





Y= tr +02 (9) | “—— = 


We shall deal now only with the first part of this 
state. Let us substitute 


K° =(K,°+K,°)/V2, K°=(K,°-K,°)/V2. 
We then have 


- Or Op o 
¥: (Kia - x. Ky, )/ v2. 


1b 








conjugation and the singlet state is even. The 
parity of the initial system is odd. 

The reaction }+p+K°+K°. It is well known 
that if the above conditions are satisfied, then 
the outgoing K°K® system must be odd under 
charge conjugation. If this is so, the outgoing 
wave function will have the factor Kj,°K2,° 
-K2q°K1p°. Hence one K° will always decay 
via the K,° mode and the other will decay via 
the K,° mode. The observation of this effect 
serves then as a check on the assumption of 
S-state capture. 

The reaction +p —+-K°+K°’ or K°’+K°, in the 
event that the K°’ has spin zero. In this case, 
the outgoing state must be an S state (since the 
relative parity of the K and K’ is odd). Hence 
the initial state must be the singlet S state. 
Therefore it is even under charge conjugation. 
Hence the final state must be of the form K°K°’ 
+K°'K®. The K°’ (or K°’) decays via the mode 





K°! » (3)"*Kq° - (3)"*K*n, 
(K”’ ‘ait (3)"K°r® = (3)*K~1*). 
If we now adopt the notation that subscript a re- 


fers to the primary K° or K° and b refers to the 
other, we can rewrite our outgoing state as 


K 


+ = 0s 
% “9 +K Kn 
a a 


v2 





o 








Hence we will observe either a two-K,° type of 
decay or a two-K,° type of decay, but no events 
in which we have both a K,° and a K,°. | 

The reaction }+p—-K°+K°’ or K°’+K°, in the 
event that the spin of the K’ is 1. In this case 
the outgoing state must be a p state. Hence we 
can have either J=0 or J=1. J=2 is forbidden 
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since there is no J=2 initial state. The J=0 part is even under charge conjugation and the J=1 part 


is odd. Hence we have 





na (12 15 ee wa) (Eu 20 


where A, refers to the J=0 part and A, refers 
to the J=1 part. 

Turning to the problem of determining the 
spin of the K’, let us first select that group of 
5 stoppings which show a K,° decay with a mo- 
mentum such that the other outgoing particle is 
a K’ (or K’) meson. This can be done above 
background to the extent that the K’ mass is 
well defined. Then we can deduce the following 
points: 

1. If the spin of the K’ is zero, then another 
K,° decay must also occur. 

2. If the spin of the K’ is 1, then the other K° 
can decay via the K,° mode some of the time, 
unless A, happens to be zero. 

3. In addition, we can unambiguously deter- 
mine the spin by considering those cases where 
both decays are observed to be of the K,° type. 
Let us consider the z axis along the direction 
in which the K’ was moving before decay. Events 
of the above sort must have J=0 regardless of 
the spin of the K’. The orbital angular momen- 
tum of the KK’ (or K’K) system must have z 
component equal to zero. Hence the K’ must 
have its spin in the m =0 state. If we now con- 
sider its decay into 7°+K° in its own center-of- 





0 0 
Ma's» ) 


mass system, the outgoing particles will have 
the following distributions about the z direction: 

(a) If spin K’=0, the distribution is obviously 
isotropic. 

(b) If spin K’=1, the distribution must have 
the form cos*6@, where @ is the polar angle be- 
tween 7° and K’. 

In the event that there is a nonresonant K,°K,°n° 
background, as there most certainly will be, 
then the distribution which one will obtain as a 
result of the above procedure will have the form 
|\>> a;P;°(cosé@)|*. In principle, one can then 
determine the coefficients a; from the distribu- 
tion and observe whether a, or a, varies signifi- 
cantly in passing through the resonance. If a, 
varies, then the spin of the K’ is zero. If a, 
varies, the spin is 1. 

The author would like to express his gratitude 
for valuable discussions with Dr. M. Nauenberg, 
Dr. G. Feinberg, and Dr. J. Steinberger. 








jAlfred P. Sloan Research Fellow. 
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‘?. B. Day, G. A. Snow, and J. Sucher, Phys. Rev. 
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PION-LAMBDA RESONANCE (Y,*)* 


J. P. Berge, P. Bastien, O. Dahl, M. Ferro-Luzzi, f J. Kirz, D. H. Miller, 
J. J. Murray, A. H. Rosenfeld, R. D. Tripp, and M. B. Watson 


Lawrence Radiation Laboratory and Department of Physics, University of California, Berkeley, California 
(Received April 27, 1961) 


Alston et al.’ have discovered an J=1 1-A reso- 
nance with a mass of 1380 Mev and a half-width 
I'/2 consistent with 32 Mev. They reported on 
141 An+n- events made by 1150-Mev/c K~ mesons 
incident on the 15-in. Lawrence Radiation Labo- 
ratory hydrogen bubble chamber which produced 
the sequence 


As *t F 
K +p-Y, +7, (1) 


followed by a strong decay, 


aE t 
Y, ~A+m +130 Mev. (2) 


In the course of a continuing study of Kp inter- 
actions, we have now explored the same reactions 
from Y,* threshold [P;(lab) =405 Mev/c] through 

=850 Mev/c. We report on ~500 Amta- and 
find m(Y,*) =1385 Mev and I’/2 closer to 20 Mev; 
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however, we are still unable to distinguish between 


spins J=1/2 and J=3/2. 


It is conventional to discuss Y* data in terms 
of a simplified model in which they are produced 


Ty+, c.m. (Mev) 


and decay “isolated,” that is, one neglects the 


interferences: (a) between Y,*’s and background 
Ant*tn~ produced in nonresonant partial waves, as 


well as (b) the influence of Bose statistics and 
final-state interactions on the m7 system. 


Isolated Y,*’s will be produced in association 
with a pion of fixed center-of-mass (c.m.) kinetic 


energy (spread, of course, by +I'/2). Thus in 
Fig. 1, which is Dalitz’ representation of the 
Antn- data at Py-(lab) = 850 Mev/c, Y,** will fall 


in a horizontal band and Y,*~ in a vertical band. 
That this model is poor is shown by the fact that 
conservation of parity requires that our isolated 
Y,* show no fore-aft asymmetry on strong decay, 


i.e., in the strong-decay distribution in the Y* 


c.m. frame, 


dn /dQ = 1 +a, (A+ ¥,*) +a,(A- ¥,*)? +eee : 


the odd coefficients must be zero. However, our 
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FIG. 1. Dalitz plot of K™ 
+p—A+ 1*+77 for 185 events 
at Px(lab) = 850 Mev/c. Note 
that within the Y,*~ band, for 
example (i.e., for constant 
T,+), Tq- is linear in (A-¥*>, 
where A is the unit vector 
along the A momentum in the 
Y* rest frame (7~=-A). Thus 
A’s which are produced for- 
ward by the strong decay of 
the Y*~ (Y*---1~ +A) fall at 
the lower edge of the ellipse, 
and vice versa. Point A cor- 
responds to maximum, and 
point B to zero, A kinetic en- 
ergy. 





data show a concentration of backwards A. The 


asymmetry coefficients a, (assuming a, =0) in 
Table I show that the “isolated” model is far 


from realistic, both for our new data and the 


older events at 1150 Mev/c.' A major purpose 


of this Letter is to point out this difficulty. 
In the following Letter,” Dalitz and Miller treat 


the reaction as a symmetrized three-body state 
which permits nonzero a,. However, some of 
the features of the isolated model are still dis- 


cernible, for example the extra population of the 


Y,* bands on the Dalitz plots. Therefore, in order 


lar distribution. 


Table I. Coefficient a; (%) in Y,* strong-decay angu- 
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+ + 
P , (lab) (Mev/c) ¥, Y, 
760 -24 +20 -16 +20 
850 -92 +26 -24 +24 
1150 -70 +26 - 2220 
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to simplify the discussion, in our Letter we still 
refer to this “isolated” model. 

The data to which the model may be most mean- 
ingfully applied are those at Px =850 Mev/c, 
where the Y,* bands do not overlap badly (see 
Fig. 1). (They cross in the middle of the Dalitz 
plot at about 700 Mev/c.) Figure 2 is an experi- 
mental histogram of the Y,* mass. The individual 
uncertainty in each mass measurements is typi- 
cally +(3 to 5) Mev. If we fit the histogram with 
an s-wave resonance curve of the form dn/dm 
«[(m - 1385)" +(I'/2)?]~ as indicated in Fig. 2, 
we find ['/2=15 to 20 Mev (depending upon whether 
or not some background is subtracted). This 
experimental value is an agreement with the “KN 
bound state” prediction of '/2=18 Mev,° but it is 
also not too far from global symmetry’s T/2 
= 25 Mev (see below). 

Figure 3 displays the excitation data for ye 
and Y,*~ and total An+n-. Data at 1150 Mev/c are 
from Alston et al.’ and at 300 and 400 Mev/c are 
combined from the present experiment and from 
Nordin.‘ The Y,* cross sections at 850 Mev/c 
can be obtained fairly unambiguously in the sense 
that the An+m- events of Fig. 1 can be interpreted 
as a flat background plus an extra population of 
Y,*’s at the bands. A background subtraction is 
also possible at 760 Mev/c. However, at 620 
Mev/c, the bands cover almost the whole ellipse, 
making it impossible to distinguish Y,* events 








from background. At 510 Mev/c, the bands no 
longer cover the middle of the ellipse, and we 
might expect any appreciable Y,* population to 
show up in the bands above background. In the 
Dalitz plot of the data for Px =510 Mev/c, we 
actually find, however, a concentration of events 
near the center of the diagram, not in the bands. 
This result could be consistent with essentially 

no Y,* production (as suggested in Fig. 3) plus a 
nonresonant background of An+m~- events tending 

to favor equal pion momenta (for reasons unknown); 
it is also conceivable and quantitatively sensible 
that there is indeed appreciable Y,* production 
with the Y,* and the production pion predominantly 
in a relative p state (or higher). This would tend 
to depopulate the ends of the ellipse where the 
production pion momenta would be low and produce 
the observed distribution of events. Therefore, 
one cannot determine unambiguously from these 
data the cross section for Y,* production at Px 
=510 Mev/c. 

Table II illustrates the apportioning of events 
between the various possible reactions. 

The Y,* production angular distributions are 
isotropic within statistics for all our beam mo- 
menta. This would suggest that Y, production 
proceeds dominantly through a single J=1/2 wave. 
However, the Amt+tzm~- cross section alone comes 
very close to 7147/2, and so other partial waves 
are probably present. 
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FIG. 2. Experimental his- 
togram of the Y,* mass, m(Y,"*), 40 
at Py (lab) = 850 Mev/c. To 
each event plotted in Fig. 1 
there correspond two effec- 
tive Y;* masses, m(Y,**) and 
m(Y,*"); we choose as m(Y,*) 
that which is closer to 1385 20 
Mev. S-wave resonance 
curves of the form dn/dm 
=[(m - 1385)?+ (I'/2)"]~* have 10 
been fitted to the data as in- 
dicated, with T/2=10, 15, 
and 20 Mev. 
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Next we discuss the Y,* spin. Global symmetry 
predicts a Am resonance with spin J=3/2.5 On 
the other hand, the Dalitz-Tuan s-wave resonance 
in the K~p system corresponds to a Y,* with J 
=1/2, anda strong decay via S,,. if the KA parity 
is odd.® If our model of an isolated Y* were cor- 
rect, the strong decay angular distribution of 
Eq. (3) should be fore-aft symmetric (a, =0), and 
the presence of any polar-equatorial anisotropy 
(a, #0) would be evidence for J>1/2. We find 
instead that a, #0 but within statistics a,=0 at 
all production angles. After looking for evidence 
for a,#0 at any (or all) angles of Y* production, 
we have confined our analysis to polar-produced 
Y*’s where, as Adair has pointed out,’ an isolated 


| 
% 100 500 500 400 500 600 700 800 900 000 100 120 


unique interpretation of 
the data. 


J=3/2 particle must decay according to 1+ 3(A-K)’ 
(the unit vector K is along the K~ beam direction) 
while J=1/2 must decay isotropically. The decay 
distribution for polar-produced Y*’s 3 (1¥*: K| 20.80) 
is given in Fig. 4. We have used i- ¥* as a meas- 
ure of the decay angle instead of A-K in order to 
display the fact that the Y* are not isolated but 
tend to undergo strong decay with the A going 
backwards. (Essentially the same result is 
achieved for these polar events by choosing +A-K 
for forward-produced Y*, and -A-K for backwards 
y*.) These Adair events of Fig. 4 again fail to 
show any evidence for 1 +3(A- Y*)? superimposed 
on the asymmetry, which is linear in (A- Y*). 
This statement applies individually to both charges 


Table II. Numbers and cross sections for V and two-prong events. Numbers in parentheses are cross sections 








in millibarns. 
An*1~ 
Py (Mev/c) Y, aa ~" Three-body L'nta An*n~1® R°p1- Total 
510 ———— 31 (1.2 +0.3) ——— 4 (0.2+0.1) 0 0 35 
620 anaes 54 (1.8 +0.3) a 10 (0.3 +0.1) 0 (0) 0 (0) 64 
760 121 (1,240.14) 122 (1.2+0,14) 56 (0.6+0.1) 55 (0.6+0.1) 20 (0.2+0.05) 2(0.03+0.02) 376 


850 56 (0.9 +0.25) 


62 (1.0+0.25) 67 (1.0+0.15) 39 (0.6 +0.1) 


7 (0.1+0.04) 3(0.03+0.02) 234 














{ 











es 














VoLUME 6, NUMBER 10 


PHYSICAL REVIEW LETTERS 


May 15, 1961 




















} _— = 4 








Number 














O J | ! | 
-10 -O6 -02 O2 O6 +10 
A-yY* 





FIG, 4. Adair analysis of 62 Y** events at P ,= 760 
and 850 Mev/c; |¥**K| 0.80. 


of Y* at each momentum and to the sum of all 
events. However, as is shown in the next Letter, 
one cannot conclude that J=1/2 since examples 

are given where Bose statistics can cause a Y* 
with J=3/2 to have a relatively isotropic angular 
distribution in the Adair analysis. 

Next we discuss the branching ratio Y,* += 

vs A. Global symmetry predicts definite ratios 
between the decay rate I'(N) of the 3/2, 3/2 pion- 
nucleon isobar (N* +N +7) and the equivalent 

rates of, for example, the three decay modes 

of Y,*~. These are I'(N): re": r(=~):T(A) =py’: 
(1/8) 50°: (1/6)p,.-*:(4/6)p,°. For py =230 Mev/c, 
bs =127, p, =207, (py /py* $=1/4, the relative 
rates are stediates to rhe 7. 37:(1/24):(1/24):(2/3). 
The half-width ['/2(N) is known to be 45 Mev; so 
global symmetry predicts a Y,* half-width I'/2 

= 3[I'(D°) +I(=~) +I (A)] of 24.7 Mev, and a branch- 
ing ratio T(D°)/T(A)=1/16. (As already mentioned, 
our observed I'/2 is 15 to 20 Mev, but the statis- 
tics are sufficiently poor and the effects of Bose 
statistics sufficiently great that '/2=25 Mev 

cannot be ruled out.) The Dalitz-Tuan resonance 

is expected to give values of I'(=°)/I'(A) ranging 
from 212% (if the KA and KZ parities are both 
odd) down to 6% (if KA odd but KE even).* We 
obtain P'(D°)/T'(A) <3 % at 760 Mev/c and <5% at 
850 Mev/c; these values are obtained by making 
a background subtraction and represent “realistic” 
upper limits. “Maximum possible” upper limits 
come from counting every £°n+n- event with an 
effective 2°r mass falling within 30 Mev of the 
Y,* mass of 1385; these are <20% at 760 Mev/c 





and <10% at 850 Mev/c. The observed limits on 
the branching ratios are consistent with the pre- 
dictions of global symmetry and not necessarily 
inconsistent with those of the Dalitz-Tuan reso- 
nance. 

The average Y* polarization (Py*) has been 
measured from the asymmetry of the A-decay 
proton assuming @ = +1 and Y*xK= up; only 
events with | ¥*- K|<0.85 have been used, giving 
Py+-= = (+114 21) % and Pys+= (-16+ 21) % at 760 
Mev/c; Py«-=(-56+ 20) % and Py++= (+124 28) % 
at 850 Mev/c. The fact that the Y*~ polarization 
seems to be large at 850 Mev/c and, by compari- 
son, surprisingly small at 760 Mev/c could be a 
statistical fluctuation or may be connected with 
the fact that the >-production channels also show 
a rapid variation with energy in this region.® 

Since no strong evidence for Y* spin J=3/2 
has yet been found either at Berkeley or else- 
where,®?° it is intriguing to assume J=1/2 and 
to try to distinguish between S,,. and P,,. Y* strong 
decay. 

We use about 60 polarized Y*~ at 850 Mev/c 
(assuming that the polarization is not a statistical 
fluctuation). In addition to the normal fi to the 
production plane, it is convenient to define another 
unit vector m=2(n-A)A-n. This vector m lies in 
the plane of nandA (the direction of flight of 
the A), and if A makes an angle 6 with n, m makes 
an angle 26.° We write P,, for the measured A 
polarization along fi and, similarly, Pi. along m. 

If the Y* polarization is P, then for strong decay 
via S,., one can show 


P =P, P =-P/3, P /P =-1/3. (5a) 
n m m n 


For P,, decay the roles of fi and m are inter- 
changed; i.e., the A has its maximum polariza- 
tion along m, and we have 


* =-P/3, P =P, P/P., =-3. (5b) 


Py find P,,=(- 56+ 20)%, P m= (+33 + 25) %, and 
P,,/P y= 0. 6. It is not meaningful to state statis- 
tical errors on a calculated ratio like P,, /P,, when 
the fractional standard deviations are large. In- 
stead we apply the x’ test; the average value of x’ 
should be 1.0. For the S,, hypothesis, we find 
Xs’ =0.3 (high probability), but for Py, xp*=4.5 
(probability~3%). Martin et al.° also report P,, 
= -0.38+ 0.25 and P,, =0.19+ 0.25, which yield 
xy =0, xp?=1. Thus if it were established (a) that 
is the Dalitz-Tuan resonance of a K+p in an 
S,,. bound state and (b) that Eqs. (5) are not badly 
perturbed by interference phenomena, then our 


561 















VoLuME 6, NuMBER 10 


PHYSICAL REVIEW LETTERS 





May 15, 1961 





data would strongly suggest odd KA parity. 
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Experimental evidence for the occurrence of 
an J=1, A-7 resonant state (Y*) as an intermed- 
iate step in the reaction 


t- + 
Y +7 


K"+p- +Atntn (1) 
y** we" 

has been presented by Alston et al.’ for K” (lab) 
momentum 1150 Mev/c. The same reaction has 
also been studied recently by Berge et al,,- and 
the related K,°-p reactions have been analyzed 
by Martin et al.° for K,° momentum 975 Mev/c. 
The interpretation of these data in terms of the 
mass (M*), half-width ([/2), spin (J), and pari- 
ty of the Y* state has been confused by the evi- 
dence that the Y* decay in these reaction se- 
quences (1) cannot be regarded as the decay of 
a free particle. In this Letter, we show that 
this evidence can be largely understood as due 
to interference effects arising from the require- 
ment of Bose statistics for the final pions, and 
we discuss the extent to which these Y* param- 
eters may ultimately be determined from data 
in this momentum range. 

From the reported estimates (10 to 30 Mev) 
for ['/2, the mean distance traveled by the pri- 
mary pion in one Y* mean lifetime is more than 
4 fermis. It is therefore a plausible assumption 
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that the successive pion-emission processes do 
not interfere dynamically to any marked degree. 
In this case we are led to an amplitude M(1, 2) 
for the reaction sequence K+N~+Y*+7,, Y*+A 
+1, of the general form* 


M(1, 2)=(,q, p,, P,)A(P,), (2) 


where q and p, denote the c.m. momenta of the 
K~ meson and the primary pion, and P, denotes 
the momentum of the secondary pion in the 7,-A 
rest frame. The final configuration A +7,+7, 
may also be specified by giving /, the orbital 
angular momentum of 7,, and L, the orbital 
angular momentum in the 7,-A system which 
corresponds to the Y* spin and parity; the form 
of @ is then determined for total angular mo- 
mentum j by the angular momentum coupling 
and the (KA) parity and by centrifugal barrier 
considerations. The role of the resonant state 
is represented by the second factor, 


A(P) = exp[i6(P) sind(P)/P?” ait (3) 


The amplitudes for the two sequences (1) must 
be added coherently, and their sum must cor- 
respond to a final state with correct symmetry 
for interchange of the two pions. Thus, for total 
isotopic spin J=0 and 1, the amplitudes M, and 














M, are given by 
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M,=M"(1,2)+M"(2, 1), (4a) 
M,=M"(1, 2) -M""(2, 1). (4b) 


For comparison with experiment, it is con- 
venient to consider the probability distribution 
P(E,,E-) on the (E,,E~) phase-space diagram. 
Since the distribution P(E,,E_) sums over all 
orientations of the An*tn™ plane, the contributions 
to P(E,,E~) from initial states of different angu- 
lar momentum and parity add incoherently. Since 
P(E,,E_) =3|M,+M,\? and P(E_, E,)=31M,-M,|’, 
the symmetrized probability distribution 
{P(E,,E.)+P(E-,E,)} is simply a superposition 
of the separate J=0 and /J=1 distributions, with- 
out interference between them. The observation 
that the experimental distributions P(E,,E-) and 
P(E.,E,) do not differ markedly from each other 
at 760 and 850 Mev/c suggests that, for these 
two momenta, the reaction may be dominated by 
a single isotopic-spin channel. 

Two An‘n™~ configurations are of particular 
interest: 

(A) E,=E_, and maximum A recoil energy. 

For this configuration, the two pions have the 
same c.m. momentum, and M(1, 2) and M(2, 1) 
are necessarily equal. For 850 Mev/c and below, 
the Y** and Y*” bands overlap strongly, and 
marked interference between these two terms 

is to be expected at A (see Fig. 1 of preceding 
Letter”)—constructive for ]=0, destructive for 
I=1. The experimental distributions at 760 and 
850 Mev/c show an especially low density of 
events in the region near A, which suggests that 
it is the J=1 channel which is dominant at these 
momenta.® The density of events expected in 
this region for the models discussed below is 
compared with the data in Table I. 

(B) E, =E-, and the A hyperon at rest. For 
this configuration, the two pions have equal and 


opposite momenta in the c.m. system. With 
I=1 (@=0), M(1,2) and M(2,1) interfere con- 
structively or destructively according as 1+ L 
is odd (even) or even (odd). The observed dis- 
tributions at 760 and 850 Mev/c indicate strong 
constructive interference in the region near B, 
leading to the conclusion that, with J=1, 1+L is 
odd.® An S,,. A-a resonance therefore requires 
a p-wave primary pion at these production en- 
ergies, a P,, or P,, A-n resonance requires an 
s- (or d-) wave primary pion. 

We now discuss in detail the calculations of 
the symmetrized (E,,E-) distributions for some 
cases of particular interest, to illustrate the 
nature of the interference effects resulting from 
Bose statistics. These cases represent only the 
simplest possibilities for fitting the available 
data, and are neither unique nor exhaustive ex- 
amples. 

S..2 A-m resonance. With the assumption of 
p-wave excitation, we have 


M(1, 2) = (aq: p, + ibo- 4xp,)(M,-M"+il/2)™, (5) 


where M, denotes the total energy in the A-7z, 
rest frame. Figure 1(a) compares the calcu- 
lated distribution of M(A-7) for M* =1385 Mev 
and several values of ['/2 with the 850-Mev/c 
data (obtained by projecting the symmetrized 
(E,,E£-) plot onto one axis). These curves dis- 
play an insensitivity to [/2, for an adequate fit’ 
to the data is obtained with any value of ['/2 
between 20 and 30 Mev. As shown in Table I, 
the destructive interference obtained for the 
I=1 pS,, case is in reasonable accord with the 
data; J=0 sS,, production, on the other hand, 
gives too high a density of events near A, as 
well as a poor fit to the mass distribution.® 
The mass distribution calculated for J=1 pS,, 
production also gives an adequate fit to the ex- 
perimental data! at 1150 Mev/c, for the same 


Table I. The number of A+2*+727 events expected in the region near A, with both E, and E_ less than the limit 
Em Stated, for the K~-p interaction data presently available. * 





Configuration (JL Py 





Events 
(sSy2)y2" (PP y2)v2° (pPy2)y2° (pSy2)y2' (sPy2)y2' observed 
850 Mev/c; Ey =150 Mev 35.8 26.0 8.5 12.3 11.0 11/262 
Em =125 Mev 29.8 32.8 8.6 9.7 10.5 18/252 
760 Mev/c; 
Em =120 Mev 17.4 21.1 4.2 3.1 2.2 8/252 





~ 
From references 2. 
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FIG. 1. The M(7-A) distribution observed at 850 
Mev/c is compared with (a) the calculated distributions 
for the J=1 (pSy») and I=1 (sPy,) configurations with 
M*=1385 Mev and various values of '/2, and (b) the 
calculated distributions for the J=0 configurations sSyy», 
(pPy2v2 and (pPy)y_ with T/2=25 Mev. 


M* and '/2=25 Mev. J=0 pS,,. production can 
also fit these data but, the interference at B 
being destructive, only with a larger width ['/2. 
Since this K~ energy lies just above the strong 
1=0 K-N resonance [at 1100-Mev/c K (lab) 
momentum], there is, of course, little reason 
to believe that this Y production should neces- 
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sarily involve the same states as are effective 
at 850 Mev/c. In fact, these data show appreci- 
able asymmetry between y** and y*~ production, 
an indication that both J=0 and J=1 production 
must be occurring from the same angular-mo- 
mentum states. 

Py, A-n resonance. For J=1 production, we 
consider s-wave excitation, 


M(1, 2) = (2q- P, -io -qxB)A(?,), (6) 


where the phase shift 5(P,) has been taken to 
have the form known for the (3,3) resonance,’ 
with suitable parameters M* andI. The fit to 
the 850-Mev/c mass distribution is quite ade- 
quate with [/2=25 Mev, as shown in Fig. 1(a); 
the same parameters again give an adequate fit 
to 1150-Mev/c distribution. Figure 1(b) shows 
the calculated distributions for J=0 production 
in the j = 1/2 and 3/2 states (pP,,). Neither of 
these configurations reproduces the data. Table 
I shows that for (pP,,.),2, the constructive inter- 
ference near A is too strong; (pP,2),. does give 
small intensity in this region because the y*- 
decay distribution (relative to production direc- 
tion) happens to have approximately (7-6 cos?é6) 
form for this case, but this angular distribution 
is incompatible with the data (see Fig. 2). Al- 
though a rough fit may be obtained by superpos- 
ing these two configurations, it is clear that this 
comparison adds support to our assignment of 
I=1 production. 

P,, A-n resonance. For s-wave excitation, 
we have 





M(1, 2)=0- P,A(P,), (7) 


where A(P,) is taken to have the form appropriate 
to a P-wave resonance, that used already with 
Expression (6). The symmetries required for 
fitting the (E,,£~) plots are similar to those for 
the P,, possibility just discussed. { 

The highly anisotropic distribution observed 
for Y* decay relative to its direction of motion 
may also be understood in terms of Bose statis- 
tics. For J=1 and1+L=odd, the configuration 
with the secondary pion emitted forwards gives 
constructive interference, so that the angular 
distribution for A emission in the Y* decay is 
peaked backwards. The decay angular distribu- 
tions (averaged over all production directions), 
calculated for (sP,,.) and (pS,,.) j for 850 Mev/c, 
are shown in Fig. 2 and compared with the ex- 
perimental data. 

The decay angular distribution for Y* pro- 
duction at 0 and 180 deg is of particular inter- 
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FIG. 2. The angular distribution of the pion from 
Y* decay, relative to the Y* direction of motion, for 
all events at 850 Mev/c with M(1-A) between 1360 and 
1410 Mev. The calculated curves for J=1 sPy, and 
I=1 pSy, production reproduce the general trend rather 
well. The curve for J=0 (pPy2)y32 production is plotted 
to show that the data strongly exclude the possibility 
that this is the dominant production state. 


est, in view of Adair’s general result for its 
dependence onJ. Again, the angular distribu- 
tions characteristic of isolated decay are strong- 
ly modified by the requirements of Bose statis- 
tics. The distortion of the idealized (1+3 cos*@) 
distribution for J =3/2 is shown in Fig. 3(a). 
The difficulty in distinguishing conclusively be- 
tween the cases J =1/2 andJ=3/2 at 850 Mev/c 
is apparent.'° The character of these distortions 
is closely linked with the requirement of con- 
structive interference at B. 

As the K~” momentum is increased further, 
the primary pion energy will exceed the energy 
of the secondary pion and Bose interference ef- 
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FIG. 3. Adair distributions are shown for forward- 
and backward-produced ig decay (integrated over the 
mass range 1360 =M(A-7) = 1410 Mev) at K” momenta 
850 Mev/c and 1150 Mev/c, where the angle @, of the 
A particle in the Y* rest frame is measured from the 
direction of motion of the Y*. The distortion from the 
expected (1+3cos*@,) distribution for J= 3 is much 
weaker at 1150 Mev/c than at 850 Mev/c. However, 
the present statistics at 1150 Mev/c do not yet allow 
the cases J=4 and J=}4 to be distinguished. 


fects will become less and less important. In 
the calculations for 1150 Mev/c, this is already 
apparent.’ As shown in Fig. 3(b), the distortion 
of the Adair distribution for J=3/2 is much less 
severe; unfortunately the paucity of data at this 
momentum does not yet allow a decision on the 
Y* spin. It is clear that, in further experiments 
at higher K~ momenta, the effects of Bose sta- 
tistics will be sufficiently reduced to allow the 
possibility of a clear-cut distinction between 
J=1/2 andJ=3/2 for the Y* state. 

For Y* production from a K~-p state of angular 
momentum j and isotopic spin J, it is important 
to note that the geometric limit on the absorption 
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cross section is (2j +1)7X?/4, which takes the 
value 2.85(j+4) mb at 850 Mev/c. At this mo- 
mentum, the cross section observed for Re- 
action (1) is 3.2+ 0.3 mb, comparable with the 
geometric limit for j=1/2. This limitation pro- 
vides a severe condition” which must be satis- 
fied by any simple model of the production proc- 
ess. For J=1, this condition is satisfied by 
(sP,,) production, but not by (sP,,) production. 
The possibility (pS,,),. [together, perhaps, with 
some (pS,,.),,. production] also satisfies this con- 
dition. At this point we note that the angular dis- 
tributions for Y* production are isotropic at 760 
and 850 Mev/c. This is consistent with (sP,,) 
production, of course. For the (pS,,),, case, * 
a calculation of the angular distribution (includ- 
ing Bose effects) shows that a strong cos*6@ term 
is to be expected at these energies; this dis- 
crepancy could be reduced by a large admixture 
of (PS, 2). excitation of suitable phase. The 
large value observed for the (A+7*+77) pro- 
duction cross section is most probably an in- 
dication that the distribution receives contri- 
butions from a number of partial waves and 
configurations, of which only the dominant terms 
have the symmetry specified in our discussion. 
In conclusion, we are glad to express our deep 
indebtedness to Mr. Joseph Schwartz for his en- 
thusiastic computational assistance in this work. 
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nuit pected as a consequence of the (3, 3) resonant interac- 


wave with j= 3. Although the configuration (pS 1»)32 












~— tion between the primary pion and the nucleon of a fits the decay angular correlations and the (E+, E_) 
gh j= . K-N bound state. This interpretation suggests plot rather well, it appears excluded by the production 
‘ naturally that the primary pion be emitted in the p angular distributions observed. 
1 by 
n. 
8 
PION- NUCLEON INTERACTIONS IN THE REGION OF THE HIGHER RESONANCES 
; P. Carruthers 
Laboratory of Atomic and Solid-State Physics, Cornell University, Ithaca, New York 
ch (Received April 17, 1961) 
| 
oa In this note we attempt to clarify certain fea- origin of the third resonance and the behavior of 
i tures of pion-nucleon interactions in the region reactions (1) and (2) are discussed. 4. The rel- 
of the “higher” resonances. We give most at- evance of the present remarks to the structure 
tention to the process of single-pion production.’ of the abrupt rise in the 7+ -p cross section at 
2 Data are most complete for the following re- 800-Mev pion lab kinetic energy*"" are ex- 
actions: plained. 
n +p—n7 +nt4n, (1) 1. The most damaging evidence against the 
explanation of the behavior of reaction (2) by 
1 +p—n-+n°+p. (2) OPE production [ Fig. 1(a)] is the fact’ that those 
ed. Recent experiments*~* near a Bev have re- events which contribute to the peak in the Go. 
3 ‘ spectrum do not occur primarily at the small 
vealed the following interesting features of re- = ; 
: values of the proton invariant momentum trans- 
actions (1) and (2). While the qualitative features f d. A? ted f th unoee 
of the usual 7-N “isobar” model’ are displayed pebnnnigiie < m9 paca — 2 
of Chew and Low” and Goebel.” Instead, the 
by reaction (1), reaction (2) behaves quite dif- a? Ginnie atthe on ok i iy flat 
ferently. For the latter process the number of pig camce Seattle adh ely sai “ig es 
: A over a wide range of A*, up to A*=35. More- 
3. protons with small recoil energy is in excess of 5 . 
10- over, ; of those events with A*<10 occur for 
the phase-space value. Moreover, at least at 


' 960 Mev, the distribution of events for (2) has a 
prominent peak at a 7-7 total center-of-mass 
energy squared, s, of about 22, in units of pion 
rest mass squared.*»*»® Neither of these features 
are displayed by reaction (1). 

The principal subjects to be discussed are as 
follows: 1. Present data for reaction (2) can- 
not be explained by the one-pion exchange (OPE) 
4 process | Fig. 1(a)]. 2. Consideration of the re- 
scattering by the nucleon of one of the pions re- 
: sulting from a m-7 collision [ Fig. 1(b)] cannot 
explain the behavior of reaction (2). 3. The 





e 


2 Re > eel ™ a eel 
— Dox 
(a) (b) 


FIG. 1. (a) Pion production by the one-pion ex- 
change process. (b) The rescattering of one of the 
pions in the 3-3 resonance. ; 





























s values less than 20, the peak occurring for 
20< s< 25. Thus the excess of events at low A? 
seems to have little to do with the peak in the 
Qan spectrum.* 

The second argument against the OPE explan- 
ation of the 960-Mev phenomena is that (for m7 
isospin ¢=1) the 7-7 effect should be stronger in 
reaction (1) than in reaction (2). As is well 
known, for Fig. 1(a) the ratio r=0(-+)/o(-0) is 
2/1 for an isospin 1 interaction. Although this 
ratio is not far from that observed’* near a Bev, 
OPE cannot possibly explain the 3-3 isobar be- 
havior in reaction (1). Even if one supposes that 
the over-all process is initiated by OPE, the 
succeeding events in reactions (1) and (2) are so 
different that the survival of the ratio 7 is not to 
be expected a priori. 

At this point it is convenient to introduce the 
isospin decomposition of the production ampli- 
tude.',*® Clearly there are four independent 
amplitudes corresponding to the reaction 71+N 
+2n+N. We specify the total isospin T and the 
m-m isospin ¢ to define the coupling scheme. The 
isospin amplitudes are written as A(T?t). The 


567 










VoLUME 6, NuMBER 10 


PHYSICAL REVIEW LETTERS 


May 15, 1961 





four amplitudes are A($0), A($1), A(31), and 
A(32). For even (odd) ¢ these amplitudes are 
symmetric (antisymmetric) under interchange 
of the momenta of the outgoing pions. The im- 
plications of this important symmetry property 
have been discussed in reference 16. Denoting 
the observed reactions by the charge of the out- 
going pions, we have’»’® 


A(-+) = -4V¥2.A(30) + 3[A(41) -A(31)]+ (1/3V 5) A (32), 
(3) 
A(-0) = -3V2[A (31) + $A ($1)] - (1/V 10) A (22). (4) 


In terms of these amplitudes the t= 1 OPE proc- 
ess gives A(31) = -2A(31), A(30) =A (32) =0. 
[This result is most simply derived by com- 
paring (3), (4), and the ratio 2/1 for o(-+)/o(-0) 
quoted above.| It was first suggested by Peierls”” 
that the data under discussion are suggestive of 
the relation A(31)~A(31). This allows the for- 
mation of a ¢=1 7-7 isobar in the final state of 
reaction (2), but not in reaction (1). This view, 
that the final state of reaction (2) is reached by 
some process more complicated than OPE, is 
certainly not in conflict with the A? distribution 
of the events in the 7-7 peak of reference 5. 

The possibility of a t=2 “di-pion” in the final 
state seems to be excluded by the predominance 
of the cross section for 1++p — 1++7°+p over 
1++p—2n++n in the energy range 900-1300 Mev. 
The experimental ratio’® is o(+0)/o(++)~=4. [In 
terms of the pure isospin cross sections,’® 
o(+0) = $0(31)+ 20(32); o(++) = 20(22).] 

The fact that in the crucial range of the 1-7 
energy the A? distribution of reaction (2) is 
qualitatively not at all like that of OPE clearly 
makes any extrapolation procedure unreliable at 
960 Mev. In our opinion the minimum require- 
ment for extrapolation is that the cross section 
in the physical region from which the extrapol- 
ation is made have the qualitative features of 
OPE. It is difficult to tell whether the data of 
Anderson et al.'® (at 890-Mev 1~ energy) satisfy 
this requirement. While the 1-7 cross section 
may well behave as indicated in reference 19, 
the procedure from which this result is derived 
is open to serious objection. 

We suggest that the following use of “extra- 
polation” data may be more informative. First, 
we consider the results of references 3-5 as 
good evidence for an isospin one m-7 resonance. 
The position and width are determined from 
these experiments. Using this information one 
can extrapolate from the pole in the unphysical 
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region to the data in the physical region (rather 
than the other way around). The difference be- 
tween this curve and that arising from OPE de- 
pends on the magnitude and phase of the non- 
OPE part of the production amplitude in a simple 
manner. The resulting information should be 
very useful in unraveling what seems to be a 
rather intricate situation. 

2. About a year ago a number of people*°~?? 
suggested that the 3-3 isobar state was reached 
predominantly through the 7-7 channel [ Fig. 1(b)]. 
“Quantitative” calculations have been made? 
taking into account essentially** the processes 
of Figs. 1(a) and 1(b), the nucleon being consid- 
ered infinitely heavy. Let ¢q(w,,w,) denote the 
appropriate partial-wave projection of the pro- 
duction amplitude, a denoting the relevant set 
of quantum numbers. w, and w, are the energies 
of the final pions; other variables are sup- 
pressed. Denoting the partial-wave projection 
of the OPE amplitude by ¢,q(w,,w,), the solu- 
tion for which pion 1 is in the 3-3 state is’ 


26 w0*” sind 
=e ef 
>, Wy Wo) e c0sd,d | (wy,wo)+ ap,) 
P .« 
x— f dw 
mJy 


5 is the 3-3 phase shift; v(p) is the usual cutoff 
function of the static-nucleon meson theory. It 
is easy to see the significance of the various 
terms in (5). Write e?5cosé6 =1+ ie? sind. The 
1 corresponds to the OPE term; ie*® sinddrg 
corresponds to 1-7 scattering followed by real 
excitation of the 3-3 isobar.** It is curious to 
note that these contributions interfere in a de- 
structive way. The integral term, proportional 
to the 3-3 scattering amplitude for meson 1, is 
due to off-the-shell scattering. This term re- 
ceives large contributions from large p; since 
it depends critically on the cutoff, its value is 
numerically unreliable. Numerical calculations 
for the D,, incident channel [leading to an s-wave 
and a p-wave (3-3 state) meson] indicate that the 
second term of (5) is about 3 times as large as 
the first. There seems no reason to doubt this 
feature for other channels, in view of the experi- 
mental evidence for the isobar model.*® 

When Eq. (5) is supplemented by the relevant 
charge state information, one finds the very in- 
teresting result’ that for a ¢=1 7-7 interaction, 
the 3-3 rescattering terms practically cancel 


pr*(p)4 (wv, wo) 


ww -w,) 





(5) 
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out for reaction (2), but not for (1). This can 
also be seen by considering the various ways in 
which Fig. 1(b) can contribute to the reaction (2) 
in a purely phenomenological way. Thus, were 

it not for the A? distribution at 960 Mev, the evi- 
dence would support the notions of a t=1 a-7 
resonance, the model of reference 1, and the 
extrapolation procedure for reaction (2). 

Neither the integral term nor the partial-wave 
projection of the OPE amplitude seems capable 
of producing such a sharp peak in the Q,, spec- 
trum as is found.*~> Moreover, the 7-N corre- 
lations introduced by symmetrization of (5) can- 
not reproduce this phenomenon. If we abandon 
the t=1 channels to some unknown process 
leading to a final-state 7-7 isobar, such that 
A(31)=A(1), then some interesting speculations 
can be made on the origin of the third resonance. 

3. Near the third (900-Mev) resonance one is 
naturally led to consider those final-state 27-N 
configurations in which both pions are in p states 
relative to the nucleon,”* one state having 7-N 
angular momentum 3/2. Of these states, those 
reached from the F,,. channel have the special 
property that the orbital momenta of both pions 
are parallel to the nucleon spin. As a conse- 
quence, the angular part of the production am- 
plitude for this channel is symmetric in the two 
pions, so that symmetrization affects only the 
energy part, Eq. (5). Consider reaction (1), and 
suppose the isospin 1 amplitudes in Eq. (3) cancel 
as discussed above. Then since the integral term 
in (5) is supposed to dominate, and since the in- 
tegral is a slowly varying function of w, and w, 
for moderate w, the 30 and $2 contributions to 
reaction (1) are approximately proportional to 


w,e°°" sind, /p,2+ we? sin6,/p,*, (6) 


for all angles. [For w,=w, the integral is the 
same for both terms in (6).] Thus the Bose sym- 
metry, which requires that either pion have a 
chance to scatter off the nucleon, leads to a 

large constructive interference for w, and w, 
both near the 3-3 resonance energy. From this 
argument one expects a peak in reaction (1) at 
the third resonance. Whether a peak occurs in 
(2) depends on the sensitivity of the unknown 
intermediate state to the initial energy. If this 
state is not sensitive, o(-0) should rise abruptly 
above the threshold for production of the di-pion 
(740 Mev for s = 22) and remain constant there- 
after. Below the third resonance both o(-+) and 
o(-0) will be rising as the resonance is ap- 








proached.”” [o(-+) is already large.] If o(-0) 
levels off and o(-+) decreases above the third 
resonance, then the ratio ry should drop. This 
picture is not incompatible with experiment; 
unfortunately, the ratio at 960 Mev is not yet 
certain.'® It should be emphasized that in the 
present picture only one pion interacts with the 
nucleon rather than both, as in the theories of 
Peierls,”® and Wong and Ross.”® 

4. It is of interest to ask whether the abrupt 
rise in the 1*-p cross section above 700 Mev 
has something to do with the 740-Mev threshold 
for producing an s = 22 “particle” in the final 
state. First we give a simple argument that 
such a process leads to a strong D,, state, as 
suggested by phenomenological analysis.”° Ac- 
cording to this picture, one is producing a spin- 
one, odd-parity di-pion. Shortly above thresh- 
old the orbital momentum of the di-pion relative 
to the nucleon is primarily s-wave. Thus we 
have J =1/2 or 3/2, odd parity, i.e., S,. and 
Dg, for the initial state.*° The other possibility 
giving D,, is the mechanism of the second res- 
onance: an s-wave pion accompanying the 3-3 
isobar.'s*” 

A measurement of the @ values for the various 
pairs of particles in the final state of the reac- 
tion 7++p ~ 7*+n°+p at about 800 Mev would be 
very useful. Should this reaction behave as 
does reaction (2), one could verify the p-wave 
character of the 7-7 resonance by measuring 
the angular distribution of the decay pions in 
their own center-of-mass system.*’ Let 6 be 
the decay angle relative to the beam direction. 
For the D,, case the distribution is 1+3cos*6. 
This does not interfere with the isotropic S,, 
component. A similar measurement could be 
made for reaction (2). Preliminary data'® at 
910, 1090, and 1260 Mev indicate the presence 
of the 3-3 7+-p state in the reaction n++p — nt+n°+p. 

In conclusion we wish to emphasize once more 
the importance of measuring the various charge- 
state ratios throughout the resonance region. 
The author would like to thank Professor H. A. 
Bethe and T. Truong for stimulating discussions. 
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I-SPACE PARITY CONNECTION IN NONLEPTONIC WEAK INTERACTIONS 


W. B. Zeleny 
School of Physics, University of Sydney, Sydney, N. S. W., Australia 
(Received April 3, 1961) 


In a recent paper by Pais,’ it was shown that 
the parity properties of nonleptonic hyperon de- 
cays suggest a doublet approximation for both 
strong and weak interactions and the existence 
of two isotopic spins, the doublet spin, f, and 
the K spin, K, with the usual isotopic spin, T, 
given by T-f+K. The nonleptonic weak inter- 
action Hamiltonian, H, is assumed to consist of 
two parts, 

H=H@+H®, (1) 
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where H™ is a scalar andH™ a vector in J space. 
H® andH™ are to be separately parity conserv- 
ing, but their sum is to be parity nonconserving. 
Thus, there exists a connection between rotations 
in J space and parity. It is the purpose of this 
note to suggest a possible basis for this connec- 
tion. 

Any connection between J space and parity im- 
mediately brings to mind the 7 meson, which is 
pseudoscalar and isovector. Clearly, the ad- 
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dition or removal of a 7-meson field to an inter- 
action will change both its parity and its Fspace 
transformation properties. In particular, if w° 
is a scalar and isoscalar neutral meson, then 
the interaction 


H =W(w°+°)+W'n’ +Herm. conj., (2) 


where W°* is either scalar or pseudoscalar and 
is a scalar inJ space, has precisely the desired 
properties of (1). This is accomplished without 
the artificiality of introducing y*® and 7 matrices 
ad hoc in some terms of H and not in others. 
These y® and 7 matrices appear instead in the 
strong interaction, where their positions are 
determined by the requirements of parity and 
isotopic spin conservation. The connection be- 
tween J space and parity in nonleptonic weak in- 
teractions would appear simply as a result of 
the J-space and parity properties of the w° and 

7 mesons. 

Now, experimental evidence for the existence 
of the isoscalar neutral meson is meager. Abash- 
ian, Booth, and Crowe? have found evidence 
for either a new neutral meson or a two-pion 
resonant state of total energy about 310 Mev in 
the reaction 


p+d— He®+w°*. 


Berkelman, Cortellessa, and Reale,° on the other 
hand, failed to observe an w® of this energy ina 
search for the reaction 


ytprp+w®. 


Sachs and Sakita* have argued against a spin- 
less w° of this energy on the grounds that it 
would be expected to yield too large a ratio for 
W(K* + w°+*)/W(K°+n*+77). However, if w° 
is scalar, they are comparing reaction of op- 
posite parities, and the decay K* + w°+n* might 
be expected to proceed at a rate comparable to 
that of direct 7 decay. These arguments, of 
course, would not rule out a spinless w° of mass 
greater than the K-7 mass difference. 

In any event, should the scalar w° exist, the 
following realization of (2) immediately suggests 
itself: 


H = gl ,y5N, + Ngy5N,)(w°+ 7°) 
- V2(N ,y°N, - Nzy®N,)n* ]+ Herm. conj., (3) 


where the notation is as in reference 1. From 


(1) and (3) we have 
H = g(N,y°N, + Nzy®N,)w° + H.c., 

H™ = g[(Niy°N, + Ngy®N,)n° 

- V2(N ,y°N, - N,y°N,)n* ]+ H.c. 


The Hamiltonian (3) has the following important 
properties: 

I. H has the same renormalizable Yukawa 
structure as the strong interactions. 

II. H™ contains both AS = 2AK, and AS = -2AK, 
terms. This is absolutely essential since H ™ 
must contribute to both 5“+n+a~ and 5+ +p+7°. 
H™ contains only AS = -2AK, terms, which is 
sufficient since it need only contribute to >* + 
n+at, Dt +p+n°, and Y°sp+n-. 

III. Since 7* decays have AS=2AK,, they may 
proceed only through H“™. On the other hand, 

a pseudoscalar interaction cannot contribute to 
tT decays. Thus H™ must be scalar, and the 
possibility (WH scalar, H™ pseudoscalar) is 
ruled out. Hence, 5++n+zn* and 5~+n+77 
are S-wave and P-wave, respectively, and 
should it turn out that the parities are the other 
way around, this model must be discarded. 

Iv. H™ is of the “allowed form”; that is, the 
AI, =0 part of H™ shares with the strong 7 inter- 
action invariance under 


N,+i€,T,N,, 1*> -en', 
N,~+1€,T.N,, 1°+-en°, 
N,+i€,T.N,, w°+ ew, 
N,~ 1€,TN3, 


where the e’s equal +1, and €€,€,=~-1. Hence, 
in the approximation where we neglect the N, -N, 
mass difference, and provided the K interactions 
do not violate this symmetry appreciably, H™ 
does not contribute to 5+ +n+7*,! and in this 
approximation, 5+ +n+7* is parity conserving, 
as required by experiment. 

The author is indebted to A. Pais for a pre- 
print of his work. 
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3k. Berkelman, G. Cortellessa, and A. Reale, Phys. 
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‘R. G. Sachs and 3. Sakita, Phys. Rev. Letters 6, 
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ERRATUM 





MICROWAVE EMISSION FROM AN OPTICALLY 
PUMPED ATOMIC SYSTEM. V. E. Derr, J. J. 
Gallagher, R. E. Johnson, and A. P. Sheppard 
[Phys. Rev. Letters 5, 316 (1960)]. 


The calculation giving an estimate of the line- 
width is incorrect, and the paragraph beginning 
at the bottom of page 316, “Second, an estimate 
of the linewidth may be made...,” should be de- 
leted from the article. On the basis of this ex- 
periment alone, no estimate of the linewidth can 
be made. 
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ABSTRACTS 








In this section are printed the abstracts of Articles 
that have been forwarded to The American Institute of 
Physics for publication in THE PHYSICAL REVIEW. 
In quoting information obtained from this section be- 
fore the appearance of the corresponding Article, ref- 
erence should be made to “Physical Review (to be 
published)” rather than to this Journal. 


DIFFUSION OF NEON, HT, AND DEUTERIUM 

IN LIQUID HYDROGEN. G. Cini-Castagnoli,* 
Istituto di Fisica dell’ Universita di Roma, Italy 
and Istituto Nazionale di Fisica Nucleare, Sezione 
di Roma, Italy, A. Giardini-Guidoni, Laboratorio 
di Chimica delle Radiazioni e Chimica Nucleare 
del Comitato Nazionale Energia Nucleare, Roma, 
Italy, and F. P. Ricci, Comitato Nazionale En- 
ergia Nucleare, Roma, Italy (Received February 
10, 1961; revised manuscript received March 

22, 1961). 


Experimental results on diffusion of Ne, HT, 
and D, in liquid H, are reported and analyzed 
with the corresponding-states principle. The 
features of the quantum deviations are indicated. 


*Now at: c/o Consiglio Nazionale delle Ricerche, 
Istituto di Metrologia-Sezione Termometrica, Torino, 
Italy. 


COMPARISON OF TWO THEORETICAL AP- 
PROACHES TO ELECTRON BEHAVIOR IN 
Ar-CO,, Ar-N,, Ar-H,, AND Ar-CO GAS 
MIXTURES. Martin A. Uman, The University 
of Arizona, Tucson, Arizona (Received March 7, 
1961). 


The electron drift velocity and electron average 
energy for low-energy electrons in binary gas 
mixtures of Ar-CO,, Ar-N,, Ar-H,, and Ar-CO 
are determined using two theoretical methods of 
approach: (1) a “distribution function” or “Boltz- 
mann equation” approach; and (2) an “average 
electron” approach. The results of the two theo- 
retical methods of approach are compared and 
discussed. 


OPTICAL ABSORPTION AND FLUORESCENCE 
OF OXYGEN IN ALKALI HALIDE CRYSTALS. 

J. Rolfe, F. R. Lipsett, and W. J. King,* Radio 
and Electrical Engineering Division, National 





Research Council, Ottawa, Canada (Received 
March 10, 1961; revised manuscript received 
April 17, 1961). 


Single crystals of NaCl, KCl, and KBr were 
grown from the melt in an oxygen atmosphere, 
and their optical absorption, fluorescence exci- 
tation, and fluorescence emission spectra were 
measured at 300°K, 77°K, and 4.2°K. The weak 
absorption band caused by oxygen was the same 
in all three cyrstals, and did not vary with tem- 
perature. The band had maximum absorption at 
5.0 ev and a half-width of 1.0 ev. All fluores- 
cence excitation spectra contained a component 
identical to this absorption band. The fluores- 
cence emission spectra consisted of a series of 
peaks in the wavelength range 4000-10000 A, 
with an approximately equal energy separation of 
1000 cm~. At 300°K, 12 to 15 peaks were re- 
solved, and at 4.2°K each of these peaks split 
into 4 to 6 components. From these optical re- 
sults and from paramagnetic resonance experi- 
ments, it is concluded that O,~ molecule ions 
located in anion sites in the crystal are respon- 
sible for the absorption and fluorescence. 


*Present address: Goodrich High Voltage Astro- 
nautics, Inc., Burlington, Massachusetts. 


CURRENT-CARRIER TRANSPORT WITH SPACE 
CHARGE IN SEMICONDUCTORS. W. van Roos- 
broeck, Bell Telephone Laboratories, Murray 
Hill, New Jersey (Received March 9, 1961). 


Differential equations are given for a general 
formulation that includes space charge. Arbitrary 
dependences of diffusivities and magnitudes of 
drift velocities on electrostatic field are con- 
sidered, and extension is made for applied mag- 
netic field. Though excess electron and hole con- 
centrations are not equal, the small-signal re- 
combination rate depends on a single lifetime, 
the “diffusion-length lifetime,” 7). The formu- 
lation is applied to one-dimensional drift with 
recombination for an injected pulse of electron- 
hole pairs. The exact electron and hole distri- 
butions are obtained in closed form for the linear 
small-signal case. The condition is given for 
linearity or substantially unperturbed applied 
field, E,. There are two principal types of so- 
lutions, essentially according to whether 7, is 
larger or smaller than the dielectric relaxation 
time, Tg. For T9> 7g, the electron and hole dis- 
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tributions in not too strongly extrinsic material 
are ultimately similar Gaussian distributions 
displaced by the “polarization distance,” xp, the 
distance electrons and holes drift apart in time 
(r°" -79 *)"*. These distributions drift at a 
velocity that differs from the ambipolar velocity 
by an amount which, besides being small for 
small Tq/ To» vanishes for equal mobilities. 

They spread, exhibiting an apparent diffusion. A 
“pseudodiffusivity,” D,,, is defined. For To> Tq 
and constant mobilities, D, is proportional to 
TqE 9 /%', with og the conductivity. The ambi- 
polar diffusivity and D, are additive. They are 
equal in intrinsic material for E, equal to kT/e 
divided by the Debye length (kTe/8mn;e”)”, or 10 
volt/cm for silicon at 300°K. An extension to a 
nonlinear case involving high-level injection is 
given; concentration-dependent D,, and velocity 
function are defined. For sufficiently strongly 
extrinsic material and T9> 7,, the minority car- 
riers drift in a delta pulse that leads the majority 
carriers distributed in an exponential tail of char- 
acteristic length x p, which may be quite large. 
For nonconstant mobilities and To> Ty ambipolar 
velocity in the majority-carrier or “reverse” di- 
rection may occur. For Tqg> To» the other princi- 
pal type of solution gives distributions that in gen- 
eral (and for constant mobilities) drift in the re- 
verse direction. Involving also regions of local 
carrier depletion and thus generation as well as 
recombination, they may persist for times long 
compared with 7,, being attenuated then with 

time constant 7 . 


SPECIFIC HEATS OF DELTA-PHASE Zr-H AND 
Zr-D. Walter J. Tomasch, Atomics International 
Division of North American Aviation, Canoga 
Park, California (Received March 6, 1961). 


The specific heats of fcc ZrH,,,, and ZrD, .. 
have been measured in the temperature interval 
30-500°C. The data are interpreted in terms of 
a harmonic oscillator model for the hydrogen and 
deuterium specific heat contributions. This model 
is in accord with recent inelastic neutron scatter - 
ing studies and predicts an isotopic depression of 
the deuteride Einstein temperature by a factor of 
1/V2 relative to the hydride value. Over the in- 
terval 30-200°C, the data and model are quanti- 
tatively consistent. At higher temperatures, the 
deuteride specific heat is somewhat smaller than 
anticipated. Quantities of noncubic y-phase ma- 
terial are known to be present in the samples 
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used, particularly in the deuteride, and this is 
thought to be a likely cause for the deviations. 
The hydride Einstein temperature calculated 
from the difference between deuteride and hydride 
specific heats at 150°C is 1500+300°K, as com- 
pared with the inelastic neutron scattering value 
of 1500+ 60°K. 


GYROMAGNETIC RATIO OF NICKEL FERRITE. 
G. G. Scott, Research Laboratories, General 
Motors Corporation, Warren, Michigan (Re- 
ceived March 10, 1961). 


The gyromagnetic ratio of the ferrite NiOFe,0, 
was determined by measurements of the Ein- 
stein— de Haas effect. The g’ value of 1.849 
+ 0.002 indicates that the magnetization is large- 
ly due to the Ni** ions as in the Néel model. 
Comparison with values of g determined by fer- 
romagnetic resonance investigations furnishes 
evidence as to the validity of the Kittel— Van 
Vieck relation for ferrites. 


ANTIFERROMAGNETIC RESONANCE IN FeF, 
AT FAR-INFRARED FREQUENCIES. R. C. 
Ohlmann”* and M. Tinkham, Department of Phys- 
ics, University of California, Berkeley, Califor- 
nia (Received March 14, 1961). 


Antiferromagnetic resonance in single crystals 
of FeF, was observed between 1.5°K and 66°K 
(Ty =78.4°K) in the far-infrared region. The reso- 
nance frequency at 7 ~=0 was found to be 7(0) = 52.7 
+0.2 cm™*. Using the antiferromagnetic resonance 
relation derived by Kittel, Nagamiya, Keffer, and 
others, and using the experimental value for the 
static susceptibility, the uniaxial anisotropy con- 
stant at absolute zero, K(0), was inferred to be 
1.1x10° ergs/cm® (40 cm~!/atom). Measurements 
of the splitting of the line caused by an external 
magnetic field gave g) =2.25+0.05. The aniso- 
tropy energy in FeF,, being primarily due to the 
crystalline field- spin-orbit interaction, may be 
described by a term )/;DS, ,? in the Hamiltonian. 
By including this interaction in the molecular- 
field treatment, we have determined D=-9+2 
cm™*, and have calculated the temperature de- 
pendences of the sublattice magnetization, the 
anisotropy constant, the resonance frequency, 
and the linewidth. The last two were compared 
with the experimental results and found to be in 
reasonable agreement. The linewidths were 
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found to follow a 7* law above 15°K. The search 
over a frequency region of 13-70 cm“! for the 
absorption lines expected in the paramagnetic 
region was unsuccessful, possibly indicating a 
relaxation time of less than 10-™” sec. 


*Now at Westinghouse Research Laboratories, Pitts- 
burgh, Pennsylvania. 


SIZE EFFECTS IN THIN SUPERCONDUCTING 
INDIUM FILMS. A. M. Toxen, International 
Business Machines Research Center, Yorktown, 
New York (Received March 8, 1961; revised 
manuscript received April 27, 1961). 


Measurements of residual resistivity, super- 
conductive critical temperature, and critical 
magnetic field have been carried out on indium 
films ranging in thickness from 650 A to 126000 
A. The thickest films had the bulk critical field 
and critical temperature. The variation of re- 
sidual resistivity with thickness is consistent 
with Fuchs’ model if one assumes an intrinsic 
resistivity p, of 1.3110~° ohm cm and an in- 
trinsic mean free path J, of 152000 A. The value 
of p,/, so obtained was 2.0x10™ ohm cm’. The 
critical temperature was found to be a systematic 
function of film thickness, increasing with de- 
creasing thickness. The magnitude of this change 
in critical temperature is in good agreement with 
a simple model relating critical temperature to 
elastic stresses in the films. The penetration 
depth, as calculated from the critical field by 
means of the London theory or the Ginzburg- 
Landau theory, was found to increase with de- 
creasing film thickness. This result is consistent 
with a nonlocal model and implies a coherence: 
length of approximately 2600 A. 


GALVANOMAGNETIC EFFECTS IN SEMICON- 
DUCTORS AT HIGH ELECTRIC FIELDS. E. M. 
Conwell, General Telephone and Electronics 
Laboratories, Incorporated, Bayside, New York 
(Received March 13, 1961). 


A treatment of magnetoconductivity is developed 
for high electric fields and general energy-band 
structure using a partial solution of the Boltz- 
mann equation in a form similar to that set up 
by McClure for low electric fields. The present 
treatment is valid when the scattering processes 
are such that the distribution function varies but 





a small amount over an entire constant-energy 
surface, or, in the case of the many-valley band 
structure, over the part of a constant-energy 
surface within each valley. In the latter case 
different distribution functions must be used for 
the different valleys. The elements of the mag- 
netoconductivity matrix that results are expressed 
in terms of carrier concentration, total or within 
each valley, and averages over the carriers ofa 
quantity involving the momentum relaxation time 
and the § tensor defined by McClure. This tensor, 
which depends on the shape of the constant-energy 
surfaces and on the magnetic field strength, is 
evaluated for the individual valleys in a nonde- 
generate many-valley semiconductor. The mag- 
netoconductivity matrix is then in a form conven- 
ient for calculation of conductivity and galvomag- 
netic effects for either low or high electric fields. 
It is used to obtain expressions for anisotropy 
voltage and Hall coefficient in high electric fields 
involving the number of carriers in each valley, 
orientation of the valleys, and valley averages 
over quantities involving relaxation time and 
energy. 


DIFFERENTIAL PARAMAGNETIC EFFECT IN 
SUPERCONDUCTORS. Robert A. Hein and Ray- 
mond L. Falge, Jr., U. S. Naval Research Lab- 
oratory, Washington, D. C. (Received February 
2, 1961; revised manuscript received April 20, 
1961). 


The magnetic moment and the differential mag- 
netic susceptibility of two spherical samples of 
tin and tantalum have been measured as a func- 
tion of magnetic field and temperature. The dif- 
ferential paramagnetic effect (DPE) is observed 
in both ac and dc mutual inductance measurements 
provided the sample exhibits a good Meissner 
effect. For a superconducting sample in which 
the infinite conductivity behavior dominates the 
Meissner effect, the DPE will not appear in the 
ac measurements but will, under certain con- 
ditions, show up in dc measurements. The results 
of the dc mutual inductance measurements are 
used to classify the DPE as reproducible or non- 
reproducible. The former is characteristic of 
ideal Meissner-type superconductors while the 
latter is more characteristic of superconductors 
whose macroscopic magnetic properties are dom- 
inated by the classicai infinite electrical conduc- 
tivity behavior. The superconducting to normal 
transitions obtained by the three techniques are 
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compared and the data discussed in view of (a) oc- 
currence of the DPE; (b) magnitude of the DPE; 
(c) ability of these data to give information about 
the volume participation of the sample; and (d) re- 
lation of the “characteristic magnetic fields” de- 
termined by these transitions and the bulk critical 
field of the sample. 


STUDY OF THE SUPERCONDUCTING TRANSI- 
TION TEMPERATURE IN DILUTE THALLIUM 
SOLID SOLUTIONS. D. J. Quinn and J. I. Bud- 

’ nick,* International Business Machines Research 
Laboratory, Poughkeepsie, New York (Received 
March 6, 1961). 


The superconducting transition temperature, 
T;, has been measured for dilute solid solutions 
of indium, bismuth, and lead in thallium. The 
transition temperature is found to increase in 
all cases, thus exhibiting a behavior opposite to 
that observed by Serin, Lynton, and co-workers 
in their studies of solid solutions of tin, indium, 
and aluminum. Some loss in residual impurity 
scattering which occurs upon annealing suggests 
the migration of solute atoms to grain boundaries 
in the dilute alloys. 


*Now at International Business Machines Watson 
Laboratory, Pougkeepsie, New York. 


INTERSTITIAL DEFECTS IN LiH AND NaCl IR- 
RADIATED AT LOW TEMPERATURES. F. E. 
Pretzel and R. L. Petty, Los Alamos Scientific 
Laboratory, Los Alamos, New Mexico (Received 
March 10, 1961). 


Expansions of the lattice parameters of both 
Li and NaCl exposed to 6 radiation at 77 and 
245°K have been observed using a vacuum x-ray 
diffractometer stage operated at 77°K. The 
source of the 8 radiation was 5.56% LiT con- 
tained in a LiH crystal which was ground and 
mixed with NaCl in the samples. The LiH lattice 
parameter increased 0.0020+ 0.0007 A in three 
samples. This result is consistent with density 
measurements and with the formation of inter- 
stitial He*® atoms in the LiH lattice. The NaCl 
lattice parameter increased 0.0061+ 0.0009 A 
even though the average 8 dose for the NaCl was 
less than 1% of that in the Li(H, T) and tritium 
8 particles with maximum energy cannot produce 
direct displacements in NaCl. A mechanism is 
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proposed for the production of Cl~ interstitials 
and vacancies (Frenkel defects) in NaCl which 
is similar to one proposed by Klick for the pro- 
duction of F andH centers. This mechanism is 
related to observations of the production of a 
centers in alkali halides irradiated at low tem- 
peratures. 


KAPITZA RESISTANCE BETWEEN HELIUM 
AND METALS IN THE NORMAL AND SUPER- 
CONDUCTING STATES. W. A. Little, Stanford 
University, Stanford, California (Received 
February 16, 1961). 


The contribution of the conduction electrons 
of a metal to the heat flow across a helium- 
metal interface has been calculated. It is found 
that the “totally” reflected phonons from the 
fluid play an important role in the transfer mech- 
anism as had been predicted previously. The 
dominant term in the heat flow is proportional 
to T°AT in agreement with the experimentally 
observed value on lead. However, the numerical 
agreement is poor. The reasons for this are 
discussed. In the superconducting state it is 


shown why this heat transfer becomes inoperative. 


Several interesting consequences of this calcula- 
tion are given. In particular, it is predicted 
that a phonon- drag effect may be observed be- 
tween the conduction electrons in the metal and 
the phonons in the fluid and vice versa. Also, 

it is shown that the variation of the Kapitza 
resistance with applied magnetic field can help 
to distinguish between the various contributions 
to the heat flow. 


GENERALIZED KOOPMANS’ THEOREM. J. C. 
Phillips, Department of Physics and Institute 
for the Study of Metals, University of Chicago, 
Chicago, Illinois (Received February 24, 1961). 


Koopmans’ theorem states that if the wave 
function of a many-electron system is approxi- 
mated by a Slater determinant of Hartree-Fock 
one-electron wave functions, with one-electron 
energies defined as the difference in energy of 
(N +1)- and N-particle systems, then these one- 
electron energies are given by the expectation 
value of the Hartree-Fock Hamiltonian with re- 
spect to the one-electron wave functions. Koop- 
mans’ theorem is here generalized to include 
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correlation effects by using Hubbard’s expres- 
sion for the total energy of a free electron gas. 
The resulting one-electron Hamiltonian contains, 
in first-order, screened exchange. Hubbard’s 
lowest polarization diagram gives, in addition, 
part of the screened second-order Coulomb 
interactions, which is small for metallic 
densities. Collective terms are also obtained. 
Comparison with the Bohm-Pines Hamiltonian 
shows a one-to-one correspondence but with dif- 
ferent cutoff functions in each term. Following 
Hubbard we extend the method to include the 
effects of a periodic potential to first order. The 
resulting one-electron Hamiltonian provides a 
convenient and accurate basis for self-consistent 
energy band calculations including exchange and 
correlation in metals and semiconductors. 


INDIRECT EXCHANGE MODEL FOR FERROMAG- 
NETIC METALS. S.H. Liu, Research Center, 
International Business Machines Corporation, 
Yorktown Heights, New York (Received March 

13, 1961). 


The ferromagnetic properties of rare earth 
metals and their alloys are discussed in terms 
of the indirect exchange model. It is shown by 
the molecular field approximation that in calcu- 
lating the Curie temperatures of these metals 
the simple theory of Fréhlich and Nabarro and 
Zener is applicable. The second-order energy 
terms calculated by Ruderman and Kittel, Kasuya, 
and Yosida are important in discussing the low- 
temperature properties. Some numerical results 
are obtained which are in good agreement with 
the experiments. 


CONTRIBUTION OF LATTICE SCATTERING 
BETWEEN NONEQUIVALENT VALLEYS TO 
FREE-CARRIER INFRARED ABSORPTION IN 
SEMICONDUCTORS. H. Risken, Deutsche Philips 
G.m.b.H., Zentrallaboratorium Aachen, Germany, 
and H. J. G. Meyer, Philips Research Laborato- 
ries, Eindhoven, Netherlands (Received March 
13, 1961). 


In a multivalley band-structure like that of the 
conduction band of germanium a contribution to 
the absorption of infrared radiation by free car- 
riers is made by a scattering process which does 
not play a role in normal transport processes but 





which may be of importance for hot-electron phe- 
nomena; this is the scattering between nonequiva- 
lent valleys of the conduction band. The absorp- 
tion induced by this extra scattering process re- 
sults in a transfer of the electrons from the (111) 
valleys to the (100) valleys or the [000] valley. 

A quantum mechanical calculation was made of 
the partial absorption constant u;* due to this 
scattering on the basis of a deformation- potential 
type theory. The final formula obtained for u;* 

is similar to that derived previously for the par- 
tial absorption constant °Pt due to optical intra- 
valley scattering. The physical significance of 
some limiting forms of ».;* at low temperatures 
is discussed. 


ULTRAVIOLET ABSORPTION SPECTRA IN RUBY. 
A. Linz, Jr., and R. E. Newnham, Massachusetts 
Institute of Technology, Cambridge, Massachu- 
setts (Received March 6, 1961). 


The optical properties of highly doped rubies 
have been investigated. Al,O, and Cr,O, forma 
complete solid-solution series; single crystals 
containing up to 5 mole % chromia were grown 
by the Verneuil technique. The optical absorp- 
tion of a group of ultraviolet crystal-field bands 
near 3400 A was studied as a function of tempera- 
ture, crystal orientation, and chemical composi- 
tion. The intensity of these absorption bands 
varies with the square of the Cr concentration, 
perhaps indicating strong chromium-pair inter- 
actions. 


LOW-TEMPERATURE SPECIFIC HEAT OF IN- 
DIUM AND TIN. C. A. Bryant* and P. H. Keesom,T 
Department of Physics, Purdue University, La- 
fayette, Indiana (Received March 9, 1961). 


The heat capacities of indium and tin were 
measured between 0.4 and 4.2°K. In the normal 
state, the specific heat could be represented by 
AT~?+yT+aT*+ T*+ pT". For Sn, in molar mil- 
lijoule units, A=0; y, the coefficient in the elec- 
tronic term, is 1.80; a =0.242, corresponding to 
a Debye temperature, 6,, of 200°K; g=0.004; and 
u=0.00014. For In, A, the coefficient of a nucle- 
ar electric quadrupole term, is calculated to be 
8.97x10~* from resonance data; y=1.61 for one 
ingot and 1.59 for another; 9,=109 and 108°K; 
and g=0.008. In the superconducting state, the 
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specific heat of Sn could be expressed as the 
normal lattice term plus an electronic term of 
the form ayT,, exp(-bT,,/T), with T, =3.70°K 
(0.02 deg lower than found in a magnetic meas- 
urement), a=7.63, and b=1.41 when 2< T,/T< 7; 
the value of b agrees with infrared measurements 
of the energy gap. This sort of analysis could 
not be applied to In, for below 0.8°K, the total 
superconducting specific heat was less than the 
normal lattice term. A possible interpretation 
is that 6, is 9% higher in the superconducting 
state than in the normal metal at 0.4°K; this is 
not supported, however, by the recent acoustic 
measurements of the elastic constants by Chan- 
drasekhar and Rayne. The anomaly is not as yet 
understood, but a few plausible explanations are 
discussed. 

*Present address: International Business Machines 
Corporation, Yorktown Heights, New York. 

tOn leave of absence at the Kamerlingh Onnes Labo- 
ratory, Leiden University, Leiden, The Netherlands. 


THEORY OF RESONANCE ABSORPTION OF EN- 
ERGY BY A ROTATING SOLID. Joseph Dreitlein, 
University of Pennsylvania, Philadelphia, Penn- 
sylvania, and Horst Kessemeier, Washington Uni- 
versity, Saint Louis, Missouri (Received Feb- 
ruary 13, 1961). 


The theory of both nuclear resonance energy 
absorption and the relaxation of a coherent mag- 
netization transverse to an applied external mag- 
netic field is developed for nuclei in a mechani- 
cally rotating solid. The consequences of the for- 
malism in the simple case of a solid composed of 
effectively isolated nuclear pairs are presented. 
For the more general nuclear lattice, a procedure 
for isolating and experimentally measuring the “ex- 
change” interaction between the nuclei in the solid 
is proposed. Finally, the moments of the energy 
absorption line shape for the rotating solid are 
investigated and the second and fourth moments 
are explicitly calculated. 


PINCH EFFECT IN INDIUM ANTIMONIDE. A. G. 
Chynoweth and A. A. Murray, Bell Telephone 
Laboratories, Murray Hill, New Jersey (Received 
March 10, 1961). 


The critical current at which pinching occurs in 
indium antimonide has been measured by three 
independent methods: (i) by noting the current at 
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which the pinched current-voltage characteristic 
deviates from the unpinched characteristic that 
is obtained in the presence of a longitudinal mag- 
netic field, H, using crystals of sufficiently high 
resistance for the avalanche breakdown current 
to be considerable before pinching sets in as the 
electric field is increased; (ii) by noting the cur- 
rent at which the magnetoresistance as a function 
of H shows a change in its behavior; and (iii) from 
a study of the critical current as a function of H. 
The three methods lead to a value for the critical 
pinching current of 4 to 5amp. This current is 
the same for both single-crystal and polycrystal- 
line samples and is insensitive to small changes 
in the donor concentration or cross-sectional 
area of the crystal. The value of the critical cur- 
rent leads to a mean carrier temperature of 0.04 
ev in avalanche breakdown. An irregular form of 
noise is observed when the crystal is operated in 
the transition region between the pinched and un- 
pinched conditions, and it is thought that this 
noise is caused by pinching-unpinching instabili- 
ties. 


SINGLE-CRYSTAL AND POLYCRYSTAL RE- 
SISTIVITY RELATIONSHIPS FOR YTTRIUM. 

J. K. Alstad, R. V. Colvin,* and S. Legvold, In- 
stitute for Atomic Research and Department of 
Physics, Iowa State University, Ames, Iowa 
(Received February 20, 1961). 


Different proposals for calculating polycrystal 
resistivities from single-crystal values are ap- 
plied to yttrium metal. It is shown that a simple 
average, yielding ppoly=3(2p, +p > gives the 
best fit to experimental data. 


*Presently at the Edgar C. Bain Laboratory of the 
United States Steel Corporation, Monroeville, Penn- 
sylvania. 


LOW-TEMPERATURE DISSIPATION PEAK IN 
NIOBIUM. P. G. Bordoni, M. Nuovo, and 
L. Verdini, Istituto Nazionale di Ultracustica, 
Rome, Italy (Received February 7, 1961). 


The energy dissipation coefficient Q~' and the 
resonant frequency of a circular plate of niobium 
have been measured as a function of temperature 
in the range 60-300°K. The measurements have 
been carried out at four different frequencies 
from 18 kc sec”* and 174 kc sec™* with strain 
amplitudes smaller than 10°’. For each vibration 
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mode, a pronounced peak is found for the dissipa- 
tion coefficient while the frequency-temperature 
curves show a corresponding inflection. The 
temperature T,, of the dissipation peak depends 
on the vibration frequency according to an Ar- 
rhenius equation which characterizes the ther- 
mally activated relaxation effects. 

The activation energy W and the limiting re- 
laxation time 7, of the process have been com- 
puted and these values [W =0.265 ev and (7,)"* 
=61x10" sec™*] are in agreement with the values 
previously found in some fcc metals for the re- 
laxation effect due to the motion of dislocations. 

The value of W found in niobium shows that 
the effect cannot be directly produced by the 
motion of interstitial atoms of hydrogen. 


DECAY OF PHOSPHORESCENCE FROM A DIS- 
TRIBUTION OF TRAPPING LEVELS. W. L. 
Medlin, Field Research Laboratory, Socony 
Mobil Oil Company, Incorporated, Dallas, Texas 
(Received March 13, 1961). 


In a previous paper it was shown that the usual 
model for second-order decay predicted the 
correct form for the decay in many thermolu- 
minescent crystals but gave the wrong behavior 
for the parameters involved. Specifically, it 
was shown that b and m in the decay expression, 
I=1,[b/(b+t)!”, should behave differently as 
functions of the decay temperature and the de- 
gree of trap filling than is observed experimen- 
tally at temperatures near or below the glow 
peak. In the present paper it is shown that the 
discrepancies can be accounted for by assuming 
a first-order decay from a distribution of trap- 
ping levels. Most of the results are based ona 
Gaussian distribution but it is shown that other 
distributions can produce similar results. The 
first-order mechanism is justified by consider- 
ing the relative magnitudes of the rate constants 
for trap emptying, retrapping, and recombina- 
tion. At temperatures well above the glow peak 
this assumption is no longer justified, but in 
this range the second-order decay predicts the 
observed results for b and m. The effects of 
retrapping and of crystal dimensions are con- 
sidered. Also, the effect on the glow peak of 
having a distribution of levels rather than a set 
of discrete levels is worked out and it is shown 
that the peak is broadened appreciably even for 
relatively narrow distributions. 










STATISTICAL ELECTRON DENSITY DISTRIBU - 
TIONS AND THOMAS-FERMI-DIRAC SCREENING 
FUNCTIONS FOR NEUTRAL ATOMS. Adolf A. 
Abrahamson, The City College of New York, New 
York, New York and Brookhaven National Labo- 
ratory, Upton, New York (Received March 3, 
1961). 


Statistical electron density distributions and 
Thomas-Fermi-Dirac (TFD) screening functions 
have been obtained for the 104 elements corre- 
sponding to integral atomic numbers Z =2 to 
Z=105, for 117 values of radial distance from 
the atomic center in each case. These results 
were calculated, in part, by using Thomas’ so- 
lutions of the TFD equation in terms of the sta- 
tistical electron distributions for nonintegral 
values of Z and Jensen’s boundary conditions. 
Values for argon and copper are given here. The 
complete set of tables for all 104 elements has 
been deposited with the American Documentation 
Institute Auxiliary Publications Project. 


TEMPERATURE VARIATION OF IONIC MOBILI- 
TIES IN HYDROGEN. Lorne M. Chanin, Honey- 
well Research Center, Hopkins, Minnesota (Re- 
ceived March 6, 1961). 


Measurements have been made of the tempera- 
ture variation of the mobilities of positive ions in 
hydrogen over the range 77-300°K. The zero-field 
mobility values are 1,=12.3 cm?/volt-sec (300°K), 
13.3 (195°K), and 13.0 (77°K). The present results 
at 300°K are in agreement with the data of Lauer, 
Bradbury, and Mitchell at low £/p, and at high 
E/p, they agree with Rose’s measurements. Only 
a single ion species was observed in the present 
studies. Reasons are given which support the be- 
lief that the ion observed in these measurements 
was H,*. 


ANALYTIC HARTREE-FOCK WAVE FUNCTIONS 
FOR THE 3p SHELL ATOMS. R. E. Watson, Avco 
Research and Advanced Development Division, 
Wilmington, Massachusetts, and A. J. Freeman, 
Materials Research Laboratory, Ordnance Mate- 
rials Research Office, Watertown, Massachusetts 
(Received March 2, 1961). 


Hartree-Fock wave functions have been ob- 
tained for the 3p row atoms, i.e., for neutral 
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Al, Si, P, S, Cl, and Ar, and for Cl~. Solutions 
were determined in analytic form using a version 
of Nesbet’s symmetry and equivalence restric- 
tions to simplify the calculations for atoms with 
both closed and unclosed shells of the same / 
value. These restrictions, the reason for their 
use, and their relation to other open-shell meth- 
ods are discussed and the calculated one-electron 
wave functions and their eigenvalues are pre- 
sented. 


ATOMIC BEAM RESONANCE EXPERIMENTS 
WITH STORED BEAMS. H. Mark Goldenberg, 
Daniel Kleppner, and Norman F. Ramsey, Har- 
vard University, Cambridge, Massachusetts (Re- 
ceived March 13, 1961). 


The atomic beam separated oscillatory field 
resonance technique has been used to study the 
hyperfine frequency of cesium, which is per- 
turbed by collisions with storage box walls. With 
a wall coating of long, straight-chain saturated 
hydrocarbons, resonances are observed after as 
many as 200 wall collisions. A theory of the ef- 
fect of wall collisions on the hyperfine frequency, 
which is in qualitative agreement with experi- 
mental results, is described. The shape of the 
resonance curve is analyzed by a detailed con- 
sideration of the statistical nature of the wall 
collision. 


CESIUM TRANSITION PROBABILITIES FOR OP- 
TICAL PUMPING. W. Bruce Hawkins, Depart- 
ment of Physics, Oberlin College, Oberlin, Ohio 
(Received March 6, 1961). 


Optical transition probabilities are presented 
between the several magnetic sublevels of the 
Cs'** ground state via the first excited state, 
given separately for the two fine-structure com- 
ponents. Graphs of the populations of represent- 
ative states are presented as a function of the 
time the atoms are illuminated, and a table is 
given of the populations of atoms absorbing inte- 
gral numbers of photons from both D lines simul- 
taneously, both polarized and unpolarized. 


ANGULAR DISTRIBUTION OF THE GROUND- 
STATE NEUTRONS FROM THE C*5( p, n)N** AND 
N**(p,n)O*® REACTIONS. C. Wong, J. D. Ander- 
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son, S. D. Bloom, J. W. McClure, and B. D. 
Walker, Lawrence Radiation Laboratory, Univer- 
sity of California, Livermore, California (Re- 
ceived March 1, 1961). 


The angular distribution of the neutrons from 
the C*5( p,”)N** and N**( p, 2)O** ground-state re- 
actions has been measured in 10° steps from 0° 
to 150° for incident proton (laboratory) energies 
of approximately 6.5, 6.9, 7.5, 8.1, 8.6, 8.9, 9.4, 
10.6, 11.4, and 12.2 Mev. Additional measure- 
ments were made for C** at 5.0, 10.2, 10.9, and 
13.3 Mev, and for N’* at 5.5, 7.7, 7.8, and 13.6 
Mev. A calibrated plastic or stilbene scintillator 
was used in order to obtain absolute differential 
cross sections. Time-of-flight techniques on the 
Livermore variable-energy cyclotron allowed 


positive identification of the ground-state neutrons. 


The targets [CO,(58 % C**) and N,(90% N**)] were 
sufficiently thick to average out the effects of 
possible compound-nucleus contributions. From 
preliminary fits to the C** angular distributions, 
Glendenning and Bloom have inferred an effective 
neutron-proton interaction inside the nucleus. 


MULTILEVEL ANALYSIS OF THE TOTAL NEU- 
TRON CROSS SECTION OF Pu? BELOW 12 ev. 
O. D. Simpson and M. S. Moore, Atomic Energy 
Division, Phillips Petroleum Company, Idaho 
Falls, Idaho (Received March 7, 1961). 


Measurements of the total cross section of Pu™ 
have been analyzed with the use of a multilevel 
formula, under the assumption that the observed 
resonance asymmetries are due to interference 
in a small number of fission channels. As was 
the case for U***, the analysis is interpreted as 
evidence that there are differences in the aver- 
age fission widths of resonances belonging to the 
two possible spin states of the compound nucleus. 


ANGULAR DISTRIBUTION AND RANGES OF N”* 
PARTICLES FROM N“* ON N™. K.S. Toth, Oak 
Ridge National Laboratory, Oak Ridge, Tennessee 
(Received March 2, 1961). 


The distribution of N’* particles from the neu- 
tron transfer reaction N'*(N"*, N**)N"* was inves- 
tigated from 3.5°-32° in the laboratory. Range 
curves for N'* particles were obtained from 16° to 
32° (laboratory system). Transfers leaving both 
residual nuclei in their ground states were dis- 
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tinguished from those in which the products are 
left in excited states. It is found that (1) the 
ground-state transfer cross section decreases as 
the bombarding energy is lowered; (2) of the N*® 
excited states only the first and/or second con- 
tribute significantly to the transfer cross section; 
(3) for a given bombarding energy the excited 
state transfer distribution peaks at an angle 
larger than that at which the ground-state trans- 
fer distribution reaches its maximum; (4) from 
the ground-state transfer distribution an 7, of 

2.2 f can be determined, while the excited state 
transfer distribution yields an r, of about 1.65 f; 
and (5) the excited state distributions are more 
consistent with the tunneling mechanism of Breit 
and Ebel than are the ground-state distributions. 


BETA-DECAY COUPLING CONSTANT AND THE 
ft VALUE OF O*. R. J. Blin-Stoyle and J. Le 
Tourneux, Clarendon Laboratory, Oxford, Eng- 
land (Received March 1, 1961). 


A calculation is made of the effect of a charge- 
dependent internucleon potential on the O“(0+, 
T =1)~N**(0*, T=1) g-decay matrix element. It 
is found that a not unreasonable strength and 
form for such a potential can lead to a reduction 
in the matrix element = 1% which is of the right 
order of magnitude to resolve the present dis- 
crepancy between the @- and p-decay polar vec- 
tor coupling constants. 


PROTON SPECTRA FROM THE NITROGEN BOM- 
BARDMENT OF FLUORINE. C. E. Hunting, Oak 
Ridge National Laboratory, Oak Ridge, Tennes- 
see (Received March 13, 1961). 


Energy spectra and absolute differential cross 
sections of protons from bombardment of F’*® 
with 21.4-Mev and 27.4-Mev N** were measured. 
The energy spectra were analyzed in terms of the 
statistical model with a level density exp(E*/T), 
where E* is the excitation energy in the residual 
nucleus of the assumed reaction F'*(N"*, p)P*. 
The nuclear temperature T increased from about 
1.9 Mev at 0 deg to about 2.5 Mev at 145 deg 
¢.m., but did not vary with bombarding energy. 
The angular distributions show minima near 90 
deg c.m., with anisotropy increasing as the 
bombarding energy and proton energy increase. 
At the lower bombarding energy there is ap- 





proximate agreement with the Ericson-Strutin- 
ski theory of angular momentum effects in com- 
pound nucleus processes, over the entire ob- 
served ranges of angle and proton energy. The 
fit is consistent with o=3 in the nuclear level 
spin distribution (2j+1) exp[-j(j+1)/(207)]. At 
the higher bombarding energy a considerable 
excess of protons are emitted into the backward 
hemisphere, especially at high proton energies, 
suggesting a direct-interaction mechanism in 
which the observed proton comes from the F’®. 


STUDY OF THE DIFFERENTIAL CROSS SEC- 
TIONS OF DEUTERON STRIPPING REACTIONS 
AS A FUNCTION OF THE INCIDENT ENERGY. 
E. W. Hamburger,* University of Pittsburgh, 
Pittsburgh, Pennsylvania and Universidade de 
Sao Paulo, Sao Paulo, Brasil (Received March l, 
1961). 


Angular distributions for the reactions C*(d, p) 
and O'*(d,) to the ground and first excited states 
of C** and O” have been obtained at deuteron en- 
ergies of 10.2, 12.4, and 14.8 Mev. These re- 
sults and those of previous experiments at other 
energies between 3 and 19 Mev are compared with 
the predictions of the usual plane-wave Born ap- 
proximation theory, due originally to Butler. The 
assumption of plane waves leads to the prediction 
that the differential cross section is a function of 
energy and scattering angle only through the trans- 
fer momentum g: Angular distributions at dif- 
ferent energies should coincide when plotted vs q. 
It is found that this is only approximately true on 
the main stripping peak of the angular distribution; 
at larger angles the cross section is a more com- 
plicated function of energy and angle. Further- 
more the stripping peak itself shifts as a function 
of g, principally in the deuteron energy range 10 
to 19 Mev. 

*Permanent address: Department de Fisica, Facul- 


dade de Filosofia, Universidade de Sio Paulo, Sao Paulo, 
Brasil. 


BETA-GAMMA DIRECTIONAL CORRELATION 
IN Eu'™. K. S. R. Sastry, R. F. Petry, and R. G. 
Wilkinson, Indiana University, Bloomington, In- 
diana (Received March 9, 1961). 


The energy dependence of the beta-gamma an- 
gular correlation between the outer beta-ray 
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group of Eu'™, W,=1.86 Mev, and the 123-kev 
gamma ray of the daughter Gd'™ has been meas- 
ured with a shaped magnetic field beta-gamma 
coincidence spectrometer. A negative correla- 
tion coefficient, accurate to about 5%, is obtained 
which ranges from -0.10 to -0.17 in the energy 
region 0.80 to 1.60 Mev. It is shown that the mod- 
ified Bij approximation must be relaxed to explain 
the data. The nuclear parameters (Kotani’s nota- 
tion) which result are: x =-0.24+0.05, u=+0.05 
+0.03, Y =+0.76+0.08. These values are com- 
pared with those which have been reported for 
Eu'™: x=u=0, Y =0.69+ 0.06. 


LIFETIMES OF ENERGY LEVELS IN Al™, Mn*®, 
Cu™*, Rh’™, AND I’* EXCITED BY SLOW-NEU- 
TRON CAPTURE. S. J. du Toit* and L. M. Bol- 
linger, Argonne National Laboratory, Argonne, 
Illinois (Received March 8, 1961). 


Measurements of delayed coincidences between 
y rays were used to determine the half-lives of 
states excited by slow-neutron capture in various 
nuclei. The half-lives found were: Mn, first 
excited level 11.472 nsec, second excited level 
5.1+0.5 nsec; Rh’, first excited level 2.64 0.2 
nsec, 96-kev transition <0.6 nsec; Cu™, first 
excited level <0.3 nsec, second excited level 
<0.3 nsec; Al”, first excited level 2.340.2 
nsec; I’, 30-kev level 8.8+1.0 nsec, 137-kev 
level 8.0+ 0.6 nsec, 90-kev transition < 0.7 nsec. 
Measurements of pulse-height spectra satisfying 
the condition of prompt and delayed coincidence 
and measurements at neutron resonances of io- 
dine were used to establish the presence of ex- 
cited states at 30 and 137 kev in ?*. A previ- 
ously unreported neutron resonance was observed 
in I’?’ at a neutron energy of 10.7 ev. 


*Permanent address: South African Atomic Energy 
Board, Pretoria, Union of South Africa. 


DIRECTIONAL CORRELATION MEASUREMENTS 
IN Hf'”*. U. Bertelsen, J. Borggren, and O. Na- 
than, Institute for Theoretical Physics, Univer- 
sity of Copenhagen, Copenhagen, Denmark (Re- 
ceived March 10, 1961). 










Measurements of the directional correlations 
of the 1335-93 kev, 1345-93 kev, and 1390-93 key 
y-y cascades in Hf'”* have been made, using a li- 
quid source of 21-day W’” equilibrium with 9.3- 
min Ta!”*. It is shown that the data support the 
tentative spin assignments of Gallagher et al. 
The 1390-kev y ray is found to be almost pure 
quadrupole radiation. If the spin and parity of the 
1430-kev level is 1+, as suggested by Gallagher 
et al., the multipolarity of the 1335-kev y ray is 
determined by our measurements to be mainly 
Mi. Using these results, the amount of E0 radi- 
ation in the 1390-kev transition has been es- 
timated. It is found that approximately 70% of 
the 1390-kev K-conversion electrons are due to 
monopole transitions. 


DETERMINATION OF M1-E2 MIXING AMPLI- 
TUDES IN Mg”’, Al”’, Si?®, AND P™. G. J. Mc- 
Callum,* Atomic Energy of Canada Limited, 
Chalk River, Ontario, Canada (Received March 
13, 1961). 





A method of analyzing gamma-ray angular dis- 
tribution and direction-polarization data to yield 
values for the relative amplitudes of mixed mul- 
tipole gamma-ray transitions is described. The 
analysis of the data does not demand a knowledge 
of the mechanism of formation of the gamma- 
emitting level. Measurements are reported on 
some mixed M1-E2 de-excitation gamma rays 
from low-lying levels in Mg**, Al*’, Si*®, and 
Pp"; these levels were excited by inelastic pro- 
ton scattering. The following results were ob- 
tained: 











Initial state Final state 
Energy Spin & Energy Spin & E2 amplitude 
Nucleus (Mev) parity (Mev) parity M1 amplitude 
Mg™ 0.98 3/2+ ground state 5/2+ 0.30+0.15 
0.98 3/2+ 0.58 Mev 1/2+ 0.15 +0.05 
1.01 3/2+ ground state 5/2+ -0.32+0.14 
1.28 3/2+ ground state 1/2+ 0.21+0.03 
1.27 3/2+ ground state 1/2+ -0.25 +0.15 














ama = oo. & of owe Ff tt mH Fh HCO ee ee ee eee ee ee ee 





- 


le 

















VoLuME 6, NUMBER 10 





PHYSICAL REVIEW LETTERS 


May 15, 1961 





CLUSTER MODEL CALCULATION ON THE 
ENERGY LEVELS OF THE LITHIUM ISOTOPES. 
y. C. Tang and K. Wildermuth, Florida State 
University, Tallahassee, Florida, and L. D. 
Pearlstein, John Jay Hopkins Laboratory of 
Pure and Applied Physics, General Atomic, San 
Diego, California (Received February 24, 1961). 


A variational procedure is adopted to determine 
the energies of the levels of the highest space- 
symmetry type in Li® and Li’. The trial wave 
functions employed take into consideration the 
existence of cluster structures in those nuclei. 
With a simple two-body force, it is shown that 
the computed energies of the various states are 
in reasonable agreement with experiment. The 
*F y. level in Li’, as yet undetermined experi- 
mentally, is found to have an excitation energy 
of about 5.6 Mev and a rather large level width. 
The calculation also indicates that to explain the 
splitting of levels in those nuclei, a constant two- 
body spin-orbit force of the pure neutral form is 
inadequate. 


YIELD RATIOS FOR THE ISOMERIC PAIR Sc**””,# 
FORMED IN (a,an) AND (a,n) REACTIONS. 

Sylvia M. Bailey,* Lawrence Radiation Laborato- 
ry, University of California, Berkeley, Califor- 
nia (Received January 23, 1961). 


The yield ratio of Sc*”” /Sc** was measured in 
Sc“(a,an)Sc* reactions with helium ions of ener- 
gies between 20 and 43 Mev and at 320 Mev. The 
measured ratio was nearly constant at a value of 
1.5 between 20 and 43 Mev and was 0.62 at 320 
Mev. The Sc**”/Sc* ratio was measured in 
K*(a,n)Sc* reactions at 10 and 43 Mev with val- 
ues of 0.3 and 0.9, respectively. The isomer ra- 
tio was calculated for the K*'(10-Mev a,n)Sc“ 
reaction by means of a compound-nucleus model 
and for the Sc**(320-Mev a, an)Sc*“ reaction by 
means of a classical knock-on model. The cal- 
culated ratios were 0.32 and 0.51, respectively. 


*Present address: National Bureau of Standards, 
Washington 25, D. C. 


PION PRODUCTION IN ELECTRON- POSITRON 
COLLISIONS. Frank Chilton, Department of 
Physics, University of Washington, Seattle, 
Washington (Received March 13, 1961). 


The properties due to the presence of sym- 





metries in pion production processes in electron- 
positron collisions are discussed. Cross sections 
are calculated for the production reactions e* +e~ 
~+7°+y, et+e” ~et+e-+7°, andet+e” ater”. 
The photon spectrum for the reaction e++e~ + 

n+ +n~+y is also calculated. The role of form 
factors and some of the possible effects of reso- 
nant strong interactions are discussed. 


SPALLATION OF URANIUM AND THORIUM 
NUCLEI WITH Bev-ENERGY PROTONS. B. D. 
Pate* and A. M. Poskanzer, Chemistry Depart- 
ment, Brookhaven National Laboratory, Upton, 
New York (Received March 6, 1961). 


Cross sections have been measured for the 
production of isotopes of uranium, protactinium, 
thorium, and actinium in the irradiation of U*™, 
U**>, and Th*” with 0.68- and 1.8-Bev protons. 

In addition, some yields were determined at other 
bombarding energies ranging up to 6.2 Bev. Cal- 
culations of the cross sections at 1.8-Bev bom- 
barding energy were made, based on recent Monte 
Carlo calculations of the knock-on phase of the 
interaction, combined with published systematics 
of nuclear evaporation, and several assumptions 
as to fission-evaporation competition. Even with- 
out fission competition the calculated yields are 
considerably lower than the experimental ones, 
indicating a failure of the model for the knock-on 
phase of the reaction to predict sufficient proba- 
bility for simple processes with low deposition 
energy. In view of the better agreement with ex- 
periment that previous workers have obtained 
with this type of calculation for 0.34-Bev protons 
on uranium, it is suggested that the present dis- 
crepancy may be due to an overemphasis of meson 
processes in the proliferation of the knock-on 
cascade in the calculation. 

*Present address: Chemistry Department, Washing- 
ton University, St. Louis, Missouri. 


INTERACTION OF 1.0-, 2.0-, AND 3.0-Bev 
PROTONS WITH Ag AND Br IN NUCLEAR EMUL- 
SION. Elizabeth W. Baker and Seymour Katcoff, 
Chemistry Department, Brookhaven National Lab- 
oratory, Upton, New York (Received March 10, 
1961). 


Stars produced in insensitive nuclear emulsions 
by 1.0-3.0 Bev protons have been classified into 
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different groups depending on whether light frag- 
ments and/or fission fragments are emitted. 
Alpha-particle spectra and angular distributions 
are presented for each of the various groups. 
The probability for light-fragment emission in- 
creases rapidly with increasing beam energy up 
to 2.0 Bev. The angular distribution of the light 
fragments is peaked forward but also shows a 
preference for emission at 90° to the beam. Fis- 
sion events increase from ~3 % of the interactions 
with Ag and Br at 1.0 Bev to ~11% at 3.0 Bev. 
Ranges and angular distributions are also given 
for the recoil and fission fragments. 


ANGULAR ASYMMETRY THEOREMS FOR DE- 
CAY PRODUCTS. Murray Peshkin, Argonne 
National Laboratory, Argonne, Illinois (Received 
February 1, 1961). 


The method of Eberhard and Good, for using 
decay angular distributions to determine the spins 
of unstable systems, is developed. Their inequal- 
ity, which expresses the impossibility of getting 
too symmetric a decay pattern from an asym- 
metric initial state, is tightened. The results 
are generalized to include spinning decay prod- 
ucts. The Adair analysis is generalized in the 
same framework. 


SEARCH FOR POSITIVE PARTICLES OF MAS- 
SES ABOUT 500mze AND 1400meg. V. Cook, 

D. Keefe, L. T. Kerth, P. G. Murphy, W. A. 
Wenzel, and T. F. Zipf, Lawrence Radiation 
Laboratory, University of California, Berkeley, 
California (Received March 2, 1961). 


The mass spectrum of positive particles in a 
secondary beam from the Bevatron has been 
measured at a distance of 90 feet from the inter- 
nal target. At a confidence level of 93 % the pro- 
portion of particles with mass between 420m, 
and 630m, in the 1-Bev/c beam is estimated to 
be <6x10°*. An upper limit of about 3 x10~° is 
obtained for the proportion of particles of mass 
approximately 1400m, in the 2.3-Bev/c beam. 


TEST OF GLOBAL SYMMETRY IN PION-BARYON 
INTERACTIONS BY K~+p REACTIONS. Jogesh C. 
Pati,* Department of Physics, University of Mary- 
land, College Park, Maryland and Norman Bridge 
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Laboratory of Physics, California Institute of 
Technology, Pasadena, California (Received 
March 8, 1961). 


Under the hypothesis that the K-meson interac- 
tions do not mask the symmetries of the pion- 
baryon interactions appreciably, the branching 
ratios of the K~ +p reactions are studied to test 
the validity of global symmetry. The T~-matrix 
formalism of Matthews and Salam is adopted to 
calculate the branching ratios. The new Dalitz- 
Tuan solutions for KN scattering lengths, which 
incorporate the (K*, K°) mass difference and the 
new branching ratios of the various K~ +p reac- 
tions, presented at Kiev, are adopted in the anal- 
ysis. The errors in the experimental branching 
ratios are so chosen as to satisfy the Amati- 
Vitale inequality. It is found that the a~ and 5* 
(also a+, though poorly) Dalitz-Tuan solutions 
can explain the branching ratios for K~ captured 
at rest. The extension of the analysis to 30-Mev 
incident K~ mesons under the zero-range approx- 
imation leads to very poor agreement with exper- 
iments. 


*Now at the Norman Bridge Laboratory of Physics, 


California Institute of Technology, Pasadena, California. 


SINGULARITIES OF SCATTERING AMPLITUDES 
ON UNPHYSICAL SHEETS AND THEIR INTER- 
PRETATION. R. Blankenbecler, M. L. Gold- 
berger, S. W. MacDowell,* and S. B. Treiman, 
Palmer Physical Laboratory, Princeton Univer- 
sity, Princeton, New Jersey (Received March 1, 
1961). 


The analytic structure of two-particle scatter- 
ing amplitudes on the unphysical sheet of the Rie- 
mann surface reached by crossing the two-par- 
ticle cut is discussed. The singularities of the 
amplitudes there are shown to be poles, and their 
physical interpretation is studied. The way in 
which bound states appear on the physical sheet 
in the Mandelstam representation, both as iso- 
lated poles and as cuts, is traced in detail. The 
properties of partial-wave amplitudes and of the 
field amplitude as a function of energy and angle 
and of energy and momentum transfer are dis- 
cussed. Finally, a few remarks are made in 
connection with unstable states. 


*Present address: Centro Brasileiro de Pesquisas 
Fisicas, av. Venceslau Braz. 71, Estado da Guanabara, 


Brasil. 
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EXCITATION SPECTRUM IN MANY-BOSON 
SYSTEMS. Fumihiko Takano,* Metcalf Chemical 
Laboratories, Brown University, Providence, 
Rhode Island (Received March 14, 1961). 


When only two-body correlations are fully taken 
into account, there appears an energy gap in the 
excitation spectrum for many-boson systems as 
shown by Girardeau and Arnowitt and as con- 
firmed by Wentzel. This energy gap is shown to 
disappear and the spectrum to become phonon- 
like again and proportional to the momentum for 
small momentum, if we construct the eigenmodes 
of excitations (collective excitations) taking into 
consideration appropriate higher order terms. 


*On leave of absence from Tokyo University of Ed- 
ucation, Tokyo, Japan. 


OPTIMAL MEASURING APPARATUS. Mutsuo M. 
Yanase, Palmer Physical Laboratory, Princeton 
University, Princeton, New Jersey (Received 
March 6, 1961). 


An upper limit for the accuracy of the meas- 
urement of a simple quantity, which does not 
commute with a conserved quantity, is obtained 
in terms of the “size” of the apparatus. The 
“size” of the apparatus is defined as the mean 
square value, #*M?, of the conserved quantity for 
the apparatus which is, in the example chosen, 
the z component of the angular momentum. The 
measured quantity is the projection of a spin in 
a perpendicular direction. It is found that the 
probability of an unsuccessful measurement is 
at least 1/8M?. 


THEORY OF SINGLE PION PRODUCTION IN 
PROTON- PROTON COLLISIONS. J. lizuka”™ and 
A. Klein, University of Pennsylvania, Philadel- 
phia, Pennsylvania (Received March 9, 1961). 


Single pion production in nucleon-nucleon col- 
lisions in the Bev energy region is studied by 
means of the picture due to Weizsacker-Williams 
(peripheral interaction model) and in the (3/2, 3/2) 
resonance approximation for the final-state pion- 








nucleon interaction. The excitation curve for the 
reaction p+p—p+n+zn*, nucleon momentum dis- 
tributions, and Q-value distributions in the center- 
of-mass system are calculated and compared with 
the published experimental data as well as with 

the work of Lindenbaum and Sternheimer. Agree- 
ment is reasonable with both, thus indicating that 
these particular distributions are insensitive to 
the isobar production mechanism (single pion ex- 
change), though more recent experiments seem 

to have indicated clearly a preference for the 
theory developed here. 


*Present address: Enrico Fermi Institute for Nu- 
clear Studies, University of Chicago, Chicago, Illinois. 


DOMAINS OF DEFINITION FOR FEYNMAN 
INTEGRALS OVER FEYNMAN PARAMETERS. 
Tai Tsun Wu, The Institute for Advanced Study, 
Princeton, New Jersey (Received March 1, 1961). 


Given a Feynman diagram, the corresponding 
integral over real Feynman parameters is mean- 
ingful and analytic in a certain domain in the 
space of the Lorentz invariants formed from the 
external momenta, each of which is on the mass 
shell. In the case where all the masses are equal, 
the intersection of these domains for all proper 
convergent diagrams is studied. For the cases 
of four and five external lines, the real inter- 
sections are explicitly found; for the case of six, 
seven, and eight external lines, procedures for 
finding the real intersections are given. A knowl- 
edge of the real intersection makes it possible 
to construct geometrically a subset of the com- 
plex intersection. Generalization to unequal ex- 
ternal masses is briefly considered. 


PROPERTIES OF NORMAL THRESHOLDS IN 
PERTURBATION THEORY. Tai Tsun Wu, In- 
stitute for Advanced Study, Princeton, New 
Jersey (Received March 1, 1961). 


Making use of the relation between a Feynman 
diagram and the corresponding electric circuit, 
several properties of the normal thresholds are 
established. 






